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Abstract 

HQ  Air  Force  Operational  Test  and  Evaluation  Center 
(AFOTEC)  is  responsible  for  the  operational  test  and 
evaluation  of  major  weapon  systems  during  their  acquisition 
phase.  In  operational  testing  at  AFOTEC,  it  is  often  the 
case  that  not  enough  data  or  samples  are  available  to  make 
high  confidence  classical  inferences  on  various  statistics. 
The  use  of  Bayesian  statistical  methods  may  be  a  solution  in 
which  higher  confidence  inferences  can  be  made  on  statistics 
by  using  other  sources  of  data. 

This  research  explored  the  applicability  of  Bayesian 
methods  to  the  problem  of  determining  the  reliability  of  a 
series  system,  a  parallel  system,  and  a  bridge  system.  The 
time  to  failure  of  each  component  in  each  system  is  assumed 
to  be  exponentially  distributed.  A  basic  exponential  case 
and  a  basic  binomial  case  are  also  evaluated. 

In  general,  using  simulated  data  and  real  data  when 
possible,  the  Bayesian  methods  produced  confidence  intervals 
that  were  much  tighter  than  the  classical  inference  methods, 
thus  allowing  a  decision  maker  to  have  higher  confidence  in 
making  a  decision.  System  level  data  and  aggregated 
component  level  data  produced  a  much  tighter  confidence 
interval  than  component  level  data  aggregated  to  system. 
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USING  BAYESIAN  STATISTICS  IN  OPERATIONAL  TESTING 

I .  Introduction 

General  Issue 

In  operational  testing  at  the  Air  Force  Operational 
Test  and  Evaluation  Center  (AFOTEC) ,  it  is  often  the  case 
that  not  enough  data  or  samples  are  available  to  make  high 
confidence  classical  inferences  on  various  statistics.  The 
use  of  Bayesian  statistical  methods  may  be  a  solution  in 
which  higher  confidence  inferences  can  be  made  on  statistics 
by  using  other  sources  of  data. 

The  testing  of  new  weapons  systems  in  the  Air  Force  is 
divided  into  two  types  of  testing:  Developmental  Test  and 
Evaluation  (DT&E)  and  Operational  Test  and  Evaluation 
(OT&E) .  In  developmental  testing,  the  contractor 

participates  directly  with  the  DT&E  test  agency  and  the  goal 
is  to  produce  a  workable  system.  However  in  the  operational 
testing  phase,  AFOTEC  performs  independent  tests  to  ensure 
the  system  meets  Air  Force  standards  and  requirements  and  is 
ready  for  production.  Traditionally  these  two  types  of 
testing  were  separate,  but  with  the  new  call  for  streamlined 
acquisition  and  reduction  of  testing  funds,  AFOTEC  has  to 
look  at  incorporating  DT&E  data  into  OT&E. 

At  present,  AFOTEC  will  use  DT&E  data  only  when  DT&E 
data  was  collected  by  sources  that  are  independent  of,  and 
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not  influenced  by,  the  contractor.  Additionally,  AFOTEC  has 
to  make  sure  the  data  from  the  DT&E  system  was  from  a  system 
configured  like  the  system  to  be  used  in  OT&E. 

Bayesian  statistics  incorporates  "prior  information", 
such  as  engineering  estimates  and  DT&E  data,  into  making 
statistical  inferences.  For  example,  let  f(y|U)  be  the 
probability  distribution  function  (pdf)  for  the  time  to 
failure  for  a  component,  where  y  is  the  data  and  U  is  a 
parameter  of  the  pdf.  To  utilize-  prior  information, 
Bayesian  statistics  assume  that  the  parameter  U  is  indeed  a 
random  variable  in  itself,  with  its  own  pdf.  Bayesian 
statistics  call  the  distribution  of  the  random  variable  U  an 
a  priori  or  a  prior  pdf.  This  prior  pdf  is  called  g{U), 
which  expresses  the  state  of  knowledge  (or  ignorance)  about 
U  before  any  new  sample  data  is  analyzed.  There  are  various 
methodologies  to  convert  prior  information,  such  as 
engineering  estimates  and  DT&E  data,  into  a  prior 
distribution.  Then  given  the  prior  distribution  g(U),  the 
probability  distribution  f(y|U),  and  the  new  data  y,  Bayes' 
theorem  can  be  used  to  calculate  the  so-called  posterior  pdf 
g(U|y) .  The  posterior  distribution  g(U|y)  is  in  a  sense, 
the  updated  prior  distribution.  The  posterior  distribution 
can  now  be  used  to  update  the  pdf,  f(ylU),  the  time  to 
failure  distribution. 
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Problem  Statement 


Can  a  system  reliability  analysis  based  on  Bayesian 
statistics  using  DT&E  data  be  used  to  reduce  the  number  of 
samples  and  test  time  needed  for  OT&E  and  provide  the  same 
or  better  confidence  levels  than  classical  statistical 
inference  methods? 

Research  Objectives 

The  objectives  of  this  research  are:  1)  to  investigate 
different  Bayesian  methods  and  various  weighting  schemes  for 
using  prior  information;  2)  develop  a  methodology  for  using 
Bayesian  statistics  with  DT&E  data  to  predict  system 
reliability  based  on  data  obtained  from  a  program  which  is 
in  the  process  of  transition  from  development  to  production 
and  compare  the  results  with  those  obtained  from  classical 
statistical  methods. 

Scope 

This  research  will  use  DT&E  and  OT&E  reliability  data 
obtained  from  a  test  program  when  possible.  Simulated 
reliability  data  will  be  used  to  verify  and  validate  the 
methodology.  This  research  will  concentrate  mainly  on  the 
exponential  case.  This  research  is  limited  to  exploring  the 
applicability  of  Bayesian  statistical  methods  to  the  problem 
of  determining  the  reliability  of  a  series  system,  a 
parallel  system,  and  a  bridge  system.  The  time  to  failure 
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of  each  component  in  each  system  is  assumed  to  be 
exponentially  distributed.  A  basic  exponential  case  and  a 
basic  binomial  case  are  also  evaluated.  The  methodology 
developed  in  this  research  is  based  on  the  assumption  that 
both  DT&E  and  OT&E  data  are  available.  The  next  chapter 
will  identify  and  discuss  some  of  the  different  approaches 
of  Bayesian  statistics  as  well  as  comparing  and  contrasting 
these  with  the  classical  inference  method. 
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II .  Literature  Review 


Overvi ew 

This  chapter  will  begin  with  the  classical  statistical 
inference  and  Bayesian  statistical  inference  approaches. 
The  chapter  concludes  with  current  thinking  on  the  use  of 
Bayesian  statistical  concepts  for  the  type  of  reliability 
analysis  found  in  the  literature. 

Two  Views  of  Probability- 

In  the  classical  method,  the  probability  of  an  event  is 
defined  as  the  relative  frequency  of  an  event  occurring 
during  repeated  runs  of  an  experiment .  Let  S  be  the  set  of 
all  possible  outcomes  of  an  experiment.  Let  A  be  one  of  the 
possible  outcomes.  The  classical  method  defines  the 
probability  of  A  (the  likelihood  of  outcome  A  happening)  as 
P[A]  =  limn-«o  (a/n) ,  where  a  is  the  number  of  times  outcome 
A  occurs  during  n  runs  of  the  experiment.  The  value  of  P[A] 
is  an  unknown  constant  that  can  be  estimated  if  enough  runs 
of  the  experiment  are  conducted.  This  classical  view  is 
also  known  as  the  "frequentist"  viewpoint  (14:50). 

The  Bayesian  method  defines  probability  subjectively. 
In  contrast  to  the  classical  method,  the  Bayesian  method 
does  not  consider  P[A]  to  be  an  objective  estimate  of  the 
likelihood  of  outcome  A.  Instead  the  Bayesian  method 
considers  P[A]  as  a  rational  expression  of  one  individual's 
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degree  of  belief  that  outcome  A  will  occur.  Different 
values  of  P[A]  based  on  different  perceptions  of  prior 
experience  can  be  assigned  by  different  individuals.  This 
does  not  mean  the  Bayesian  method  abandons  experimentation. 
Experimentation  provides  an  important  tool  for  an  individual 
to  refine  their  assessment  of  P[A]  (14:50-51). 

Consider  the  following  experiment:  a  coin  is  flipped, 
with  the  outcome  to  be  "head"  or  "tail"  (S  =  {head,  tail })  . 
Let  the  event  A  =  {head}.  To  find  P[A],  the  classical 
method  would  require  n  flips  of  the  coin  (n  large),  count 
the  m  occurrences  of  A  =  {head},  and  compute  P[A]  to  be 
approximately  m/n.  In  contrast,  the  Bayesian  method  might 
declare  that  P[A]  =  0.5  based  on  experience  with  similar 
coins  that  a  "head"  is  as  likely  to  occur  as  a  "tail"  in  any 
flip.  To  validate  and  refine  the  estimate,  a  number  of 
experimental  flips  of  the  coin  would  then  be  made. 

There  are  distinctive  differences  between  classical 
statistical  inference  and  Bayesian  statistical  inference 
methods.  The  classical  statistical  inference  uses  inductive 
reasoning  while  the  Bayesian  method  uses  deductive 
reasoning.  The  Bayesian  method  takes  into  account  prior 
information  while  the  classical  method  only  considers  the 
prior  information  in  an  informal  manner,  if  at  all.  The 
classical  method  is  more  restrictive  due  to  the  exclusive 
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use  of  sample  data.  To  achieve  the  same  level  of 
inferences,  the  Bayesian  method  usually  requires  less  sample 
data  as  compared  to  the  classical  method.  When  sample  data 
is  expensive  or  difficult  to  obtain,  as  reliability  data 
often  is,  the  Bayesian  method  has  the  advantage  of  being 
able  to  use  prior  information.  Table  1  summarizes  the 
characteristics  of  the  two  approaches. 


Table  1-  A  Summary  of  Certain  Characteristics  of  the 
Sampling  Theory  and  Bayesian  Methods  of  Statistical 
Inference  (6:169) 

^  II  mil  'I  I  . . . . . .  .  — —  . nil  mil  ■ii-i.ii  1. 


Characteristic 

Sampling  Theory 

Bayesian  j 

Parameter (s)  of 
Interest 

Unknown 
constant (s) 

Random  Variable (s) 

Prior  Distribution 

Does  not  exist 

Exists  and 
explicitly  assumed 

Sampling  Model 

Assumed 

Assumed 

Posterior 

Distribution 

Does  not  exist 

Explicitly  derived 

Method  of 

Reasoning 

Inductive 

Deductive 

Type  of  Interval 
Estimate 

Confidence 

interval 

Probability 

interval 

Role  of  Past 
Experience 

Not  applicable 

Applicable 

Purpose  of 

Sampling 

Experience 

Supply  data  for 
making  inferences 

Confirm  or  deny 
expected 
performance  as 
predicted  from 
past  experience 

Quality  of 
Inferences 

More  restricted 
than  Bayes' 
because  of 
exclusive  use  of 
sample  data 

Depends  on  ability 
to  quantitatively 
relate  past 
experience  to  the 
sample  data 

Quantity  of  Sample 
Data 

Bayes'  approach 
usually  requires 
less  because  it 
utilizes  relevant 
past  data 
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Martz  and  Waller  identify  four  advantages  for  using  the 
Bayesian  method:  1)  increased  quality  of  the  inferences, 
provided  the  prior  information  accurately  reflects  the  true 
variation  in  the  parameters;  2)  reduction  in  testing 
requirements  (test  time  and/or  sample  size);  3)  inferences 
that  are  unacceptable  must  come  from  incorrect  assumptions 
and  not  from  inadequacies  of  the  method  used  to  provide  the 
inferences;  and  4)  the  rules  for  manipulating  probability 
statements  on  components  into  corresponding  statements  on 
system  reliability  are  well  known,  whereas  equivalent  rules 
for  manipulation  of  confidence  statements  are  not  (6:172-3). 
Classical  Statistical  Inference. 

Let  the  population  be  defined  as  the  total  number  of 
systems  in  use  for  a  particular  system.  When  it  is  not 
practical  to  test  all  members  of  a  population,  a  subset  or 
sample  of  systems  can  be  selected  from  the  population  for 
testing.  The  sample  is  presumed  to  be  randomly  selected 
such  that  each  system  having  an  equal  chance  of  being  in  the 
sample . 

Two  types  of  inferences  can  be  drawn  from  the  sample. 
A  point  estimate  approximates  the  true  value  of  the 
population  parameter  with  a  single  value.  However,  this 
single  value  of  a  point  estimate  does  not  provide  any 
information  about  the  uncertainty  associated  with  that 
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estimate.  To  account  for  uncertainty,  an  interval  estimate 
provides  an  estimate  that  the  sample  statistic  is  within  a 
certain  interval  about  the  true  population  parameter. 
Consider  the  mean  of  a  sample,  ~x,  which  is  a  point  estimate 
of  the  population  mean.  Compared  with  the  true  population 

mean  n,  ~x  will  differ  by  some  unknown  sampling  error  amount 
8.  Thus  ir  =  |.i  +  8or‘x=|x-8.  For  large  sample  sizes, 
both  the  sampling  distribution  and  sampling  error 
distribution  are  approximately  normal  by  the  Central  Limit 
Theorem.  From  this,  a  confidence  interval  based  on  the 
population  standard  error,  ex,  is  developed.  The  population 
standard  error  is  the  deviations  or  distances  of  the 
population  measurements  from  their  mean. 

fx  -  z*<t)  <  fl  <  fx  +  Z*(T)  (1) 
where  z  is  the  critical  value  associated  with  a  specified 
probability.  The  population  standard  error  is  usually  not 
known  in  most  cases.  However,  it  can  be  estimated  by  using 
the  sample  standard  error,  s  with 

s  =  (i)  tiXi-Hy  (2) 

and  T  =  (i)  fxi  .  (3) 
Thus  an  approximate  confidence  interval  is  obtained  using  s. 
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Bayesian  Statistical  Inference 

Bayes'  Theorem.  Bayes'  theorem  (6:174-5)  can  be 
written  as: 

f(x\e)gi0) 

g(e|x)  =  jF(x) 

where  0  is  the  population  parameter 

g(0|x)  is  the  posterior  probability  distribution  of  0 
given  sample  data  x 

g(0)  is  the  prior  probability  distribution  of  0 

f (x|0)  is  the  probability  of  observing  sample  data  x 
given  the  true  parameter  0 

f(x)  is  the  marginal  probability  distribution  of  x. 
f(x)  can  be  obtained  by 

i  f(x\0)g(0)d6,  6  continuous 

f(x)  = 

T.fix\0)g{9),  0  discrete. 

The  Prior  Distribution.  Bayesian  inference  is 
appropriate  only  when  the  prior  assessment  of  0  is  accurate. 


Great  care  must 

be 

exercised  when 

selecting 

the 

prior 

density  g(0). 

The 

prior  density 

can  take 

any 

form. 

However,  natural  conjugate  priors  and  noninformative  priors 
(also  called  priors  of  ignorance)  are  used  in  most  practical 
applications.  When  using  a  natural  conjugate  prior,  the 
posterior  distribution  is  in  the  same  family  as  the  prior 
distribution  and  can  be  obtained  simply  by  using  analytical 
tools.  When  there  is  little  or  no  information  on  the  value 
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of  a  parameter,  a  noninformative  prior  is  used.  The  most 
commonly  used  is  the  uniform  distribution.  For  a  detailed 
discussion  of  the  noninformative  prior,  see  Press  (8:46-51). 
For  the  continuous  prior  distribution,  discretization  is 
often  used  to  allow  for  more  tractable  data  manipulation. 
Chay  (1)  argues  that  since  numerical  integration  can  be  done 
using  computer  calculations,  it  is  not  necessary  to 
approximate  continuous  prior  distributions  using 
discretization.  The  process  of  determining  the  form  of  g(6) 
is  difficult  and  controversial.  When  selecting  a  prior,  it 
must  be  both  reasonable  and  justifiable. 

Bayesian  Reliability  Testing  Approaches.  In  this  section, 
current  approaches  in  applying  Bayesian  methods  in 
reliability  testing  will  be  addressed. 

Chen  and  Papadopoulos  (2)  proposed  a  Bayesian  method 
for  a  generalized  exponential  failure  model  for  components 
in  a  parallel  or  series  system.  The  gamma  or  uniform 
distribution  is  used  as  a  prior.  The  generalized  failure 
models  under  consideration  are  the  exponential,  Weibull, 
Rayleigh,  and  as  a  special  case,  the  extreme  value  failure 
model.  Sharma  and  Bhutan!  (11)  considered  a  parallel  system 
composed  of  independent  and  identical  components  whose 
failure  distribution  was  exponentially  distributed.  The 
exponential  distribution  is  used  as  a  prior.  Villacourt  and 
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Mahaney  (12)  used  the  Bayesian  method  to  design  a 
reliability  demonstration  test  on  a  lithography  expose  tool. 
The  failure  times  are  assumed  to  be  exponentially 
distributed.  The  gamma  distribution  is  used  as  a  prior. 

Martz,  Waller,  and  Fickas  (7)  developed  a  Bayesian 
procedure  for  a  series  system  of  binomial  subsystems  and 
components.  Test  or  prior  data  (perhaps  both  or  neither)  at 
the  system,  subsystem,  and  component  level  are  considered. 
The  beta  distribution  is  used  as  a  prior.  The  method  was 
used  to  estimate  the  overall  reliability  of  an  air-to-air 
heat  seeking  missile  system  with  five  major  subsystems  and 
up  to  nine  components  per  subsystem.  Martz  and  Waller  (5) 
extended  the  method  to  cover  the  case  of  a  complex 
series/parallel  systems  of  binomial  subsystems  and 
components.  These  methods  are  very  complex  and  difficult  to 
implement . 

Willits,  Dietz,  and  Moore  (15)  show  that  for  a 
series-system,  reliability  estimates  using  very  small  sets 
(less  than  10  data  points)  of  binomial  test-data,  there  is 
no  clear  advantage  of  using  Bayes  interval  estimation  unless 
the  prior  mean  system  reliability  is  believed  to  be  within 
20  percent  of  the  true  system  reliability.  A  classical 
method  for  very  small  samples,  the  Linstrom-Madden 
estimator,  should  be  used. 
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When  both  prior  and  current  data  are  available,  which 
is  the  case  in  this  research,  the  methods  mentioned  above  in 
general  are  not  applicable.  However,  studying  these  methods 
provided  the  general  knowledge  of  understanding  how  Bayesian 
statistical  methods  worked. 

In  this  research,  the  applicability  of  Bayesian 
statistical  methods  to  the  problem  of  determining  the 
reliability  of  a  series  system,  a  parallel  system,  and  a 
bridge  system  is  being  studied.  The  time  to  failure  of  each 
component  in  each  system  is  assumed  to  be  exponentially 
distributed.  A  basic  exponential  case  and  a  basic  binomial 
case  are  also  evaluated.  There  is  very  little  information 
concerning  the  use  of  Bayesian  methods  to  the  systems 
mentioned  above  in  the  literature.  The  methodology  being 
developed  in  Chapter  III  is  based  on  the  works  of  the 
following  individuals. 

Lemaster  (4)  studied  the  applicability  of  Bayesian 
method  to  the  problem  of  determining  cruise  missile 
component  reliability.  The  exponential  and  binomial 
distributions  were  used  for  the  failure  rate.  For  the 
exponential  case,  a  gamma  prior  is  used.  For  the  binomial 
case,  a  beta  prior  is  used. 

Sharma  and  Bhutani  (10)  considered  the  case  for  a 
series  system  when  the  failure  time  for  the  system  is 
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exponentially  distributed.  The  gaitma  distribution  is  used 
as  a  prior. 

Frickenstein  (3)  studied  the  case  for  a  simple  system 
(one  component)  when  the  failure  time  for  the  system  is 
assumed  to  be  exponentially  distributed.  The  gamma 
distribution  is  used  as  a  prior. 

The  text  book,  Bayesian  Reliability  Analysis,  written 
by  Martz  and  Waller  (6)  is  the  standard  reference  work  on 
Bayesian  statistical  methods  and  will  be  used  extensively  in 
this  research  effort.  The  binomial  model  developed  in 
Chapter  III  is  based  on  the  methods  developed  by  Lemaster 
(4)  with  modifications  based  on  Martz  and  Waller  (6)  .  The 
exponential  model  developed  in  Chapter  III  is  based  on  the 
methods  developed  by  Lemaster  (4)  and  by  Frickenstein  (3) 
with  modifications  based  on  Martz  and  Waller  (6) . 

Reliability  analysis  methods  for  a  series  system,  a 
parallel  system,  and  a  bridge  system  (with  the  time  to 
failure  of  each  component  in  each  system  is  assumed  to  be 
exponentially  distributed)  are  developed  in  Chapter  III 
utilizing  the  exponential  model.  For  the  series  system,  the 
Bayesian  methods  developed  by  Sharma  and  Bhutan!  (10)  will 
also  be  modified  for  use  in  this  research.  The  next  chapter 
will  discuss  the  methodology  used  in  this  research  to 
accomplish  the  research  objectives  stated  in  Chapter  I. 
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III.  Methodology 


Overview . 

This  chapter  details  the  methodology  used  to  accomplish 
the  research  objectives  stated  in  Chapter  I.  The 
methodology  developed  in  this  chapter  is  based  on  the 
assumption  that  both  DT&E  (prior)  and  OT&E  (current)  data 
are  available.  A  basic  binomial  model  and  a  basic 
exponential  model  are  developed  first.  The  basic 

exponential  model  is  then  applied  to  the  problem  of 
determining  the  reliability  of  a  series  system,  a  parallel 
system,  and  a  bridge  system.  The  time  to  failure  of  each 
component  in  each  system  is  assumed  to  be  exponentially 
distributed.  Classical  and  Bayesian  point  estimates  and 
confidence  intervals  obtained  from  each  system  will  be 
compared  in  Chapter  IV. 

Data  Source. 

For  the  basic  binomial  model,  a  hypothetical  set  of 
missile  test  data  will  be  used  to  validate  the  methodology. 
For  the  basic  exponential  model,  the  model  is  applied  to 
data  obtained  by  AFOTEC  during  a  reliability  testing  of  a 
selected  simple  system  (system  is  treated  as  one  component) 
with  the  time  to  failure  of  the  system  is  assumed  to  be 
exponentially  distributed.  This  is  the  only  set  of  true 
data  used  in  this  research.  Simulated  data  will  be  used  to 
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validate  the  methodology  for  a  series  system,  a  parallel 
system,  and  a  bridge  system.  This  research  uses  a  computer 
program  called  RAPTOR  (Rapid  Availability  Prototyping  for 
Testing  Operational  Readiness),  which  is  developed  by 
AFOTEC,  to  simulate  the  true  system  (series,  parallel, 
bridge)  mentioned  above.  RAPTOR  is  verified  and  validated 
by  AFOTEC  and  is  an  excellent  tool  for  use  in  this  research 
effort . 

The  Binomial  Model. 

This  model  is  appropriate  for  systems  or  components 
that  are  assumed  to  have  binomial  failure  distribution.  The 
classical  binomial  method  is  given  first.  Then  the  Bayesian 
binomial  method  is  given  based  on  the  assumption  that  both 
prior  data  (DT&E  data)  and  current  data  (OT&E  data)  are 
available.  The  model  is  then  applied  to  a  hypothetical 
missile  launch  problem  in  Chapter  IV. 

The  Classical  Binomial  Method. 

The  binomial  distribution  is  one  of  the  commonly  used 
distribution  in  reliability.  Consider  a  missile  system  that 
is  being  tested  with  the  outcome  of  each  test  as  the  missile 
success  or  failure  to  launch.  Assuming  a  missile  system  has 
the  binomial  distribution,  the  point  estimate  for  the 
probability  of  success  p  (which  is  the  reliability  of  a 
missile  system)  and  80%  confidence  interval  about  p  will  be 
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computed.  For  the  binomial  distribution,  the  point  estimate 
using  the  maximum  likelihood  estimator  is 

P.(^)  (5) 

where  X  is  the  number  of  successes  from  a  set  of  n  missile 
units  placed  into  a  test  of  given  duration  (6:53).  The 
confidence  limits  can  be  approximated  by  using  the  F 
distribution.  The  100(1  -  y)%  two-sided  confidence  interval 
(TCI)  is  (6:56) 

f _ X _  (X+l)Fi-y/2(2X+2 , 2  n-2X)  \ 

V  X4<n-X+l)Fi_j,/2(2n-2X+2,2X) '  (n-X)4<X+l)Fi-y/2(2X+2,2n-2X)y  • 

The  Bayesian  Binomial  Method. 

Since  they  are  relatively  easy  to  compute  and  can  be 
adapted  to  a  number  of  situations,  the  family  of  beta 
distributions,  B(x,n),  are  used  for  the  binomial  case  to 
calculate  both  the  prior  and  posterior  distributions.  For 
example,  the  analyst's  knowledge  and  experience  with  the 
system  is  utilized  to  select  x  and  n.  Thus  the  analyst  can 
increase  or  decrease  the  importance  of  the  prior  on  the 
resulting  posterior  distribution  (4:39) . 

Bayesian  Binomial  Prior.  For  the  binomial 
sampling,  a  B(xo,no)  distribution  can  be  used  to  maximize 
flexibility  and  make  it  easier  to  calculate  the  posterior 
distribution  since  the  B(xo,no)  distribution  is  the  natural 
conjugate  prior  distribution.  The  parameters  Xo  and  no  can 
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be  interpreted  as  the  pseudo  number  of  survivals  and  pseudo 
sample  size  for  the  prior  life  test.  Here  "pseudo"  should 
be  thought  of  as  meaning  "pretended"  (6:265) .  When  the  true 
prior  distribution  is  not  of  the  beta  family,  Weiler  (13) 
showed  that  assuming  a  B(xo,no)  distribution  is  acceptable 
in  many  practical  applications.  Weiler  was  able  to  show 
that  severe  deviations  from  the  beta  prior  parameter  values 
will  result  in  minor  changes  in  the  corresponding  posterior 
distributions.  DT&E  data  will  be  used  to  approximate  xo  and 
no  that  fit  the  B{xo,no)  distribution.  For  the  missile 
launch  problem  in  Chapter  IV  where  the  outcome  of  each  test 
is  a  missile  success  or  failure  to  launch,  Xo  is  the  number 
of  successes  during  DT&E  test  and  similarly  no  is  the  total 
number  tested. 

Bayesian  Binomial  Posterior.  For  the  B(xo,no) 
prior  distribution,  the  resulting  posterior  distribution  is 
also  a  beta  distribution  of  the  form  B(x  +  xo,n  +  no),  where 
X  is  the  observed  successes  and  n  is  the  sample  size  from 
the  test  being  used  to  update  the  prior.  The  probability 
density  function  (6:266)  is  given  as 

g(p|x;x.,n.).  [  r(x-^xo)?(nt”°o-x-x,))  -p)«— fV 
where  0  <  p  <  1. 
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Bayesian  Binomial  Point  Estimate.  The  mean  of  the 
posterior  distribution  is  the  Bayesian  point  estimator. 
Under  the  squared  error  loss  function,  the  Bayesian  point 
estimator  is  (6:267) 

p  =  E(p|x;xo,no)  =  ( •  (8) 

The  loss  function  measures  the  error  sustained  by 
estimating  p  by  p.  The  squared  error  loss  function  often 
used  in  practical  application  is  (p  -  p)^  (14:56). 

Bayesian  Binomial  Probability  Interval.  The 
Bayesian  methods  use  "probability  intervals"  instead  of 
"confidence  intervals".  Recall  that  the  classical  method 
assumes  that  the  parameter  p  estimated  is  an  unknown 
constant.  Thus  the  estimators  of  the  end  points  of  the 
confidence  interval  ‘associated  with  p  estimated  are  random 
variables.  A  classical  confidence  interval  is  not  an 
explicit  probability  statement  about  p.  However,  it  is  "the 
probability  that  the  interval  estimator  will  generate  an 
interval  that  will  contain  the  true  value  of  p"  (14:57).  In 
contrast,  Bayesian  method  considers  p  to  be  a  random 
variable.  Thus  a  Bayesian  probability  interval  is  an 

explicit  probability  statement  about  the  value  of  p  (6:208) . 

Using  the  F(ni,n2)  distribution,  the  Total  Bayesian 
Probability  Interval  (TBPI)  can  be  calculated.  A  80% 
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confidence  interval  will  be  calculated.  For  the  100(1  - 


TBPI,  the  upper  and  lower  interval  endpoints  are  (6:270) 


_ X-fXQ _ _ 

x+xo+(n+no-x-xo)Fi_y/2(2n+2no-2x-2xo/  2x+2xo) ' 
(x+xo)Fi-3>/2(2x+2xo  ,  2n+2no-2x-2xo) 
n+no-x-xo+(x+xo)Fi-y/2(2x+2xo,2n+2no-2x-2xo)  j 


(9) 


The  Exponential  Model. 

This  model  is  appropriate  for  systems  or  components 
that  are  assumed  to  have  exponential  failure  distribution. 
The  classical  exponential  method  is  given  first.  Then  the 
Bayesian  exponential  method  is  given  based  on  the  assumption 
that  both  prior  data  (DT&E  data)  and  posterior  data  (OT&E 
data)  are  available.  In  Chapter  IV,  the  model  is  applied  to 
data  obtained  by  AFOTEC  during  a  reliability  testing  of  a 
selected  simple  system  with  the  time  to  failure  of  the 
system  is  assumed  to  be  exponentially  distributed. 

The  Classical  Exponential  Method. 

The  exponential  distribution  is  the  most  widely  used 
distribution  in  reliability.  Assuming  the  time  to  failure 
of  a  system  or  component  has  the  exponential  distribution, 
point  estimates  for  the  failure  rate  X  and  80%  confidence 
intervals  on  X  will  be  computed.  For  the  exponential 
distribution,  the  point  estimate  of  X  using  the  maximum 
likelihood  estimator  is 
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where  f  is  the  number  of  failures  and  T  is  the  total  test 


time.  The  number  of  failures,  f,  is  presumed  to  be  a  random 
variable.  The  confidence  limits  can  be  approximated  by 
using  the  chi  square  distribution.  The  100(1  -  y)  % 

two-sided  confidence  interval  (TCI)  is  (6:122) 
f  Xy/2(?f) 

^  22’  '  2T  J  .  (11) 


For  the  reliability  function  R(t)  =  exp(-A,t),  the 
reliability  estimator  of  R(t)  for  a  mission  time  t  is 

R(t)  =  exp(-lt)  (12) 

where  t  is  the  mission  length. 

Equation  12  is  also  used  to  calculate  an  interval 
estimate  of  R{t)  by  using  the  lower  and  upper  bound  of  A 

obtained  from  Equation  11  in  place  of  A. 

The  Bayesian  Exponential  Method. 

When  examining  a  continuously  operating  electronic 
system,  the  failure  time  is  often  assumed  to  be 
exponentially  distributed.  The  distribution  of  the  number 
of  failures,  f,  in  fixed  total  test  time,  T,  can  be 


where  A,,T  >0,  f  =  0,  1,  2,  ...  and  assuming  that  failed 

items  are  replaced.  This  is  known  as  Poisson  sampling. 

Bayesian  Exponential  Prior.  For  the  Poisson 
sampling,  a  gamma  prior  distribution,  G(a,P),  may  be  used  to 
maximize  flexibility  and  make  it  easier  to  calculate  the 
posterior  distribution  since  the  G(a,P)  distribution  is  the 
natural  conjugate  prior  distribution  (6:289).  The  G(a,p) 
prior  distribution  is  probably  the  most  widely  used  and  has 
the  following  probability  density  function: 

g(>.;a,p)  =  [  J  (14) 


where  X,  a,  and  P  >  0  and  r(a)  is  the  gamma  function  of  a. 

The  shape  parameter  a  can  be  interpreted  as  the  pseudo 
number  of  failures  in  a  prior  life  test  of  duration  P  pseudo 
time  units  (6:289).  The  mean  and  variance  of  the  gamma 
distribution  are 


E(;.)  = 

and  V(X)  =(^)  . 


(15) 

(16) 


E(X)  and  V(X)  can  be  estimated  using  the  sample  mean 


and  variance  of  X  respectively  with  E(X)  =  A  and  V(X)  =  s| . 
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Equations  2  and  3  (page  9)  are  used  to  compute  A  and  s| . 
Thus  a  and  P  can  be  computed  as 


(17) 


(18) 


Bayesian  Exponential  Posterior.  For  the  G(a,P) 
prior  distribution,  the  resulting  posterior  distribution  is 
a  G{f  +  a,T  +  p) ,  where  f  is  the  number  of  failures  and  T  is 
the  total  test  time.  The  probability  density  function 
(6:290)  is  given  as 


g(A,|f;a,p) 


I  r(jF+a)  J 


(19) 


where  X  >  0.  The  parameter  {f  +  a)  is  referred  as  the 
combined  number  of  failures  and  (T  +  P)  is  the  combined 
total  test  time. 

Bayesian  Exponential  Point  Estimate.  The  mean  of 
the  posterior  distribution  is  the  Bayesian  point  estimator. 
Under  the  squared  error  loss  function,  the  Bayesian  point 
estimator  is 


A  =  E(Xlf;a,p)  =  (20) 

where  (/  +  a)  is  the  combined  number  of  failures  and  (T  +  P) 
is  the  combined  total  test  time  (6:292) . 
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Bayesian  Exponential  Probability  Interval. 


For 


the  100(1  -  Y) %  TBPI,  the  upper  and  lower  interval  endpoints 
are  (6:294) 

f  4/2(2 -f+Sa)  ;ff_j,/2(2r+2a)^ 

2(t+)3)  '  2(t+fi)  )  (21) 


For  the  reliability  function  R(t)  =  expi-Xt),  the 

Bayesian  estimator  of  R(t)  for  a  mission  time  t  is  (3:8) 

00 

R(t)  =  E[R(t)|f]  =  lexp(-Xt)g(2.\f;a,p)dX 

0 


where  f  is  the  number  of  failures,  T  is  the  total  test  time 
and  t  is  the  mission  time. 

Another  way  of  computing  the  Bayesian  estimator  of  R(t) 

is  R{t)  =  exp(-At).  Both  methods  should  yield  very  similar 
result.  Mathematically  speaking.  Equation  22  is  selected 
for  use  since  it  is  directly  derived.  The  reliability 
function  R(t)  =  exp(-X,t)  is  also  used  to  calculate  an 
interval  estimate  of  R{t)  by  using  the  lower  and  upper 
bound  of  7.  obtained  from  Equation  21  in  place  of  7. 

Bayesian  Neibull  Method.  If  the  true  distribution 
of  the  time  to  failure  is  not  exponentially  distributed,  the 
Bayesian  method  can  still  be  developed  but  it  is  a  lot 
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harder  to  do  the  number  calculation.  Another  frequently 
used  time  to  failure  distribution  is  the  Weibull 
distribution.  The  Bayesian  Weibull  method  described  below 
shows  how  difficult  it  can  be  to  do  this  calculation. 

If  the  prior  distribution  of  X  is  the  two-parameter 
Weibull  distribution  with  probability  density  function 
(6:300) 

g(>,;a,P)  =  (  a )  (  o’)  ^  ®xp( "(  a)  >  0  (23) 

where  a  is  the  scale  parameter  and  P  is  the  shape  parameter. 

For  the  Weibull  prior  distribution,  the  resulting 
posterior  probability  density  function  (6:300)  is  given  as 


g(X,|f;a,P)  =  "ss 

\  0 


X  >  0 


where  f  is  the  number  of  failures  and  T  is  total  test  time. 

For  the  Weibull  prior  distribution,  the  resulting 
posterior  mean  (6:300)  is  given  as 

1  =  E(X|f;a,p)  =  -  X  >  0.  (25) 

f  Q-XT-iX/a)f  (jX, 

Vo  J 


Both  equations  given  above  must  be  numerically  evaluated. 
Due  to  the  difficulty  associated  with  numerical  evaluation, 
the  Weibull  prior  will  not  be  used  in  this  thesis. 
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Reliability  Analysis  of  a  Series  System. 

Consider  a  system  composed  of  k  independent  components 
in  series  configuration.  The  time  to  failure  of  each 
component  is  assumed  to  be  exponentially  distributed.  Let 
Xi,  "kz,  .  . . ,  Xk  be  the  respective  constant  failure  rates  of 
each  component.  The  failure  rate  for  the  system  is  (10:761) 

k 

k  =  LAi  .  (26) 

i=l 

This  shows  that  for  a  system  composed  of  k  independent 
exponentially  distributed  components  in  series 
configuration,  the  distribution  of  the  failure  time  of  the 
system  is  also  exponential  with  failure  rate  k. 

The  Bayesian  Method  1  (10:761-763) . 

Bayesian  Prior  Analysis .  Assuming  that  the  prior 
belief  about  the  failure  rates  Xi  (i  =  1,  2,  ...,  k)  is 


exponentially  distributed  with  parameter  6.  Then  the  prior 

k 

belief  about  k  =  Jlki  is  defined  by  a  gamma  distribution 


with  probability  density  function  (10:761) 

g(>.)  =  >  0.  (27) 

The  mean  and  the  variance  for  the  failure  rate  of  the  system 
are  E(X)  =  k/0  and  V{k)  =  k/G^. 
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Bayesian  Posterior  Analysis .  For  a  system  with 
failure  time  being  exponentially  distributed  with  failure 
rate  X,  then  for  the  accumulated  test  time  T,  the  number  of 
system  failures  f  is  Poisson  distributed  with  pdf  (10:761) 

P{f\X)  =  - J  ^  =  0,  1,  2,  ...  (28) 

Thus  the  posterior  probability  density  function  for  X 
given  f  failures  over  a  time  interval  (0,T)  with  respect  to 
a  gamma  prior  given  above  is  (10:761-762) 

P(>.|f)  =  -  nf+k) - J  (T+0),  (f+k)  >  0  (29) 


which  is  also  a  gamma  with  parameters  (T+0)  and  (f+k) .  The 


Bayesian  estimator  of  X  given  f  failure 
testing  is  (10:762) 

3=E(>.lf)  =  (^)  . 


in  T  hours  of 


(30) 


For  the  reliability  function  R(t)  =  exp(-Xt),  the 
Bayesian  estimator  of  R(t)  for  a  mission  time  t  given  f 
failures  in  T  hours  of  testing  is  (10:762) 

r  \ 


R{t)  =  E[R(t)  \f] 


1  + 


(3i: 


■(  T+eJ  J 

Using  The  Exponential  Model.  This  method  uses  the 
Exponential  Model  given  above  to  calculate  the  classical  and 
Bayesian  reliability  point  estimates,  Ri  (i  =  1  ..  k) ,  and 
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the  classical  and  Bayesian  component  reliability  confidence 
intervals,  (LBi,UBi)  .  LBi  is  the  interval  lower  bound  and  UBi 
is  the  interval  upper  bound  for  each  component.  The  next 
step  is  to  aggregate  the  appropriate  numbers  together  to 
obtain  the  classical  and  Bayesian  system  reliability  point 
estimates  and  the  classical  and  Bayesian  system  reliability 
confidence  interval,  (LBs,UBs) .  The  system  reliability 
point  estimate  is 

Jr 

R3=niRi.  (32) 

The  system  reliability  confidence  interval  is  calculated  as 

LBs  =  (33) 

Jr 

and  UBs  =  •  (34) 

The  Exponential  Model  is  also  used  to  calculate  the 
classical  and  Bayesian  reliability  point  estimates  and  the 
classical  and  Bayesian  confidence  intervals  for  the  system 
using  component  level  data  aggregated  to  system  level  data. 
This  is  possible  since  for  a  series  system  with  exponential 
components,  the  failure  rate  for  the  system  is  also 
exponential  and  is  the  sum  of  the  components'  failure  rates. 
Reliability  Analysis  of  a  Parallel  System. 

Consider  a  system  composed  of  k  independent  components 
in  parallel  configuration.  The  time  to  failure  of  each 
component  is  assumed  to  be  exponentially  distributed.  Let 
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Xi  (i  =  1  ..  k)  be  the  respective  constant  failure  rate  of 
each  component.  This  method  uses  the  Exponential  Model 
given  above  to  calculate  the  classical  and  Bayesian 
reliability  point  estimates,  Ri  (i  =  1  ..  k) ,  and  the 

classical  and  Bayesian  component  reliability  confidence 
intervals,  (LBi,UBi)  .  LBi  is  the  interval  lower  bound  and  UBi 
is  the  interval  upper  bound  for  each  component.  The  next 
step  is  to  aggregate  the  appropriate  numbers  together  to 
obtain  the  classical  and  Bayesian  system  reliability  point 
estimates  and  the  classical  and  Bayesian  system  reliability 
confidence  intervals,  (LBs,UBs)  -  The  system  reliability 
point  estimate  is 

Rs  =  1  -  ni(l--Ri).  (35) 

The  system  reliability  confidence  interval  is  calculated  as 

LBs  =  1  -  n(l-iBi)  (36) 

k 

and  UBs  =  1  -  n(l-(7Bi).  (37) 

For  the  case  of  independent  and  identical  components. 
The  Exponential  Model  is  also  used  to  calculate  the 
classical  and  Bayesian  reliability  point  estimates,  Rc,  and 
the  classical  and  Bayesian  confidence  intervals,  (LBc,UBc) 
for  the  component  using  aggregated  component  level  data. 
This  is  possible  since  the  failure  rate  for  each  component 
is  the  same.  The  system  reliability  point  estimate  is 
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R,  =  1  -  (1  -  Rc)’'. 


(38) 


The  system  reliability  confidence  interval  is  calculated  as 
LBs  =  1  -  (1  -  LBc)’'  (39) 

and  UBs  =  1  -  (1  -  UBc)"'.  (40) 

Reliability  Analysis  of  a  Bridge  System. 

Consider  a  system  composed  of  5  independent  components 
given  in  Figure  1  below.  The  time  to  failure  of  each 
component  is  assumed  to  be  exponentially  distributed. 


Figure  1.  Five  Components  Bridge  System. 

Let  (i  =  1  ..  5)  be  the  respective  constant  failure 
rate  of  each  component.  This  method  used  the  Exponential 
Model  given  above  to  calculate  the  classical  and  Bayesian 
reliability  point  estimates,  Ri  (i  =  1  ..  5),  and  the 
classical  and  Bayesian  component  reliability  confidence 
intervals,  (LBi,UBi)  .  LBi  is  the  interval  lower  bound  and  UBi 
is  the  interval  upper  bound  for  each  component.  The  next 
step  is  to  aggregate  the  appropriate  numbers  together  to 
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obtain  the  classical  and  Bayesian  system  reliability  point 
estimates  and  the  classical  and  Bayesian  system  reliability 
confidence  interval,  (LBs,UBs)  .  The  system  reliability 
point  estimate  is  (9:422) 

Rs  =  RiR4  +  R2R5  +  R1R3R5  +  R2R3R4  ~  R1R2R3R4  “  R1R2R3R5  “ 
RiR2R4R5  “  R1R3R4R5  “  R2R3R4R5  +  2RiR2R3R4R5«  (41) 

The  system  reliability  confidence  interval  is  calculated  as 

LBs  ~  LB1LB4  +  LB2IJB5  +  LB1LB3LB5  +  LB2LB3LB4  “  LB1LB2LB3IJB4 
“  LB1LB2LB3LB5  ~  LB1LB2LB4LB5  ~  LB1LB3LB4LB5  ~  LB2IJB3LB4LB5  + 
2LB1LB2LB3LB4LB5  (42) 

and  UBs  =  UB1UB4  +  UB2UB5  +  UB1UB3UB5  +  UB2UB3UB4  - 
UB1UB2UB3UB4  -  UB1UB2UB3UB5  -  UB1UB2UB4UB5  -  UB1UB3UB4UB5  - 
UB2UB3UB4UB5  +  2UB1UB2UB3UB4UB5.  (43) 

For  the  case  of  independent  and  identical  components. 
The  Exponential  Model  is  also  used  to  calculate  the 
classical  and  Bayesian  reliability  point  estimates,  Ro,  and 
the  classical  and  Bayesian  probability  intervals,  (LBc,UBc) 

for  the  component  using  aggregated  component  level  data. 
This  is  possible  since  the  failure  rate  for  each  component 
is  the  same.  The  system  reliability  point  estimate  is 
calculated  using  Equation  41  above  with  Rc  replacing  Ri  (i  = 
1  ..  5).  The  system  reliability  confidence  interval  is 

calculated  using  Equations  42  and  43  with  LBc  replacing  LBi 
(i  =  1  ..  5)  and  UBc  replacing  UBi  (i  =  1  ..  5). 
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IV.  Results 


Overview. 

This  chapter  describes  and  discusses  the  results  of  the 
calculations  using  the  methodology  in  Chapter  III.  Results 
from  a  basic  binomial  case  and  a  basic  exponential  case  are 
given  first.  Then  results  from  a  series  system,  a  parallel 
system  with  independent  and  identical  components,  and  a 
parallel  system  with  different  independent  components  are 
presented.  Finally,  results  from  a  bridge  system  with 
independent  and  identical  components  and  a  bridge  system 
with  different  independent  components  are  given.  The  time 
to  failure  of  each  component  in  each  system  is  assumed  to  be 
exponentially  distributed.  For  detailed  calculations,  see 
the  appropriate  appendix.  Calculations  will  be  shown  in 
detail  the  first  time  presented  in  the  appropriate  appendix. 
Classical  and  Bayesian  point  estimates  and  confidence 
intervals  obtained  from  each  system  will  be  presented  and 
discussed  in  this  chapter. 

The  Binomial  Model. 

Suppose  the  following  data  are  obtained  from  a  missile 
test  program.  For  the  DT&E  portion  of  the  test,  two 
missiles  failed  out  of  20  missiles  launched.  For  the  OT&E 
portion  of  the  test,  one  missile  failed  out  of  12  missiles 
launched.  Assuming  the  missile  system  failures  may  be 
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modeled  with  the  binomial  distribution,  point  estimates  for 
the  probability  of  success  p  (which  is  the  reliability  of 
the  missile  system)  and  80%  confidence  interval  about  p  will 
be  computed.  From  the  data  above,  the  following  results  are 
obtained  (for  detailed  calculation,  see  Appendix  A) : 


Table  2:  Summary  of  Binomial  Results 


Classical 

Bayesian 

DT&E 

P 

80%  Interval 
Interval  Width 

90.00% 

(75.52%<p<97.31%) 

21.79% 

OT&E 

P 

80%  Interval 
Interval  Width 

91.67% 

(71.25%^p^99.13%) 

27.88% 

90.63% 

(83.73%^p£96.39%) 

12.66% 

Table  2  shows  that  the  point  estimates  for  the 
probability  of  success,  p,  are  all  very  similar.  The  bounds 
using  classical  method  are  quite  large.  You  are  80% 
confident  that  the  true  p  is  between  75.52%  and  97.31%  for 
DT&E  data  and  between  71.25%  and  99.13%  for  OT&E  data.  The 
strength  of  Bayesian  analysis  lies  in  the  posterior  analysis 
as  the  table  shows.  Although  the  estimates  of  p  is  very 
close  to  the  estimates  of  p  using  classical  method,  the 
Bayesian  method  yields  a  much  tighter  range  for  the  point 
estimate.  The  probability  is  80%  that  the  true  p  is  between 
83.73%  and  96.39%.  Even  though  the  width  of  the  probability 
interval  is  still  large,  it  is  a  significant  improvement 
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compared  to  the  width  of  the  intervals  obtained  using 
classical  method. 

The  Exponential  Model. 

The  following  data  are  obtained  from  a  reliability  test 
program  of  a  selected  simple  system  by  AFOTEC  (assuming  the 
time  to  failure  of  the  terminals  being  tested  is 
exponentially  distributed) .  For  the  OT&E  portion  of  the 
test,  five  terminals  were  tested  for  a  period  of  seven 
months.  The  number  of  operational  hours  (Ti)  and  the  number 
of  critical  failures  per  terminal  (Xi)  were  recorded  and  are 
given  in  Table  3. 


Table  3:  OT&E  Data  (2 

):6) 

Terminal  i 

Ti  Hours 

Xi  Failures 

1 

184.3 

12 

2 

232.6 

9 

3 

172.8 

9 

4 

284.4 

17 

5 

264.8 

18 

For  the  operational  assessment  (OA)  portion  of  the 
test,  two  terminals  were  tested.  Again  the  number  of 
operational  hours  (Ti)  and  the  number  of  critical  failures 
per  terminal  (Xi)  were  recorded  and  are  given  in  Table  4. 


Ta] 

ble  4 :  OA  Data  ( 3 : 

7) 

Terminal  i 

Ti  Hours 

Xi  Failures 

1 

68.5 

2 

2 

60.9 

2 
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The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  Since  the  time  to  failure 
of  the  terminals  being  tested  is  assumed  to  be  exponentially 
distributed,  a  gamma  prior  distribution  is  assumed.  From 
the  data  above,  the  following  results  are  obtained  (for 
detailed  calculation,  see  Appendix  B) : 
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Table  5:  Summary  of  Exponential  Results 


36 
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Figure  2  shows  that  the  OA  reliability  estimates  using 
the  Bayesian  method  are  lower  (more  conservative)  than  the 
OA  classical  estimates  (which  consider  only  the  4  failures 
in  129.4  hours)  for  the  six  different  values  of  t. 

Figure  3  shows  that  the  strength  of  Bayesian  analysis 
lies  in  the  posterior  analysis.  In  all  six  cases,  the 
Bayesian  method  yields  a  tighter  range  for  the  point 
estimate.  The  difference  between  the  OA  Bayesian 
probability  interval  width  and  OA  classical  confidence 
interval  width  ranges  from  1.34%  to  4.66%  for  the  mission 
times  being  considered.  As  the  mission  time  increases,  the 
widths  of  the  confidence  interval  and  the  probability 
interval  both  increase  and  the  gaps  between  the  two 
intervals  get  larger. 

Reliability  Analysis  of  a  Series  System. 

Consider  a  system  composed  of  k  =  3  independent 
components  in  series  configuration.  The  time  to  failure  of 
each  component  is  assumed  to  be  exponentially  distributed. 
Let  X.1  =  0.05  failure  per  hour,  Xz  =  0.067  failure  per  hour, 
and  ^3  =  0.05  failure  per  hour  be  the  respective  constant 
failure  rates  of  each  component.  The  failure  rate  for  the 
system  is 

3 

X  =  ^■^i=  0.05  +  0.067  +  0.05  =  0.167  failure  per  hour. 
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Thus  the  mean  failure  time  for  the  system  is 


'  (i)'  (oTist)  - « 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  and  3  are 
Ri  =  Rs  =  exp  (-0.05)  =  95.12% 

R2  =  exp (-0.067)  =  93.52%. 

The  reliability  of  the  system  can  be  calculated  as 
Rs  =  R1XR2XR3  =  95.12%x93.52%x95.12%  =  84.62% 
or  Rs  =  exp(-7,t)  =  exp (-0.167)  =  84.62%. 

Using  a  computer  program  called  RAPTOR  (Rapid 
Availability  Prototyping  for  Testing  Operational  Readiness), 
which  is  developed  by  AFOTEC,  to  simulate  the  true  system, 
the  following  data  are  obtained  (to  review  the  simulated 
data,  see  Appendix  C) ; 


Table 

6:  System  Level  Dl 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

499.86 

80 

2 

449.53 

78 

3 

384.76 

46 

4 

543.47 

97 

5 

591.88 

75 

Table 

7 :  System  Level  01 

?  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

87.86 

13 

2 

122.92 

27 

39 


Table  8 :  Component  Level  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 
Component  1 

Xi  Failures 
Component  2 

Xi  Failures 
Component  3 

1 

499.86 

24 

28 

28 

2 

449.53 

20 

37 

21 

3 

384.76 

13 

20 

13 

4 

543.47 

29 

39 

29 

5 

591.88 

21 

33 

21 

Table  9:  Component  Level  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 
Component  1 

Xi  Failures 
Component  3 

1 

87.86 

4 

5 

4 

2 

122.92 

9 

10 

8 

The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  From  the  data  above,  the 
following  results  are  obtained  (for  detailed  calculation, 
see  Appendix  D) : 
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Table  10:  System  Reliability  Using  System  Level  Data 


*  The  Value  in  ( )  in  Bayesianl  collumn  is  combined  DT  & 
OT. 
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Table  11:  System  Reliability  Using  Component  Level  Data 
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Figure  4:  OT&E  Reliability  using  Bayesian  Method  1 


OT&E  Reliability  Using  Bayesian  Method  1 


Figure  4  shows  that  using  Bayesian  method  1  (Bayesian 
method  1  only  gives  a  reliability  estimate) ,  the  reliability 
estimates  using  only  OT&E  data  are  lower  than  the  true 
reliability.  The  gaps  between  the  true  reliability  and  the 
OT&E  reliability  estimates  for  the  mission  times  being 
considered  range  from  3.01%  to  6.01%.  Since  the  true 
reliability  of  the  system  is  known,  you  would  want  the 
reliability  estimates  of  the  system  to  be  as  close  to  the 
true  reliability  as  possible.  The  results  obtained  above 
using  only  OT&E  data  are  not  very  good. 

The  reliability  estimates  using  combined  OT&E  and  DT&E 
data  are  higher  than  the  true  reliability.  The  gaps  between 
the  true  reliability  and  the  reliability  estimates  for  the 
mission  times  being  considered  range  from  0.91%  to  2%.  The 
deviations  from  the  true  reliability  are  still  significant. 
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Howeve'r,  the  results  are  much  better  than  just  using  OT&E 
data.  The  rest  of  the  the  results  below  used  the 
methodology  for  the  exponential  model  given  in  Chapter  III. 


Figure  5:  Classical  &  Bayesian  Reliability  Estimates 


Classical  &  Bayesian  Reliability  Estimates 


X  OT  System  Ciassicai 
X  OT  System  Bayesian 
O  OT  Conrponent  Classica! 
o  OT  Component  Bayesian 
A  TrueR 

O  DT  System  Classical 
*  DT  Conrponent  Qasstcal 


Figure  5  shows  that  the  DT&E  classical  reliability 
estimates  are  higher  than  the  true  reliability.  The 
deviations  range  from  1.26%  to  2.74%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be 
significant.  The  OT&E  reliability  estimates  using  the 
Bayesian  method  are  higher  than  the  OT&E  classical  estimates 
and  closer  to  the  true  reliability  for  the  five  different 
values  of  t  at  both  the  system  level  and  component  level. 
The  OT&E  reliability  estimates  using  the  classical  method 
are  lower  than  the  true  reliability.  The  deviations  range 
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from  1.9%  to  4%  for  the  mission  times  being  considered. 
Again,  the  deviations  are  considered  to  be  significant. 

The  OT&E  reliability  estimates  using  the  Bayesian 
method  are  all  very  close  (less  than  1%)  to  the  true 
reliability.  The  OT&E  reliability  estimates  using  system 
level  data  underestimated  the  true  reliability  with 
deviations  ranging  from  0.31%  to  0.58%  for  the  mission  times 
being  considered.  Using  component  level  data,  the  OT&E 
reliability  estimates  overestimated  the  true  reliability 
with  deviations  ranging  from  0.36%  to  0.83%  for  the  mission 
times  being  considered. 

Figure  6:  Comparison  of  OT&E  Interval  Widths 

(System  Level) 


Comparison  of  OT&E  Interval  Widths 


M  System  OT  Classical 
X  System  OT  Bayesian 
o  System  Different 


0.5  1  1.6  2  2.5  3  3.5 

Mssk)n  Time 
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Figure  7:  Comparison  of  OT&E  Interval  Widths 
(Component  Level) 


Comparison  of  OT&E  Interval  Widths 


Figures  6  and  7  show  the  strength  of  Bayesian  analysis 
lies  in  the  posterior  analysis.  In  all  five  cases,  the 
Bayesian  method  yields  much  tighter  ranges  for  the  point 
estimate  at  both  the  system  and  component  level.  At  the 
system  level,  the  difference  between  the  OT&E  Bayesian 
probability  interval  width  and  OT&E  classical  confidence 
interval  width  ranges  from  1.97%  to  3.72%  for  the  mission 
times  being  considered.  At  the  component  level,  the 
difference  ranges  from  5.07%  to  9.71%.  System  level 
reliability  estimates  yield  tighter  intervals  than  component 
level  aggregated  to  system.  For  the  mission  times  being 
considered,  the  difference  ranges  from  4.53%  to  9.28%  for 
OT&E  classical  confidence  interval  widths  and  from  1.43%  to 
3.29%  for  OT&E  Bayesian  probability  interval  widths. 
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Reliability  Analysis  of  a  Parallel  System  with  Independent 
and  Identical  Components. 

Consider  a  system  composed  of  3  independent  and 
identical  components  arranged  in  a  parallel  configuration. 
The  time  to  failure  of  each  component  is  assumed  to  be 
exponentially  distributed.  Let  k  =  0.05  failure  per  hour  be 
the  constant  failure  rate  of  each  component.  For  the  true 
system,  using  a  mission  time,  t  of  one  hour,  the  reliability 
of  components  1,  2,  and  3  are 

Rc  =  Ri  =  R2  =  R3  =  exp  (-0.05)  =  95.1229%. 

The  reliability  of  the  system  can  be  calculated  as 
Rs  =  1  -  (1  -  Rc)^  =  1  -  (1  -  95. 12%)^=  99.9884%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  E) : 


Table  12 :  Agg 

fregated  Component  Level  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

1,493.82 

79 

2 

1,463.13 

75 

3 

1,492.59 

66 

4 

1,490.85 

70 

5 

1,490.58 

81 

Table  13:  Agg 

'regated  Component  : 

Level  OT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

373.26 

22 

2 

346.62 

15 

47 


Table  14:  Component  Level  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 
Component  1 

Xi  Failures 
Component  2 

Xi  Failures 
Component  3 

1 

497.94 

24 

35 

20 

2 

487.71 

26 

30 

19 

3 

497.53 

22 

21 

23 

4 

496.95 

20 

19 

31 

5 

496.86 

23 

28 

30 

Table  15:  Component 

:  Level  OT  Da' 

ta 

Run  i 

Ti  Hours 

BBiiiiiii 

Xi  Failures 
Component  3 

1 

124.42 

7 

8 

7 

2 

115.54 

3 

6 

6 

The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  From  the  data  above,  the 
following  results  are  obtained  (for  detailed  calculation, 
see  Appendix  F) : 
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Table  16:  System  ReliabilH 

:y  Using  Aggregated  Data 

Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.989% 

(99. 986%, 99. 991%) 
0.005% 

99.988% 

OT  R(l) 
80%  CI/PI 
Width 

99.987% 

(99. 978%, 99. 994%) 
0.016% 

99.988% 

(99.985%, 99.991%) 
0.006% 

99.988% 

DT  R{5) 
80%  CI/PI 
Width 

98.922% 

(98.720%, 99.100%) 
0.38% 

98.918% 

OT  R(5) 
80%  CI/PI 
Width 

98.837% 

(98.069%, 99.367%) 
1.298% 

98.906% 

(98.604%, 99.161%) 
0.557% 

98.918% 

DT  R{10) 
80%  CI/PI 
Width 

93.930% 

(92. 956%, 94.823%) 
1.867% 

93.910% 

OT  R(10) 
80%  CI/PI 
Width 

93.508% 

(89. 972%, 96.217%) 
6.245% 

93.865% 

(92.394%, 95.133%) 
2.739% 

93.910% 

DT  R(15) 
80%  CI/PI 
Width 

85.355% 

(83.335%, 87.259%) 
3.924% 

85.314% 

OT  R(15) 
80%  CI/PI 
Width 

84.480% 

(77.491%, 90.355%) 
12.864% 

85.238% 

(82. 202%, 87. 937%) 
5.735% 

85.314% 
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Table  17:  System  Reliability  Using  Component  Level  Data 


Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.989% 

(99. 984%, 99. 992%) 
0.008% 

99.988% 

or  R(l) 
80%  CI/PI 
Width 

99.988% 

(99.969%, 99.997%) 
0.028% 

99.988% 

(99.980%, 99.994%) 
0.014% 

99.988% 

DT  R(5) 
80%  CI/PI 
Width 

98.927% 

(98.567%, 99.222%) 
0.655% 

98.918% 

or  R(5) 
80%  CI/PI 
Width 

98.871% 

(97.412%, 99.641%) 
2.229% 

98.880% 

(98.225%, 99.369%) 
1.144% 

98.918% 

DT  R(10) 
80%  CI/PI 
Width 

93.962% 

(92.228%, 95.452%) 
3.224% 

93.910% 

or  R(10) 
80%  CI/PI 
Width 

93.699% 

(87.227%, 97.748%) 
10.521% 

93.772% 

(90.690%, 96.225%) 
5.535% 

93.910% 

DT  R(15) 
80%  CI/PI 
Width 

85.429% 

(81.881%, 88.644%) 
6.763% 

85.314% 

or  R(15) 
80%  CI/PI 
Width 

84.920% 

(72.587%, 93.991%) 
21.404% 

85.129% 

(78. 908%, 90.377%) 

1  11.469% 

85.314% 

Figure  8  shows  that  the  DT&E  classical  reliability 
estimates  slightly  overestimate  the  true  reliability.  The 
deviations  range  from  0.001%  to  0.12%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be 
insignificant.  The  OT&E  reliability  estimates  using  the 
classical  method  underestimated  the  true  reliability.  The 
deviations  range  from  0%  to  0.83%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be  not 
too  bad. 

The  OT&E  reliability  estimates  using  the  Bayesian 
method  are  all  very  close  (less  than  0.2%)  to  the  true 
reliability.  The  OT&E  reliability  estimates  using 
aggregated  component  level  data  underestimated  the  true 
reliability  with  deviations  ranging  from  0%  to  0.08%  for  the 
mission  times  being  considered.  Using  component  level 
data,  the  OT&E  reliability  estimates  also  underestimated  the 
true  reliability  with  deviations  ranging  from  0%  to  0.19% 
for  the  mission  times  being  considered.  The  OT&E 
reliability  estimates  using  the  Bayesian  method  are  closer 
to  the  true  reliability  than  the  classical  OT&E  reliability 
estimates. 
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Figure  9:  Comparison  of  OT&E  Interval  Widths 
(Aggregated  Level) 


Figure  10:  Comparison  of  OT&E  Interval  Widths 
(Component  Level) 
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Figures  9  and  10  again  show  the  strength  of  Bayesian 
analysis.  In  all  six  cases,  the  Bayesian  method  yields  much 
tighter  ranges  for  the  point  estimate  at  both  the  aggregated 
level  and  component  level.  At  the  aggregated  component 
level,  the  difference  between  the  OT&E  Bayesian  probability 
interval  width  and  OT&E  classical  confidence  interval  width 
ranges  from  0.01%  to  7.13%  for  the  mission  times  being 
considered.  At  the  component  level,  the  difference  ranges 
from  0.01%  to  9.94%.  Aggregated  component  level  reliability 
estimates  yield  tighter  intervals  than  component  levels 
aggregated  to  system.  For  the  mission  times  being 

considered,  the  difference  ranges  from  0.01%  to  8.54%  for 
OT&E  classical  confidence  interval  widths  and  from  0.01%  to 
5.73%  for  OT&E  Bayesian  probability  interval  widths. 
Reliability  Analysis  of  a  Parallel  System  with  Different 
Independen t  Componen ts . 

Consider  a  system  composed  of  3  different  independent 
components  arranged  in  a  parallel  configuration.  The  time 
to  failure  of  each  component  is  assumed  to  be  exponentially 
distributed.  Let  A,i  =  0.1  failure  per  hour,  X2  =  0.05 
failure  per  hour,  and  X3  =  0.025  failure  per  hour  be  the 
respective  constant  failure  rates  of  each  component. 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  and  3  are 
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Ri  =  exp(-O.lOO)  =  90.48% 

Rz  =  exp (-0.050)  =  95.12% 

Rs  =  exp(-0.025)  =  97.53%. 

The  reliability  of  the  system  can  be  calculated  as 
R,  =  1  -  (1  -  Ri)  (1  -  Rz)  (1  -  R3) 

=  1  -  (1  -  0.9048)  (1  -  0.9512) (1  -  0.9753)  =  99.9885%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  G) : 


Table  18:  Component 

Level  DT  Da' 

ta 

Run  i 

Ti  Hours 

Xi  Failures 
Component  1 

Xi  Failures 
Component  2 

Xi  Failures 
Component  3 

1 

488.03 

55 

27 

10 

2 

498.15 

53 

30 

8 

3 

498.39 

44 

24 

12 

4 

496.55 

36 

26 

9 

5 

496.36 

60 

33 

9 

Table  19:  Component  Level  OT  Da' 

ta 

Run  i 

Ti  Hours 

Xi  Failures 
Component  1 

Xi  Failures 
Component  2 

1 

122.22 

13 

9 

3 

2 

113.29 

8 

8 

4 

The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  From  the  data  above,  the 
following  results  are  obtained  (for  detailed  calculation, 
see  Appendix  H) : 
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Table  20:  System  Reliability  Using  Component  Level  Data 


Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.99% 

(99. 986%, 99. 993%) 
0.007% 

99.989% 

OT  R(l) 
80%  CI/PI 
Width 

99.983% 

(99.957%, 99.995%) 
0.038% 

99.989% 

(99. 983%, 99. 994%) 
0.011% 

99.989% 

DT  R(5) 
80%  CI/PI 
Width 

99.106% 

(98.776%, 99.371%) 
0.595% 

98.977% 

OT  R(5) 
80%  CI/PI 
Width 

98.495% 

(96.682%, 99.503%) 
2.821% 

99.062% 

(98.584%, 99.428%) 
0.844% 

98.977% 

DT  R(10) 
80%  CI/PI 
Width 

95.189% 

(93.661%, 96.496%) 
2.835% 

94.498% 

OT  R(I0) 
80%  CI/PI 
Width 

92.202% 

(84.868%, 97.087%) 
12.219% 

94.974% 

(92.813%, 96.722%) 
3.909% 

94.498% 

DT  R(15) 
80%  CI/PI 
Width 

88.792% 

(85.7%, 91.578%) 
5.878% 

87.183% 

OT  R(15) 
80%  CI/PI 
Width 

82.457% 

(69.464%, 92.678%) 
23.214% 

88.323% 

(84. 066%, 91. 935%) 
7.869% 

87.183% 

Figure  11  shows  that  the  DT&E  classical  reliability 
estimates  overestimated  the  true  reliability.  The  deviations 
range  from  0.001%  to  1.61%  for  the  mission  times  being 
considered.  The  deviations  are  considered  to  be 
significant.  The  OT&E  reliability  estimates  using  the 
classical  method  underestimated  the  true  reliability.  The 
deviations  range  from  0.006%  to  5.87%  for  the  mission  times 
being  considered.  The  deviations  are  again  considered  to  be 
significant. 
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Using  component  level  data,  the  OT&E  reliability 
estimates  overestimated  the  true  reliability  with  deviations 
ranging  from  0%  to  1.14%  for  the  mission  times  being 
considered.  While  the  deviations  are  still  large,  the  OT&E 
reliability  estimates  using  the  Bayesian  method  are  much 
closer  to  the  true  reliability  than  the  classical  OT&E 
reliability  estimates. 


Figure  11:  Classical  and  Bayesian  Reliability  Estimates 


Classical  &  Bayesian  Reliability  Estimations 


Figure  12  once  again  shows  the  strength  of  Bayesian 
analysis.  In  all  six  cases,  the  Bayesian  method  yields  much 
tighter  ranges  for  the  point  estimate.  The  difference 
between  the  OT&E  Bayesian  probability  interval  width  and 
OT&E  classical  confidence  interval  width  ranges  from  0.03% 
to  15.35%  for  the  mission  times  being  considered. 
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Figure  12:  Comparison  of  OT&E  Interval  Widths 


Comparison  of  OT&E  Interval  Widths 


OT  Connponent  Classical 
X  OT  ConnponerTt  Bayesian 
O  Oorrponent  Different 


0  5  10  15  20 

Mfesion  Time 

Reliability  Analysis  of  a  Bridge  System  with  Independent  and 
Identical  Components. 

Consider  a  system  composed  of  5  independent  and 
identical  components  given  in  Figure  1  (page  30) .  The  time 
to  failure  of  each  component  is  assumed  to  be  exponentially 
distributed.  Let  X  =  0.05  failure  per  hour  be  the  constant 
failure  rate  of  each  component.  For  the  true  system,  using 
a  mission  time,  t  of  one  hour,  the  reliability  of  components 
1,  2,  3,  4,  and  5  are 

Rc  =  Ri  =  R2  =  Rs  =  R4  =  Rs  =  exp  (-0.05)  =  95.12%. 

The  reliability  of  the  system  can  be  calculated  as 

Rs  =  R1R4  t  R2R5  t  R1R3R5  +  R2R3R4  ~  R1R2R3R4  “  R1R2R3R5  ~ 
R1R2R4R5  —  R1R3R4R5  —  R2R3R4R5  t  2R1R2R3R4R5  ~  99.50%. 
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Using,  a  computer  program  called  RAPTOR  to  simulate  the 
true  system,  the  following  data  are  obtained  (to  review  the 


simulated  data,  see  Appendix  I) : 


Table  21 


Run  i 


Ti  Hours 

Xi  Failures 

2,496.85 

136 

2,486.35 

123 

2461.30 

124 

2493.05 

125 

2498.40 

130 

Table  22 :  Aggregated  Component  Level  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

618.85 

30 

2 

604.10 

26 

Table  23:  Component  Level  DT  Data 


Run  i 

Ti  Hours 

C2  Xi 

C5  Xi 

1 

499.37 

32 

21 

35 

25 

23 

2 

497.27 

28 

17 

24 

27 

27 

3 

492.26 

27 

28 

25 

20 

24 

4 

498.61 

26 

27 

22 

23 

27 

5 

499.68 

27 

23 

33 

25 

22 

Table  24:  Component  Level  OT 

Data 

Run  i 

Ti  Hours 

Cl  Xi 

C2  Xi 

C3  Xi 

C4  Xi 

C5  Xi 

1 

123.77 

7 

6 

8 

3 

6 

2 

120.82 

6 

5 

6 

4 

5 

The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  From  the  data  above,  the 
following  results  are  obtained  (for  detailed  calculation, 
see  Appendix  J) : 
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Table  25;  System  Reliability  Using  Aggregated  Data 


Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.48% 

(99.42%, 99.53%) 
0.11% 

99.50% 

OT  R(l) 
80%  CI/PI 
Width 

99.58% 

(99.43%, 99.71%) 
0.28% 

99.49% 

(99.44%, 99.54%) 
0.1% 

99.50% 

DT  R(4) 
80%  CI/PI 
Width 

92.39% 

(91.67%, 93.07%) 
1.4% 

92.74% 

OT  R(4) 
80%  CI/PI 
Width 

93.82% 

(91.70%, 95.62%) 
3.92% 

92.52% 

(91. 85%, 93. 16%) 
1.31% 

92.74% 

DT  R(7) 
80%  CI/PI 
Width 

79.95% 

(78.29%, 81.55%) 
3.26% 

80.76% 

OT  R{7) 
80%  CI/PI 
Width 

83.36% 

(78.37%, 87.80%) 
9.43% 

80.25% 

(78. 70%, 81. 75%) 
3.05% 

80.76% 

DT  R(10) 
80%  CI/PI 
Width 

65.75% 

(63.35%, 68.13%) 
4.78% 

66.95% 

OT  R(10) 
80%  CI/PI 
Width 

70.86% 

(63. 47%, 77. 91%) 
14.44% 

66.21% 

(63.94%, 68.43%) 
4.49%  : 

66.95% 

59 


Table  26:  System  Reliability  Using  Component  Level  Data 


Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.50% 

(99. 38%, 99. 61%) 
0.23% 

99.50% 

or  R(l) 
80%  CI/PI 
Width 

99.63% 

(99.24%, 99.85%) 
0.61% 

99.52% 

(99.38%, 99.65%) 
0.27% 

99.50% 

DT  R(4) 
80%  CI/PI 
Width 

92.64% 

(91.01%, 94.10%) 
3.09% 

92.74% 

or  R(4) 
80%  CI/PI 
Width 

94% 

(89.18%, 97.66%) 
8.48% 

92.87% 

(90.91%, 94.67%) 
3.76% 

92.74% 

DT  R(7) 
80%  CI/PI 
Width 

80.48% 

(76.76%, 83. 98%) 
7.22% 

80.76% 

or  R{7) 
80%  CI/PI 
Width 

84.65% 

(72.73%, 93.18%) 
20.45% 

81.05% 

(76.50%, 85.40%) 
8.90% 

80.76% 

DT  R(10) 
80%  CI/PI 
Width 

66.52% 

(61. 16%, 71. 79%) 
10.63% 

66.95% 

or  R(i0) 
80%  CI/PI 
Width 

72.81% 

(55.66%, 87.02%) 
31.36% 

67.39% 

(60. 78%, 74. 03%) 
13.25% 

66.95% 

Figure  13  shows  that  the  DT&E  classical  reliability 
estimates  underestimated  the  true  reliability.  The 
deviations  range  from  0.02%  to  1.2%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be 
significant.  The  OT&E  reliability  estimates  using  the 
classical  method  overestimated  the  true  reliability.  The 
deviations  range  from  0.08%  to  5.86%  for  the  mission  times 
being  considered.  The  deviations  are  again  considered  to  be 
significant . 

The  OT&E  reliability  estimates  using  the  Bayesian 


method  are 

all 

very 

close 

(less  than 

1% )  to 

the 

true 

reliability. 

The 

OT&E 

reliability 

estimates 

using 

aggregated 

component 

level 

data  underestimated 

the 

true 

reliability 

with 

deviations 

ranging  from 

0.01%  to 

0.74%  for 

the  mission  times  being  considered.  Using  component  level 
data,  the  OT&E  reliability  estimates  overestimated  the  true 
reliability  with  deviations  ranging  from  0.02%  to  0.44%  for 
the  mission  times  being  considered.  While  the  deviations 
are  still  significant,  the  OT&E  reliability  estimates  using 
the  Bayesian  method  are  much  closer  to  the  true  reliability 
than  the  classical  OT&E  reliability  estimates. 
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Figure  14:  Comparison  of  OT&E  Interval  Widths 
(Aggregated  Level) 


Comparison  of  OT&E  Interval  Widths 


X  OT  Aggregated  Classical 
X  OT  Aggregated  Bayesian 
O  Aggregated  Dfferent 


Figure  15:  Comparison  of  OT&E  Interval  Widths 
(Component  Level) 
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Figures  14  and  15  again  show  the  strength  of  Bayesian 
analysis.  In  all  seven  cases,  the  Bayesian  method  yields 
much  tighter  ranges  for  the  point  estimate  at  both  the 
aggregated  level  and  component  level.  At  the  aggregated 
component  level,  the  difference  between  the  OT&E  Bayesian 
probability  interval  width  and  OT&E  classical  confidence 
interval  width  ranges  from  0.18%  to  9.95%  for  the  mission 
times  being  considered.  At  the  component  level,  the 
difference  ranges  from  0.34%  to  18.11%.  Aggregated 
component  level  reliability  estimates  yield  tighter 
intervals  than  component  level  aggregated  to  system.  For 
the  mission  times  being  considered,  the  difference  ranges 
from  0.33%  to  16.92%  for  OT&E  classical  confidence  interval 
widths  and  from  0.17%  to  8.76%  for  OT&E  Bayesian  probability 
interval  widths. 

Reliability  Analysis  of  a  Bridge  System  with  Different 
Independen t  Componen ts , 

Consider  a  system  composed  of  5  independent  components 
given  in  Figure  1  (page  30)  .  The  time  to  failure  of  each 
component  is  assumed  to  be  exponentially  distributed.  Let 
=  0.1  failure  per  hour,  Xz  =  0.05  failure  per  hour,  Xz  = 
0.033  failure  per  hour,  A,4  =  0.025  failure  per  hour,  and  A,5 

=  0.02  failure  per  hour  be  the  respective  constant  failure 
rates  of  each  component. 
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For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  3,  4,  and  5  are 


Ri  = 

=  exp( 

-0.10) 

= 

90.48% 

R2  = 

=  exp( 

-0.05) 

= 

95.12% 

R3  = 

=  exp( 

-0.033) 

=: 

96.75% 

R4  = 

=  exp( 

-0.025) 

= 

97.53% 

R5  = 

=  exp( 

-0.20) 

= 

98.02%. 

The 

reliability 

Of 

the  system  can  be  calculated  as 

Rs 

-  R1R4 

+  R2R5 

+ 

R1R3R5  +  R2R3R4  —  RaR2R3R4  “  R1R2R3R5 

R1R2R4R5  ~  R1R3R4R5  ~  R2R3R4R5  2R1R2R3R4R5  ~  99.48%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  K) : 


Table  27 ;  Component  Level  DT  Data 


Run  i 

Ti  Hours 

Cl  Xi 

C2  Xi 

C3  Xi 

C4  Xi 

C5  Xi 

1 

499.09 

51 

19 

22 

15 

13 

2 

490.55 

41 

23 

21 

12 

9 

3 

497.97 

44 

30 

17 

6 

8 

4 

495.31 

52 

26 

18 

12 

10 

5 

498.76 

48 

25 

22 

18 

15 

Table  28:  Component  Level  OT  Data 
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The  reliability  of  the  system  for  a  mission  time  t  will 
be  estimated  and  80%  confidence  intervals  associated  with 
the  estimated  system  reliability  will  be  computed  using  the 
methodology  given  in  Chapter  III.  From  the  data  above,  the 
following  results  are  obtained  (for  detailed  calculation, 
see  Appendix  L) : 


Table  29:  System  Reliability  Using  Component  Level  Data 


Classical  Method 

Bayesian  Method 

True  R 

DT  R(l) 
80%  CI/PI 
Width 

99.50% 

(99.38%, 99.59%) 
0.21% 

99.48% 

OT  R(l) 
80%  CI/PI 
Width 

99.51% 

(99.14%, 99.76%) 
0.62% 

99.50% 

(99.33%, 99.65%) 
0.32% 

99.48% 

DT  R(4) 
80%  CI/PI 
Width 

93.08% 

(91.64%, 94.37%) 
2.73% 

92.97% 

OT  R(4) 
80%  CI/PI 
Width 

93.52% 

(89. 08%, 96. 70%) 
7.62% 

93.24% 

(91. 06%, 95. 12%) 
4.06% 

92.97% 

DT  R(7) 
80%  CI/PI 
Width 

82.21% 

(78. 92%, 85. 27%) 
6.35% 

82.13% 

OT  R(7) 
80%  CI/PI 
Width 

83.58% 

(73.88%, 91.23%) 
17.35% 

82.71% 

(77.73%, 87.16%) 
9.43% 

82.13% 

DT  R(10) 
80%  CI/PI 
Width 

69.82% 

(64. 97%, 74. 52%) 
9.55% 

i 

69.86% 

OT  R(10) 
80%  CI/PI 
Width 

72.29% 

(58.41%, 84.44%) 
26.03% 

70.71% 

(63. 34%, 77. 61%) 
14.27% 

69.86% 

Figure  16:  Classical  and  Bayesian  Reliability  Estimates 
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Figure  16  shows  that  the  DT&E  classical  reliability 
estimates  are  very  close  to  the  true  reliability.  The 
deviations  range  from  -0.04%  to  0.11%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be 
insignificant.  The  OT&E  reliability  estimates  using  the 
classical  method  overestimated  the  true  reliability.  The 
deviations  range  from  0.03%  to  2.43%  for  the  mission  times 
being  considered.  The  deviations  are  considered  to  be 
significant . 

Using  component  level  data,  the  OT&E  reliability 
estimates  overestimated  the  true  reliability  with  deviations 
ranging  from  0.02%  to  0.85%  for  the  mission  times  being 
considered.  While  the  deviations  are  still  large,  the  OT&E 
reliability  estimates  using  the  Bayesian  method  are  much 
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closer  to  the  true  reliability  than  the  classical  OT&E 
reliability  estimates. 

Figure  17:  Comparison  of  OT&E  Interval  Widths 

Comparison  of  OT&E  Interval  Widths 
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Figure  17  once  again  shows  the  strength  of  Bayesian 
analysis.  In  all  seven  cases,  the  Bayesian  method  yields 
much  tighter  ranges  for  the  point  estimate.  The  difference 
between  the  OT&E  Bayesian  probability  interval  width  and 
OT&E  classical  confidence  interval  width  ranges  from  0.3%  to 
11.76%  for  the  mission  times  being  considered. 

Summary. 

The  overall  trend  from  the  results  obtained  in  this 
chapter  is  that  the  Bayesian  methods  used  in  this  research 
produced  much  shorter  confidence  intervals  compared  to  the 
classical  methods.  Another  trend  is  that  system  level  data 
and  aggregated  component  level  data  yield  tighter  confidence 
intervals  than  component  level  data  aggregated  to  system. 
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For  the  simulated  systems  considered  in  this  research,  the 
Bayesian  reliability  point  estimates  are  much  closer  to  the 
true  reliability  of  the  system  as  compared  to  the  classical 
reliability  estimates. 
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V .  Conclus ions 


The  results  in  Chapter  IV  show  that  the  Bayesian 
methods  used  in  this  research  produced  much  shorter 
confidence  intervals  about  the  reliability  point  estimates 
for  a  variety  of  systems  than  the  traditional  classical 
methods.  This  research  shows  that  meaningful  results  (in 
the  form  of  shorter  confidence  intervals)  can  be  obtained  by 
simply  applying  Bayesian  techniques  to  consecutive  test 
phases. 

The  binomial  model  shows  that  while  the  OT&E  point 
estimates  are  similar  (90.63%  vs.  91.67%),  the  Bayesian 
method  produced  a  much  tighter  confidence  interval  (12.66% 
vs.  27.88%),  thus  making  it  easier  for  managerial 
decision-making.  The  exponential  model  applied  to  a  real 
set  of  data  again  shows  that  the  Bayesian  method  produced 
smaller  confidence  intervals. 

For  a  series  system  with  exponential  components,  using 
simulated  data  for  a  3-component  system,  the  Bayesian  method 
produced  tighter  confidence  intervals.  The  reliability 
point  estimates  using  the  Bayesian  method  are  closer  to  the 
true  reliability  of  the  system  as  compared  to  the  classical 
reliability  point  estimates.  System  level  reliability 
estimates  yield  much  tighter  confidence  intervals  than 
component  level  aggregated  to  system.  From  the  results  in 
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Chapter  IV,  note  that  the  Bayesian  Method  1  only  produced 
the  reliability  point  estimates.  Thus  this  method  might  not 
useful  for  a  decision  maker. 

For  a  parallel  system  with  independent  and  identical 
components,  using  simulated  data  for  a  3-component  system, 
the  Bayesian  method  again  produced  tighter  confidence 
intervals.  The  reliability  point  estimates  using  the 
Bayesian  method  are  again  closer  to  the  true  reliability  of 
the  system  as  compared  to  the  classical  reliability  point 
estimates.  Aggregated  component  level  reliability  estimates 
yield  much  tighter  confidence  intervals  than  component  level 
aggregated  to  system. 

For  a  parallel  system  with  different  independent 
components,  using  simulated  data  for  a  3-component  system, 
the  Bayesian  method  once  again  produced  tighter  confidence 
intervals.  The  reliability  point  estimates  using  the 
Bayesian  method  are  again  closer  to  the  true  reliability  of 
the  system  as  compared  to  the  classical  reliability  point 
estimates. 

For  a  bridge  system  given  in  Figure  1  (page  30)  with 
independent  and  identical  components,  using  simulated  data, 
the  Bayesian  method  again  produced  much  tighter  confidence 
intervals.  The  reliability  point  estimates  using  the 
Bayesian  method  are  closer  to  the  true  reliability  of  the 
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system  as  compared  to  the  classical  reliability  point 
estimates.  Aggregated  component  level  reliability  estimates 
yield  a  much  tighter  confidence  intervals  than  component 
level  aggregated  to  system. 

For  a  bridge  system  given  in  Figure  1  (page  30)  with 
different  independent  components,  using  simulated  data,  the 
Bayesian  method  once  again  produced  much  tighter  confidence 
intervals.  The  reliability  point  estimates  using  the 
Bayesian  method  are  much  closer  to  the  true  reliability  of 
the  system  as  compared  to  the  classical  reliability  point 
estimates. 

In  summary,  the  results  of  this  research  show  that  the 
Bayesian  method  produced  a  much  tighter  confidence  interval 
as  compared  to  the  classical  method  for  OT&E  data.  System 
level  data  and  aggregated  component  level  data  will  yield  a 
tighter  confidence  interval  than  component  level  data 
aggregated  to  system.  For  multiphase  test  process  commonly 
implemented  by  agencies  such  as  AFOTEC,  the  Bayesian  method 
will  work  with  using  all  the  previous  test  data  as  prior 
information  for  the  current  test.  However,  it  is  up  to  the 
analysts  to  determine  if  the  previous  test  data  are 
representative  of  the  current  test.  Otherwise,  you  are 
comparing  apples  to  oranges  and  the  results  obtained  are 
meaningless . 
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J^pendix  A:  Calcula'kion  For  The  Binomial  Model. 

The  Classical  Binomial  Method. 

Suppose  the  following  data  are  obtained  from  a  missile 
program.  For  the  DT&E  portion  of  the  test,  two  missiles 
failed  out  of  20  missiles  launched.  For  the  OT&E  portion  of 
the  test,  one  missile  failed  out  of  12  missiles  launched. 
Assuming  the  data  have  the  binomial  distribution,  point 
estimates  for  the  probability  of  success  p  and  80% 
confidence  interval  about  p  will  be  computed.  The  point 
estimate  of  p  for  DT&E  and  OT&E  are 

DT&E:  p  =  (^)  =  90.00% 

OT&E:  p  =  (-^)  =  91.67%. 


The  80%  TCI  are 
DT&E;  (75. 52%, 97. 31%) 


Lower  bound: 


18 


0 


18+(20-18  +  1)Fi-o.2/2(2x20-2x18  +  2,  2x18), 

,  36))  ^  (tS+^xIT^)  ""  0-7552 


18-I-3Fo.9(6 

_  (  (18  +  l)Fi-o. 2/2(2x18  +  2,  2x20-2x18)  ^ 

Upper  bound:(^  (20-18)+(18  +  l)Fi-o. 2/2(2x18  +  2,  2x20-2x18) j 

=  f  ■  19^^0.9(38,4)  ^  ^  f  19x3.81  ^  ^  ^  9-3- 
V  2  +  19Fo.9(38,  4)j  V  2  +  19x3. 81  j 

OT&E;  (71. 25%, 99. 13%) 

[  11+(12-11  +  1)Fi-o.2/2(2x12-2x11  +  2, 

(  ll  +  2Fo.l(4,  22))  ""  (  11  +  2^x2.22)  ""  0.7125 


Lower  bound: 


:,  2x11)) 
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(  (ll  +  l)Fi-o. 2/2(2x11  +  2, 2x12-2x11)  ^ 

Upper  bound:j^  (12-ll)+(ll  +  l)Fi-o.2/2(2xll  +  2,  2x12-2x11) j 

^  (  12Fo.s(24,2)  ^  ^  (-_12X9.45  ^  ^ 

V  1  +  12Fo.9(24,  2)j  V  1  +  12x9. 45^  «•»»■«•■» 

The  Bayesian  Binomial  Method. 

There  are  two  failures  out  of  20  missiles  launched 

during  DT&E,  thus  Xo  =  18  and  no  =  20  and  the  B(18,20)  prior 
distribution  is  used.  There  is  one  failure  out  of  12 
missiles  launched  during  OT&E,  thus  x  =  11  and  n  =  12. 

The  mean  of  the  posterior  distribution  is  the  Bayesian 
point  estimator.  Under  the  squared  error  loss  function,  the 
Bayesian  point  estimator  is 

P  =  E(p|x;xo,no)  =  (Irfr)  =  (12+20)  =  ( H)  =  50. 63%. 
For  the  80%  TBPI,  the  upper  and  lower  interval 
endpoints  are  (83.73%, 96.39%) . 

Lower  bound: 


11  +  18 


11  +  18+(12+20-11-18)Fi_o.2/2(2x12+2x20-2x  11-2x18,  2x11+2x18) 


) 


f  - .  29  )  1 

(  2  9  ^ 

V  29+3Fo.9(6,  58)j  ' 

k  29  +  3x1. 88>^ 

)  -  0. 


8373 


Upper  bound: 


( 


_ (11  +  18)Fi-o.2/2(2x  11+2x18,  2x12+2x20-2x11-2x18) _ ^ 

12+20-11-18+(11  +  18)Fi_o.2/2(2x  11+2x18,  2x12+2x20-2x11-2x18); 


=  f  29Fo.9(58,  6)  ^  f  29x2. IG  \  _  . 

3+29Fo.9(58,  6)j  V  3  +  29x2. 76>^  “  ”* 


9639 
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Appendix  B:  Calculation  For  The  Exponential  Model. 

The  Classical  Exponential  Method. 

The  following  data  are  obtained  from  a  reliability  test 
program. 


Table  Bl:  OT&E  Data  (3;6) 


Terminal  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

184.3 

12 

0.06511 

2 

232.6 

9 

0.03869 

3 

172.8 

9 

0.05208 

4 

284.4 

17 

0.05977 

5 

264.8 

18 

0.06798 

Table  B2:  OA  Data  (3:7) 


Terminal  i 

Ti  Hours 

Xi  Failures 

1 

68.5 

2 

2 

60.9 

2 

Assuming  the  time  to  failure  of  the  terminals  have  the 
exponential  distribution,  point  estimates  for  the  failure 
rate  X  and  80%  confidence  interval  about  X  will  be  computed. 


OT&E: 

1  = 

‘  ( 1138.9) 

=  0.05707 

OA: 

X  = 

It 

■  (  129.4) 

=  0.03091 

where  = 

65 

is  the 

number  of 

failures 

and  Ti  =  1138.9  is 

the  total 

test 

time 

from  Table 

Bl  and  fz 

=  4  is  the  number 

of  failures  and  T2  =  129.4  is  the  total  test  time  from  Table 
B2. 

The  80%  two-sided  confidence  intervals  (TCI)  are 

f  Xl2/2(2X65)  Xl-o. 2/2(2x65)  \  (  xli(130)  9(130)  ^ 

OT&E:  ^  2x1138.9  '  2x1138.9  7  =  V  2211  •  2277.8  ) 
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=  (2277^'  227^)  =  (0.04821,0.06631) 

(  ;ifo.2/2(2x4)  yi-0.2/2(2x4)^  f  Xo.li^)  /o.9(8)^ 

QA:  y  2x129.4  '  2x129.4  J  ^  \  258.8  '  258.8  y 

f  3.49  13. 36^  ,, 

“  1^  258.8  '  258.8  J  ~  (0.01349,0.05163) 

The  reliability  estimators  of  R(t)  for  a  mission  time  t 
=  1  hour  are 

OT&E;  R{t)  =  94.45% 

R{t]  =  exp(-At)  =  exp{-0. 05707X1)  =  0.9445 

OA:  R(t)  =  96.96% 

R{t]  =.  exp(-2t)  =  exp(-0. 03091X1)  =  0.9696. 

The  80%  two-sided  confidence  intervals  (TCI)  are 

OT&E:  (93. 58%, 95. 29%) 

exp(-0. 06631X1)  R{1)  <  exp  (-0 . 04821X1 ) 

OA:  (94. 97%, 98. 66%) 

exp(-0. 05163X1)  ^  R{1}  <  exp(-0. 01349X1)  . 

The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior.  OT&E  data  are  used  to 
obtain  a  gamma  prior.  The  mean  and  variance  of  the  gamma 
distribution  are  E(^)  =  a/p  and  V(X,)  =a/p^.  E(A,)  and  V(X) 

can  be  estimated  using  the  sample  mean  and  variance  of  X 
respectively.  Thus  a  and  P  can  be  computed  as  a  = 
E^{X)/V(X)  and  p  =  E{X)/V{X). 
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From  Table  Bl: 


E{1)  =  X  =  0.05673  and  V{X)  =  =  0.00014. 


These  yielded  the  parameter  estimates 


a 


0.05673^ 

0.00014 


=  23.28  and 


f  0.05673^ 
0.000147 


405.21. 


Bayesian  Exponential  Posterior  Analysis.  The  mean  of 
the  posterior  distribution  is  the  Bayesian  point  estimator. 
Under  the  squared  error  loss  function,  the  Bayesian  point 
estimator  is 

1  =  E(X|f;a,P)  =  =  (  12 +^40^5^.21)  "  0*05103. 

The  80%  two-sided  probability  interval  (TBPI)  is 

f  . 28)  \ 

2(129.4+405.21)  '  2(129.4+405.21)  ) 

(  41.18  67_:_0_^ 

“  ^  1069.22  '  1069. 22J  “  (0.03852,0.06329). 


The  Bayesian  estimator  of  R(t)  for  t  =  1  hour  is 


95.03% 


The  80%  two-sided  probability  interval  (TBPI)  is 
exp(-0. 06329X1)  <  R(l)  <  exp  (-0. 03852X1) 

(93.87%  <  R(l)  <  96.22%)  . 
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Appendix  C.  Simulated  Data  For  Series  System. 


DT  Data  Run  #1 


Time 

Component 

Failure 

Time 

Component 

Failure 

4.06 

C3 

244.08 

Cl 

6.20 

C2 

245.78 

C3 

12.19 

C2 

247.35 

C3 

15.64 

Cl 

262.18 

C2 

16.39 

Cl 

278.38 

Cl 

22.66 

C2 

279.89 

Cl 

39.58 

Cl 

280.96 

C3 

43.63 

C3 

282.47 

C2 

44.47 

C3 

288.81 

C2 

49.98 

C2 

296.84 

C3 

53.87 

C2 

308.01 

Cl 

57.15 

Cl 

308.60 

Cl 

71.94 

Cl 

310.25 

C2 

72.13 

Cl 

316.96 

C3 

73.51 

C3 

329.51 

Cl 

77.89 

C3 

333.17 

C2 

82.85 

C2 

336.92 

C3 

91.34 

Cl 

340.08 

C2 

93.36 

C3 

343.54 

Cl 

101.92 

C3 

362.98 

C2 

106.49 

C2 

369.64 

C2 

113.87 

Cl 

373.05 

C2 

114.93 

C2 

379.03 

C2 

127.01 

Cl 

395.63 

C3 

139.00 

C3 

396.23 

C3 

141.46 

C2 

400.42 

Cl 

145.83 

Cl 

407.82 

C2 

149.97 

C3 

411.09 

C2 

151.10 

C2 

415.78 

C2 

169.35 

C2 

421.76 

C3 

178.76 

Cl 

426.12 

Cl 

197.32 

C3 

430.44 

Cl 

197.58 

C2 

433.82 

Cl 

198.19 

C3 

444.57 

C3 

199.98 

C3 

461.05 

C2 

200.04 

Cl 

485.22 

C2 

226.10 

C3 

495.04 

C3 

227.85 

C2 

496.78 

C3 

228.57 

C3 

497.25 

Cl 

239.47 

C3 

499.86 

C3 
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DT  Data  Run  #2 


Time 

Component 

Failure 

Time 

Component 

Failure 

12.04 

C3 

227.72 

C3 

12.07 

Cl 

255.18 

C3 

18.71 

C3 

255.70 

C2 

24.93 

C2 

261.48 

C3 

31.60 

C2 

267.05 

C2 

32.47 

C3 

274.55 

Cl 

32.66 

C3 

276.86 

Cl 

37.16 

Cl 

281.57 

Cl 

37.68 

C3 

281.73 

C2 

41.33 

C3 

286.19 

C3 

53.45 

C2 

286.56 

Cl 

64.13 

Cl 

296.05 

C3 

73.42 

Cl 

296.35 

C2 

75.56 

Cl 

298.56 

C2 

87.54 

C2 

301.29 

C3 

87.78 

C2 

301.38 

Cl 

88.24 

C2 

308.87 

Cl 

93.96 

C2 

313.78 

C2 

97.05 

C2 

318.71 

Cl 

99.47 

C2 

323.07 

C3 

135.04 

C2 

323.64 

C3 

139.54 

C2 

336.32 

Cl 

154.84 

Cl 

337.83 

Cl 

159.06 

C3 

351.36 

Cl 

161.17 

C2 

353.18 

Cl 

165.41 

Cl 

353.51 

C2 

173.80 

C2 

363.70 

C3 

183.58 

C2 

366.73 

C2 

185.59 

C2 

374.64 

C2 

187.10 

C3 

376.48 

C2 

189.16 

C2 

388.35 

C2 

195.88 

C3 

396.67 

C2 

198.87 

C3 

397.31 

Cl 

199.44 

C2 

398.38 

C2 

204.14 

C2 

398.81 

C2 

206.88 

C2 

403.77 

C2 

207.29 

C2 

404.07 

C3 

224.62 

Cl 

410.53 

C3 

225.61 

C2 

449.53 

C2 

DT  Data  Run  #3 


Time 

Component 

Failure 

Time 

Component 

Failure 

12.73 

C2 

171.78 

Cl 

16.97 

Cl 

176.02 

C3 

21.21 

C3 

207.53 

C2 

51.23 

C2 

220.65 

Cl 

55.53 

C2 

224.89 

C3 

68.30 

Cl 

225.93 

C2 

12. bA 

C3 

272.35 

C2 

13. IQ 

C2 

276.70 

Cl 

74.04 

Cl 

280.94 

C3 

78.29 

C3 

285.39 

C2 

98.37 

Cl 

301.25 

Cl 

101.92 

C2 

305.49 

C3 

102.61 

C3 

308.09 

C2 

103.04 

C2 

324.04 

C2 

119.60 

C2 

330.31 

C2 

128.83 

C2 

330,52 

C2 

135.90 

Cl 

348,70 

C2 

137.39 

Cl 

363.13 

Cl 

140.14 

C3 

367.37 

C3 

141.63 

C3 

372.01 

C2 

159.47 

Cl 

377.28 

C2 

163.71 

C3 

380.51 

Cl 

165.49 

C2 

384.76 

C3 

79 


DT  Data  Run  #4 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

2.18 

C2 

137.72 

C3 

334.87 

C2 

2.39 

C2 

142.78 

C2 

338.17 

C2 

2.91 

Cl 

155.78 

Cl 

366.00 

Cl 

3.12 

C3 

155.98 

C3 

366.21 

C3 

3.18 

Cl 

187.69 

C2 

400.51 

C2 

3.39 

C3 

190.38 

Cl 

405.11 

C2 

15.20 

C2 

190.58 

C3 

406.97 

C2 

20.27 

Cl 

201.10 

C2 

407.45 

C2 

20.95 

C2 

206.01 

C2 

421.23 

C2 

27.89 

C2 

214.66 

C2 

432.10 

C2 

27.93 

Cl  _ 

224.39 

C2 

435.24 

C2 

28.13 

C3 

227.45 

C2 

435.88 

Cl 

37.18 

Cl 

242.22 

C2 

436.09 

C3 

37.39 

C3 

248.19 

C2 

446.39 

C2 

38.87 

C2 

250.25 

Cl 

446.49 

Cl 

42.58 

C2 

250.45 

C3 

446.69 

C3 

51.82 

Cl 

268.14 

Cl 

450.89 

Cl 

52.03 

C3 

268.34 

C3 

451.09 

C3 

53.20 

C2 

274.50 

C2 

461.14 

C2 

56.77 

Cl 

274.68 

Cl 

473.75 

C2 

56.98 

C3 

274.88 

C3 

475.42 

C2 

62.50 

C2 

286.21 

Cl 

517.25 

C2 

70.94 

Cl 

286.41 

C3 

532.99 

C2 

71.14 

C3 

299.19 

Cl 

534.02 

Cl 

82.07 

C2 

299.39 

C3 

534.22 

C3 

83.33 

Cl 

303.27 

Cl 

534.39 

C2 

83.54 

C3 

303.47 

C3 

540.14 

Cl 

103.14 

C2 

322.96 

Cl 

540.35 

C3 

109.42 

Cl 

323.17 

C3 

542.63 

Cl 

109.63 

C3 

326.91 

C2 

542.83 

C3 

116.83 

C2 

330.92 

Cl 

543.27 

Cl 

137.52 

Cl 

331.12 

C3 

543.47 

C3 
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DT  Data  Run  #5 


Time 

Component 

Failure 

2.04 

C2 

2.72 

Cl 

3.40 

C3 

7.41 

C2 

9.88 

Cl 

10.56 

C3 

27.19 

C2 

36.25 

Cl 

36.93 

C3 

133.50 

C2 

145.89 

C2 

146.74 

C2 

148.39 

C2 

178.00 

Cl 

178.68 

C3 

188.27 

C2 

194.52 

Cl 

195.20 

C3 

195.65 

Cl 

196.33 

C3 

197.85 

Cl 

198.53 

C3 

234.22 

C2 

234.27 

C2 

243.58 

C2 

251.03 

Cl 

251.71 

C3 

284.65 

C2 

297.02 

C2 

312.30 

Cl 

312.37 

Cl 

312.98 

C3 

313.04 

C3 

324.77 

Cl 

325.45 

C3 

325.61 

C2 

327.78 

C2 

342.93 

C2 

Time 


349.24 


377.55 


379.54 


380.22 


387.73 


396.03 


396.71 


All .2^ 


434.15 


434.83 


437.04 


437.72 


443.40 


456.58 


457.24 


457.92 


459.69 


465.66 


466.34 


483.97 


499.04 


503.40 


504.08 


509.51 


516.97 


517.65 


521.73 


525.92 


551.90 


555.80 


556.38 


557.06 


563.60 


572.05 


590.64 


591.20 


591.88 


Component 

Failure 


C2 


C2 


Cl 


C3 


C2 


Cl 


C3 


C2 


Cl 


C3 


Cl 


C3 


C2 


C2 


Cl 


C3 


C2 


Cl 


C3 


C2 


C2 


Cl 


C3 


C2 


Cl 


C3 


C2 


C2 


C2 


C2 


Cl 


C3 


C2 


C2 


C2 


Cl 


C3 
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OT  Data  Runs  #1  and  #2 


OT  Data  Run  #1 

OT  Data  Run  #2 

Time 

Component 

Failure 

Time 

Component 

Failure 

11.88 

C2 

27.30 

C2 

15.85 

Cl 

31.95 

C2 

19.81 

C3 

36.40 

Cl 

24.52 

C2 

41.05 

C3 

30.77 

C2 

42.60 

Cl 

32.69 

Cl 

45.94 

C2 

36.65 

C3 

47.25 

C3 

40.69 

C2 

47.72 

C2 

41.03 

Cl 

58.60 

C2 

44.99 

C3 

61.25 

Cl 

54.26 

Cl 

63.62 

Cl 

58.22 

C3 

65.90 

C3 

87.86 

C2 

68.27 

C3 

76.87 

C2 

78.14 

Cl 

82.79 

C3 

87.02 

C2 

88.69 

C2 

92.19 

C2 

102.49 

Cl 

105.31 

C2 

107.14 

C3 

116.03 

Cl 

118.26 

Cl 

120.68 

C3 

122.91 

C3 

122.92 

Cl 

82 


J!^>pendix  D.  Calculation  For  Series  System  Reliability. 

Consider  a  system  composed  of  k  =  3  independent 

components  in  series  configuration.  The  time  to  failure  of 
each  component  is  assumed  to  be  exponentially  distributed. 
Let  Xi  =  0.05,  ^2  =  0.067,  and  A,3  =  0.05  be  the  respective 
constant  failure  rates  (failure  per  hour)  of  each  component. 
The  failure  rate  for  the  system  is 

3 

X-  0.05  +  0.067  +  0.05  =  0.167  failure  per  hour. 

i=l 

Thus  the  mean  failure  for  the  system  is 

■  (t)'  (0:167)  = « 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  and  3  are 

Ri  =  Rs  =  exp(-0.05)  =  95.12% 

R2  =  exp(-0.067)  =  93.52%. 

The  reliability  of  the  system  can  be  calculated  as 

Rs  =  R1XR2XR3  =  95.12%x93.52%x95.12%  =  84.62% 

or  Rs  =  exp{-Xt)  =  exp  (-0.167)  =  84.62%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  C) : 
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Table  D1 :  System  Level  DT  Data 


Table  D2;  System  Level  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

87.86 

13 

0.14796 

2 

122.92 

27 

0.21966 

Table  D3;  Component  1  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

499.86 

24 

0.04801 

2 

449.53 

20 

0.04449 

3 

384.76 

13 

0.03379 

4 

543.47 

29 

0.05336 

5 

591.88 

21 

0.03548 

Table  D4 ;  Component  1  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

87.86 

4 

0.04553 

2 

122.92 

9 

0.07322 

Table  D5:  Component  2  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

499.86 

28 

0.05602 

2 

449.53 

37 

0.08231 

3 

384.76 

20 

0.05198 

4 

543.47 

39 

0.07176 

5 

591.88 

33 

0.05576 

Run  i 


Ti  Hours 

Xi  Failures 

87.86 

5 

122.92 

10 

0.05691 


0.08135 


Table  D7: 


Run  i 

Ti  Hours 

Xi  Failures 

1 

499.86 

28 

2 

449.53 

21 

3 

384.76 

13 

4 

543.47 

29 

5 

591.88 

21 

0.05602 


0.04672 


0.03379 


0.05336 


0.03548 
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Table  D8 :  Component  3  01 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

h  =  Xi/Ti 

1 

87.86 

4 

0.04553 

2 

122.92 

8 

0.06508 

The  Bayesian  Method  1. 

Bayesian  Prior  Analysis .  Assuming  that  the  prior 
belief  about  the  failure  rates  Xi  (i  =  1,  2,  3)  is 
exponentially  distributed  with  parameter  6.  Then  the  prior 

3 

belief  about  X  =  is  defined  by  a  gamma  distribution 

with  pdf  given  in  Equation  27.  The  mean  of  the  system  is 
E{X)  =  k/0.  With  E(X)  =  6  from  above,  solving  for  0  to 


obtain  0 


(bW)) 


0.5. 


Bayesian  Posterior  Analysis. 

At  the  system  level  OT  data,  there  are  40  failures  (f? 
=  40)  observed  during  the  total  test  time  T2  =  210.78  hours, 
k  =  3  and  X-  0.167  failure  per  hour.  Thus 

2  =  -  (21o!?8%^.5)  =  0-20352. 


At  the  system  level  OT  data,  there  are  40  failures  {£2 
=  40)  observed  during  the  total  test  time  Ta  =  210.78  hours, 

k  =  3  and  X  =  0.167  failure  per  hour.  Thus 

^  \  40-t-3 


^(t) 


1  + 


210.78+0.57  J 


=  0.8162. 
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Computations  Using  The  Exponential  Model. 

The  Classical  Exponential  Method.  Assuming  the  time  to 
failure  of  each  component  has  the  exponential  distribution, 
point  estimates  for  the  failure  rate  K  and  80%  confidence 


intervals  about  X,  will  be 

computed 

• 

System  Level  DT: 

_  f 

376  ^ 

0.15226 

l  r  J 

2469.57 

System  Level  OT: 

(I2.) 

_  f 

40  ^ 

0.18977 

1  Tz) 

210.787 

Component  1  DT: 

_  f 

107  ^ 

0.04333 

V  Tz) 

2469. 5>^ 

Component  1  OT: 

_  f 

. 13_.._.^ 

0.06168 

KTa) 

210.787 

Component  2  DT: 

1 

_  f 

157  ^ 

0.06358 

{TsJ 

2469.57 

Component  2  OT: 

(Is.) 

_  f 

15  >1 

0.07116 

V  Te) 

-  1 

210.787 

Component  3  DT: 

112  ^ 

0.04535 

1  r,; 

-  1 

2469.57 

Component  3  OT: 

2 

(I£) 

_  f 

12  ^ 

0.05693 

VTa) 

-  1 

210.787 

where  fi  =  376  is  the  number  of  failures  and  Ti  =  2469.5  is 
the  total  test  time  from  Table  D1  and  similarly  fa  =  40  and 
Ta  =  210.78  from  Table  D2,  fs  =  107  and  Ts  =  2469.5  from 
Table  D3,  ft  =  13  and  T4  =  210.78  from  Table  D4,  fs  =  157  and 
Ts  =  2469.5  from  Table  D5,  f&  =  15  and  Te  =  210.78  from  Table 
D6,  f?  =  112  and  T7  =  2469.5  from  Table  D7,  fe  =  12  and  Ts  = 
210.78  from  Table  D8. 
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The  80%  two-sided  confidence  intervals  (TCI)  are 


(  Zo.2/2(2x376) 

^1-0.2/2(2x376) 

V  2X2469.5 

r  2x2469.5 

System  Level 

OT:  (0.15248 

(  /o. 2/2(2x40) 

^1-0.2/2(2x40) 

2x210.78  f 

2x210.78  J 

Comoonent  1 

DT:  (0.03806, 

(  ;fg.2/2(2xl07) 

^1-0.2/2(2x107) 

\  2x2469.5 

'  2x2469.5 

Component  1 

OT:  (0.04102, 

f  /o. 2/2(2x13) 

^1-0.2/2(2x13) 

^  2x210.78  ' 

2X210.78  J 

Component  2 

DT:  (0.05717,1 

(  ;f§.2/2(2xl57) 

^1-0.2/2(2x157) 

2x2469.5 

'  2x2469.5 

Component  2 

OT:  (0.04886,1 

/  Jfo. 2/2(2x15) 

1-0. 2/2  (2x1 5)  'N 

^  2x210.78  ' 

2x210.78  J 

Comoonent  3 

DT:  (0.03996,1 

/  ^0.2/2(2x112) 

/i-o.2/2(2x112) 

V  2x2469.5 

'  2x2469.5 

Comoonent  3 

OT:  (0.03715,1 

f  702.75  802.11^ 

“  V  '  4939  J 


(  64.28  96.58  ^ 

“  421.56  r  421.56  ) 


187.95  240.9 

4939  '  4939 


) 


1  _  r  17.29  35.56  ^ 

J  ~  \  421.56  r  421.56  J 

7016) 

(  282.34  346.51^ 

~  V  ^939  '  4939  J 


(  Xo.2/2(^^l^)  ^1-0.2/2(2^12) ^  15.66  33.2 

2x210.78  '  2x210.78  J  =  421.56  '  421.56 ) 

The  reliability  estimators  of  R(t)  for  t  =  1  hour  are 
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System  Level  DT:  R{t)  =  exp{-0. 15226X1)  =  85.88% 
System  Level  OT:  ■R(t)  =  exp  (-0, 18977X1)  =  82.72% 
Component  1  DT :  Ri{t)  =  exp (-0 . 04333X1 )  =  95.76% 

Component  1  OT :  Ri(t)  =  exp (-0 . 06168X1 )  =  94.02% 

Component  2  DT:  Rz  (t)  =  exp (-0. 06358X1)  =  93.84% 

Component  2  OT:  Rzit)  -  exp(-0. 07116X1)  =  93.13% 

Component  3  DT;  Rsit)  =  exp(-0. 04535X1)  =  95.57% 

Component  3  OT;  Rsit)  =  exp(-0. 05693X1)  =  94.47% 

Component  Level  Aggregated  to  System  DT:  R(t)  =  85.88% 
R{t)  =  Ri*R2*R3  =  0.9576*0.9384*0.9557  =  0.8588 
Component  Level  Aggregated  to  System  OT;  R(t)  =  82.72% 

R{t)  =  Ri*R2*R3  =  0.9402*0.9313*0.9447  =  0.8272. 

The  80%  two-sided  confidence  intervals  (TCI)  are 
System  Level  DT:  (85. 01%, 86. 74%) 

exp(-0. 16240X1)  <  R{1)  <  exp  (-0. 14229X1) 

System  Level  OT:  (79. 53%, 85. 86%) 
exp(-0. 22910X1)  ^  R{1)  ^  exp (-0 . 15248X1 ) 

Component  1  DT;  (95. 24%, 96. 27%) 
exp{-0. 04878X1)  <  R(l)  <  exp (-0. 03806X1) 

Component  1  OT:  (91. 91%, 95. 98%) 
exp(-0. 08436X1)  <  R(l)  ^  exp  (-0. 04102X1) 


exp{-0. 07016X1)  <  R(l)  <  exp  (-0. 05717X1) 

Component  2  OT;  (90. 89%, 95. 23%) 

exp(-0. 09549X1)  ^  ;R(1)  <  exp  (-0. 04886X1) 

Component  3  DT:  (95. 03%, 96. 08%) 

exp(-0. 05093X1)  <  ^(1)  <  exp (-0 . 03996X1 ) 

Component  3  OT:  (92 .43%, 96.35%) 

exp(-0. 07875X1)  <  ^(1)  <  exp(-0. 03715X1) 

Component  Level  Aggregated  to  System  DT: 

(84. 37%, 87. 35%) 

Lower  Bound  =  0.9524*0.9322*0.9503  =  0.8437 
Upper  Bound  =  0.9627*0.9444*0.9608  =  0.8735 
Component  Level  Aggregated  to  System  OT: 

(77. 21%, 88. 07%) 

Lower  Bound  =  0.9191*0.9089*0.9243  =  0.7721 
Upper  Bound  =  0.9598*0.9523*0.9635  =  0.8807. 

The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior. 

For  the  system  level  data,  from  Table  Dl: 

£(>1)  =  T  =  0.15166  and  V(X)  =  sx"  =  0.00073. 

These  yielded  the  parameter  estimates 

“  '  (  o'.WotI')  =  p  =  -  207.76. 

For  Component  1  data,  from  Table  D3: 
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£(>.)  =  i  =  0.04303  and  V{X)  =  Sx"  =  0.00007. 

These  yielded  the  parameter  estimates 

=  623.58, 


_  f  0.04303^^  _ 
^  -  V  0.00007  J  - 


26.83  and  P  =  (  Q-Q43Q.3^ 


00007; 

For  Component  2  data,  from  Table  D5: 

E{X)  =  T  =  0.06356  and  V(X)  =  Sx^  =  0.00017. 

These  yielded  the  parameter  estimates 

„  _  f  0.06356^^  _  --  __  ^  _  f  0.06356^ 

“  “  I  0.00017  J  -  23.76  and  p  -  0.OOOI7J 

For  Component  3  data,  from  Table  D7: 

EiX)  =  T  =  0.04507  and  V(X)  =  Sx^  =  0.0001. 

These  yielded  the  parameter  estimates 

f  0.04507^^  _  ,  ^  _  f  0.04507^  _ 

“  “  V  0.0001  J  ~  20.31  and  P  -  0.0001  )  ~ 

Bayesian  Exponential  Posterior  Analysis . 

51 


=  373.88. 


450.7. 


Component  1: 
Component  2: 
Component  3; 


210.78+207. 

76. 

13+26.83 

210.78+623. 

58, 

15+23.76 

210.78+373. 

88. 

12+20.31 

)- 


0.17086, 


^  -  o.OAses. 


210.78- 

The  80%  two-sided  probability  intervals  (TBPI)  are 


System  Level: 


f  ;(ro  2/2(2x40+2x31.51)  y?_o  2/2(2x40+2x31 . 51)  ^ 
\  2(210.78+207.76)  '  2(210.78+207.76)  J 


f  121.8  165. 06  A  _ 

“  837.08  '  837. 08y  ~  <0.14551,0.19718). 
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f  ;fo. 2/2(2x13+2x26. 83)  2/2(2x13+2x26 . 83) 

Component  1:  2(210.78+623.58)  '  2(210.78+623.58) 


(0.03798,0.05722) 


'  (  1668??2  >  166i^?2')  *  (0.03798,0.05722). 

f  2/2(2x15+2x23. 76)  2/2(2x15+2x23 . 76) 

Component  2;  2(210.78+373.88)  '  2(210.78+373.88) 


f  61.59 
1169.32  f 


Component  3 : 


93.27 

1169.32 


(0.05267,0.07976) . 


0.2/2(2x12+2x20 . 31)  2/2(2x12+2x20 . 31)  ^ 

2(210.78+450.7)  '  2(210.78+450.7)  J 


132^2'.  96V  1322:  96  )  =  (0.03779,0.05961). 


The  Bayesian  estimators  of  R(t)  for  t  =  1  hour  are 


System  Level:  R(t}  =  84.31% 


R{t) 


r  210^/7£^  ,  T 

I  201. leJ 

f  210. ..78.^  ( _ 1 _ ^ 

V  201.16)  207.76; 


31.51+40 


Component  1 :  R{t)  =  95.34% 


R{t)  = 


f  210.78^ 
V  623.58; 

210^08.^  j.f _ 

I  623.58;  62: 


1 

623.58 


)+l 


26.83+13 


Component  2:  R[t)  =  93.59% 


R{t)  = 


f  210.78^  .  T 

V  373.88; 

(  210.78^  J  _l_\ 
\  373.88;  373.88; 


23.76+15 


Component  3;  R(t)  =  95.24% 


R{t) 


I  450.7  ; 

f  2.1.0. 78.S  +r_x_')  +■ 

V  450.7  ;  450.1  J 


20.31+12 
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Component  Level  Aggregated  to  System:  ■R(t)  =  84.98% 

R{t)  =  0.9534*0.9359*0.9524  =  0.8498 

The  80%  two-sided  probability  intervals  (TBPI)  are 

System  Level;  (82.10%  <  ^(1)  <  86.46%) 

exp(-0. 19718X1)  <  ^(1)  <  exp  (-0. 14551X1)  . 

Component  1;  (94.44%  <  ^(1)  <  96.27%) 
exp(-0. 05722X1)  <  :R(1)  <  exp (-0. 03798X1)  . 

Component  2:  (92.33%  <  ^Rd)  ^  94.87%) 
exp(-0. 07976X1)  <  R{1)  ^  exp  (-0 . 05267X1)  . 

Component  3;  (94.21%  <  R(l)  <  96.29%) 

exp(-0. 05961X1)  <  R(l)  <  exp  (-0 . 03779X1 )  . 

Component  Level  Aggregated  to  System; 

(82.15%  <  R{1)  <  87.94%) 

0.9444*0.9233*0.9421  ^  R{1)  ^  0.9627*0.9487*0.9629. 
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J^pendix  E.  Simulated  Data  For  Parallel  System  With 
Independent  And  Identical  Con^onents. 

DT  Data  Run  #1 


Time 

Component  Failure 

Time 

Component  Failure 

12.04 

C3 

282.71 

C2 

12.06 

Cl 

289.01 

C2 

18.71 

C3 

291.78 

C3 

32.47 

C3 

296.17 

Cl 

33.23 

C2 

304.14 

C2 

37.16 

Cl 

306.46 

C2 

37.35 

Cl 

311.16 

C2 

41.36 

C3 

315.75 

Cl 

42.37 

Cl 

316.15 

C2 

46.02 

Cl 

316.49 

C3 

62.38 

C2 

326.02 

C2 

68.33 

C3 

331.25 

C2 

77.62 

C3 

331.31 

C3 

79.76 

C3 

334.20 

C2 

107.82 

C2 

335.24 

Cl 

108.15 

C2 

342.72 

Cl 

108.76 

C2 

351.61 

C3 

116.38 

C2 

352.56 

Cl 

120.51 

C2 

355.99 

C2 

123.73 

C2 

356.56 

C2 

159.04 

C3 

370.17 

Cl 

163.75 

Cl 

371.68 

Cl 

165.04 

C3 

385.22 

Cl 

171.15 

C2 

387.04 

Cl 

175.61 

C3 

396.62 

C2 

192.46 

C3 

404.57 

C3 

192.59 

Cl 

414.25 

C2 

199.20 

C2 

424,79 

C2 

201.88 

C2 

427.24 

C2 

205.62 

Cl 

431.17 

Cl 

206.64 

C2 

442.25 

Cl 

214.41 

Cl 

443.08 

C2 

217.39 

Cl 

444.53 

Cl 

220.35 

C2 

444.94 

C3 

226.62 

C2 

445.11 

Cl 

230.28 

C2 

451.56 

C3 

230.81 

C2 

458.01 

C3 

246.24 

Cl 

495.84 

C2 

251.66 

C3 

497.94 

C3 

255.25 

C2 

497.94 

C2 
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DT  Data  Run  #2 


Time 

Component 

Failure 

4.20 

C3 

4.58 

Cl 

6.07 

C2 

17.53 

C2 

18.43 

Cl 

34.77 

C3 

46.48 

C3 

57.70 

C2 

62.53 

C3 

68.69 

Cl 

75.71 

C2 

90.38 

C3 

90.50 

Cl 

91.35 

C2 

93.99 

C3 

94.21 

Cl 

105.29 

C3 

107.87 

Cl 

110.82 

C3 

114.35 

C3 

135.40 

C2 

142.04 

Cl 

142.34 

Cl 

153.25 

Cl 

159.69 

C2 

161.91 

Cl 

168.85 

C2 

170.29 

C2 

177.01 

C3 

208.26 

Cl 

216.50 

C3 

218.39 

Cl 

224.21 

C2 

230.03 

C3 

230.64 

Cl 

248.50 

C2 

270.19 

C3 

274.10 

C2 

Time 


21b. bO 


281.06 


290.49 


295.71 


304.23 


306.66 


313.05 


336.41 


336.81 


337.86 


344.63 


351.90 


361.93 


363.04 


365.05 


368.64 


389.02 


391.29 


398.75 


403.27 


415.59 


425.97 


427.24 


428.28 


429.95 


431.16 


441.17 


449.38 


451.03 


455.09 


458.15 


459.00 


459.80 


463.57 


464.11 


475.82 


487.71 


Component 

Failure 


C2 


Cl 


C3 


C2 


Cl 


Cl 


C2 


C2 


Cl 


C3 


Cl 


C2 


Cl 


Cl 


C2 


C2 


C3 


Cl 


C2 


C2 


C2 


C3 


Cl 


C2 


Cl 


Cl 


C2 


C2 


C3 


-C2 


C2 


C2 


Cl 


C2 


Cl 


C3 


C2 


DT  Data  Run  #3 


Time 

Component 

Failure 

Time 

Component 

Failure 

9.08 

C2  ^ 

255.75 

C3 

12.18 

C3 

259.70 

Cl 

26.30 

Cl 

261.80 

C3 

40.09 

C3 

266.50 

Cl 

40.41 

C3 

288.06 

C2 

61.50 

C2 

288.64 

C3 

61.75 

C3 

292.05 

C2 

64.21 

C2 

294.44 

C3 

69.80 

C3 

295.71 

Cl 

71.06 

C3 

334.44 

C2 

71.24 

Cl 

334.85 

C2 

74.24 

C2 

336.79 

C2 

105.08 

C2 

356.86 

Cl 

110.15 

C2 

357.14 

Cl 

^  110.28 

C3 

366.86 

C3 

122.39 

Cl 

367.70 

C3 

123.14 

C2 

371.11 

Cl 

134.03 

C3 

376.81 

C3 

138.38 

C3 

378.97 

C2 

142.06 

C3 

380.16 

Cl 

168.67 

Cl 

381.50 

Cl 

169.65 

C3 

385.48 

Cl 

172.87 

Cl 

389.79 

Cl 

179.59 

C2 

392.03 

Cl 

189.06 

C3 

401.88 

C2 

189.37 

Cl 

436.54 

C2 

192.59 

C2 

448.39 

C3 

217.46 

C2 

448.48 

Cl 

220.53 

C2 

460.88 

Cl 

221.95 

C2 

494.32 

C3 

233.27 

Cl 

494.72 

Cl 

237.80 

C2 

495.58 

C3 

251.11 

C3 

497.53 

Cl 
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DT  Data  Run  #4 


Time 

Component 

Failure 

Time 

Component 

Failure 

7.12 

C3 

250.70 

C3 

22.95 

Cl 

256.25 

Cl 

36.55 

Cl 

262.69 

Cl 

37.98 

C3 

278.48 

C2 

55.79 

C3 

280.28 

C3 

59.30 

C2 

281.44 

C2 

70.99 

C3 

286.02 

C3 

77.46 

Cl 

291.79 

C3 

96.21 

C2 

299.82 

C2 

100.76 

C3 

300.73 

C3 

104.27 

C2 

302.40 

C3 

105.69 

Cl 

321.65 

C3 

112.94 

C2 

324.77 

C3 

120.84 

C3 

327.33 

C3 

124.10 

C2 

338.50 

Cl 

125.88 

Cl 

349.24 

C2 

130.52 

C2 

357.36 

Cl 

130.76 

C3 

360.37 

C3 

144.16 

Cl 

365.21 

C2 

150.61 

C3 

366.03 

C2 

154.18 

C3 

366.04 

Cl 

156.54 

Cl 

380.64 

C2 

163.03 

Cl 

389.64 

C2 

166.15 

Cl 

418.84 

C3 

167.21 

Cl 

427.39 

Cl 

179.93 

C2 

429.46 

C3 

195.27 

C3 

437.45 

C3 

201.07 

C3 

449.55 

Cl 

203.93 

C2 

455.18 

^  C3 

204.56 

C2 

455.99 

Cl 

205.88 

C3 

475.72 

C3 

215.39 

C3 

479.32 

C3 

228.63 

C2 

491.46 

Cl 

243.40 

Cl 

494.00 

C3 

245.96 

C3 

496.95 

C2 
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DT  Data  Run  #5 


Time 

Component 

Failure 

10.07 

Cl 

38.94 

Cl 

51.73 

C3 

57.02 

C3 

59.39 

Cl 

73.67 

C2 

77.17 

C2 

79.56 

C2 

81.77 

Cl 

83.79  ! 

C2 

102.76 

C2 

111.46 

Cl 

113.76 

C2 

117.40 

C3 

120.92 

C2 

140.39 

Cl 

140.93 

C3 

143.52 

C3 

143.67 

C3 

150.00 

Cl 

153.58 

Cl 

157.12 

Cl 

160.97 

C2 . 

164.65 

C3 

185.65 

C3 

186.87 

C3  ! 

188.08 

Cl 

190.84 

C2 

192.22 

C3 

197.47 

C3 

199.83 

C2 

200.71 

Cl 

202.67 

C2 

203.83 

Cl 

215.12 

C3 

222.97 

C3 

243.53 

C3 

266.85 

C3 

275.75 

C2 

276.66 

C3 

302.19 

C2 

Time 


311.09 


311.42 


312.23 


317.40 


321.72 


325.49 


338.23 


343.49 


346.16 


351.11 


359.99 


364.48 


367.25 


373.64 


374.40 


376.93 


384.09 


394.25 


401.81 


402.77 


417.35 


418.39 


422.69 


425.16 


426.57 


437.34 


438.51 


442.89 


443.85 


444.85 


447.39 


449.22 


455.79 


457.86 


460.98 


462.83 


479.77 


486.68 


490.60 


496.86 


Component 

Failure 


Cl 


Cl 


C3 


Cl 


C2 


C2 


Cl 


C2 


C3 


C3 


C2 


C3 


C2 


Cl 


C3 


C2 


C2 


Cl 


C2 


C3 


C3 


C3 


C3 


C2 


C2 


Cl 


Cl 


C2 


C3 


Cl 


C3 


C2 


C3 


Cl 


C3 


Cl 


C2 


C2 


C2 


C3 
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OT  Data  Runs  #1  and  #2 


OT  Data  Run  #1 

OT  Data  Run  #2 

Time 

Component 

Failure 

Time 

Component 

Failure 

4.06 

C3 

11.34 

C3 

8.27 

C2 

18.10 

C3 

15.64 

Cl 

29.26 

C3 

16.26 

C2 

35.38 

C2 

17.01 

C2 

36.41 

C3 

29.60 

Cl 

51.17 

C3 

40.19 

C2 

55.38 

«  Cl 

43.63 

C3 

75.49 

Cl 

44.47 

C3 

79.78 

C2 

57.77 

C2 

84.86 

Cl 

62.95 

C2 

88.63 

C2 

66.03 

Cl 

96.55 

C3 

73.51 

C3 

103.25 

C2 

73.69 

C3 

111.41 

C2 

80.82 

Cl 

115.54 

C2 

85.21 

Cl 

92.91 

C3 

100.67 

Cl 

101.59 

C2 

110.15 

C2 

123.20 

Cl 

124.42 

C3 
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Appendix  F.  Calculation  For  Parallel  System  Reliability 
With  Independent  and  Identical  Components. 

Consider  a  system  composed  of  k  =  3  independent  and 
identical  components  in  parallel  configuration.  The  time  to 
failure  of  each  component  is  assumed  to  be  exponentially 
distributed.  Let  X  -  0.05  failure  per  hour  be  the 

respective  constant  failure  rates  of  each  component. 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  and  3  are 
Rc  =  Ri  =  Ra  =  Rs  =  exp  (-0.05)  =  95.12%. 

The  reliability  of  the  system  can  be  calculated  as 
Rs  =  1  -  (1  -  Rc)'  =  1  -  (1  -  95.12%)^  =  99.9884%. 

Using  a  computer  program  called  BUIPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  E) : 


Table  FI;  Aggregated  Components  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

1,493.82 

79 

0.05289 

2 

1,463.13 

75 

0.05126 

3 

1,492.59 

66 

0.04422 

4 

1,490.85 

70 

0.04695 

5 

1,490.58 

81 

0.05434 

Table  F2;  Aggregated  Components  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

373.26 

22 

0.05894 

2 

346.62 

15 

0.04328 
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Table  F3:  Component  1  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

497.94 

24 

0.04820 

2 

487.71 

26 

0.05331 

3 

497.53 

22 

0.04422 

4 

496.95 

20 

0.04025 

5 

496.86 

23 

0.04629 

Table  F4 :  Component  1  01 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

Ti  =  Xi/Ti 

1 

124.42 

7 

0.05626 

2 

115.54 

3 

0.02597 

Table  F5:  Component  2  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

497.94 

35 

0.07029 

2 

487.71 

30 

0.06151 

3 

497.53 

21 

0.04221 

4 

496.95 

19 

0.03823 

5 

496.86 

28 

0.05635 

Table  F6:  Component  2  01 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

Ti  =  Xi/Ti 

1 

124.42 

8 

0.06430 

2 

115.54 

6 

0.05193 

Table  F7 :  Component  3  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

497.94 

20 

0.04017 

2 

487.71 

19 

0.03896 

3 

497.53 

23 

0.04623 

4 

496.95 

31 

0.06238 

5 

496.86 

30 

0.06038 

Table  F8:  Component  3  OT 

Data 

Run  i 

Ti  Hours 

Ti  =  Xi/Ti 

1 

124.42 

7 

0.05626 

2 

115.54 

6 

0.05193 

100 


The  Classical  Exponential  Method.  Assuming  the  time  to 
failure  of  each  component  has  the  exponential  distribution, 
point  estimates  for  the  failure  rate  X  and  80%  confidence 
interval  about  X  will  be  computed. 

Component  Aggregated  DT:  X  =  (  7 4 30^ '97')  “  0.04993 
Component  Aggregated  OT:  X  =  (  71 9^88 )  ~  0.05140 
Component  1  DT:  X  =  ~  (  2476^99)  ~  0.04643 
Component  1  OT:  X  =  ~  (  239~T6)  ~  0.04167 
Component  2  DT:  X  -  =  (  2476^99)  0.05369 
Component  2  OT:  X  =  =  (  239^96)  ^  0.05834 
Component  3  DT:  1  =  (  2476^^)  ""  0.04966 

Component  3  OT:  2  =  (^)  =  ("239^^)  ^  0.05418 

where  fi  =  371  is  the  number  of  failures  and  Ti  =  7430.97  is 
the  total  test  time  from  Table  FI  and  similarly  fa  =  37  and 
Ta  =  719.88  from  Table  F2,  fa  =  115  and  T3  =  2476.99  from 
Table  F3,  f4  =  10  and  T4  =  239.96  from  Table  F4,  fs  =  133  and 
Ts  =  2476.99  from  Table  F5,  fe  =  14  and  Ts  =  239.96  from 
Table  F6,  fa  =  123  and  T7  =  2476.99  from  Table  F7,  fg  =  13 
and  Ts  =  239.96  from  Table  F8. 
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The  80%  two-sided  confidence  intervals  (TCI)  are 


Component  Aggregated  DT:  (0.04664,0.05328) 


( 


/0.2/2(2x371)  /to, 2/2(2x37 
2x7430.97  '  2x7430.97  J 


=  ( 


693.08  791 

14861.94  f  14861 


.78  ^ 
1.94  ^ 


Component  Aggregated  OT;  (0.04091,0.06248) 

f  .yp. 2/2(2x37)  JCi-o, 2/2(2x37)^  /  58 . 9  89.96  ^ 

2x719.88  '  2x719.88  y  “  1^  1439.76  '  1439. 76^ 


Component  1  DT ;  (0.04097,0.05206) 


( 


Xl2/2i2xll5)  /!-o.2/2(2X115)>| 
2x2476.99  '  2x2476.99  ^ 


=  ( 


202.98 


257 


4953.98 '  4953 


88  ^ 
.98  j 


Component  1  OT:  (0.02593,0.05920) 


( 


/o. 2/2(2x10)  Jfi-o. 2/2(2x10)  ^ 
2x239.96  '  2X239.96  J 


(  12.44  28.41 

~  479.92  '  479.92  J 


Component  2  DT :  (0.04782,0.05974) 


(4 


2/2(2x133)  ;fLo.2/2(2Xl33)>> 

2x2476.99  '  2x2476.99  y 


236.90 


295 


4953.98  '  4953 


.95  ^ 
1.98  y 


Component  2  OT:  (0.03946,0.07901) 


( 


/o.2/2(2x14)  Afl-0. 2/2(2x14)  ^ 
2x239.96  f  2x239.96  j 


f  18.94  37.92  ^ 

“  479.92  '  479.92  y 


Component  3  DT:  (0.04401,0.05548) 


( 


^0.2/2(2x123)  jfi-o. 2/2(2x12 3) ^ 


2x2476.99  '  2x2476 


X123)  A  / 

•99  j  =  V 


274 


218.04 
4953.98  '  4953 


.82  ^ 

1.98  y 


Component  3  OT:  (0.03603,0.07410) 

f  ;fo. 2/2(2x13)  x!-o. 2/2(2x13) "X  /  17.29189  35.56317  ^ 
\  2x239.96  '  2x239.96  J  ~  \  479.92  '  479.92  y 
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The  reliability  estimators  of  R(t)  for  t  =  1  hour  are 
Component  Aggregated  DT; 


Rc it)  =  exp(-0. 04993X1)  =  95.13% 

Component  Aggregated  OT: 

^(t)  =  exp(-0. 05140X1)  =  94.99% 

Component  1  DT;  Ri {t)  =  exp (-0 . 04643X1 )  = 
Component  1  OT;  i?i(t)  =  exp(-0. 04167X1)  = 
Component  2  DT:  Rz  (t)  =  exp(-0. 05369X1)  = 

Component  2  OT:  Rz{t)  =  exp (-0 . 05834X1 )  = 

Component  3  DT:  Rsit)  =  exp(-0. 04966X1)  = 

Component  3  OT:  Rz  {t)  =  exp (-0. 05418X1 )  = 

Aggregated  Component  to  System  DT:  R{t)  = 
R{t)  =  1  -  (1  -  Rc)^  =  1  -  (1  -  95.13%)^  = 
Aggregated  Component  to  System  OT:  R{t)  = 

R(t)  =  1  -  (1  -  Rc)^  =  1  -  (1  -  94.99%)3  = 
Component  Level  Aggregated  to  System  DT: 
R{t)  =  1  -  (1  -  Ri)  (1  -  Ra)  {1  -  R3) 

=  1  -  (1  -  0.9546) (1  -  0.9477) (1  -  0.9516) 
Component  Level  Aggregated  to  System  OT; 
R{t)  =  1  -  (1  -  Ri)  (1  -  R2)  (1  -  R3) 

=  1  -  (1  -  0.9592) (1  -  0.9433) (1  -  0.9473) 


95.46% 

95.92% 

9A.n\ 

94.33% 

95.16% 

94.73% 

99.9885% 

0.999885 

99.9874% 

0.999874 

=  99.9885% 

=  99.9878%. 
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The  80%  two-sided  confidence  intervals  (TCI)  are 


Aggregated  Component  DT;  (94, 81%, 95. 44%) 

exp(-0. 05328X1)  <  :R(1)  <  exp(-0. 04664X1) 

Aggregated  Component  OT:  (93 . 94%, 95 . 99%) 

exp(-0. 06248X1)  <  ^Rd)  <  exp (-0 . 04091X1 ) 

Component  1  DT:  (94. 93%, 95. 99%) 

exp(-0. 05206X1)  <  R(l)  <  exp (-0 . 04097X1 ) 

Component  1  OT:  (94. 25%, 97. 44%) 

exp(-0. 05920X1)  <  R{1)  <  exp  (-0 . 02593X1 ) 

Component  2  DT:  (94. 20%, 95. 33%) 

exp(-0. 05974X1)  <  R(l)  <  exp (-0 . 04782X1 ) 

Component  2  OT:  (92 . 40% , 96 . 13%) 

exp(-0. 07901X1)  <  R(l)  <  exp  (-0.03946X1) 

Component  3  DT;  (94.60%, 95. 69%) 

exp(-0. 05548X1)  <  R(l)  <  exp (-0 . 04401X1) 

Component  3  OT:  (92. 86%, 96. 46%) 

exp(-0. 07410X1)  <  R(l)  <  exp (-0 . 03603X1 ) 

Aggregated  Component  to  System  DT: 


(99. 9860%, 99. 9905%) 

Lower  Bound  =  1  -  (1  -  0.9481)^  =  0.999860 
Upper  Bound  =  1  -  (1  -  0.9544)^  =  0.999905 
Aggregated  Component  to  System  OT: 

(99. 9778%, 99. 9936%) 
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Lower  Bound  =  1 


(1  -  0.9394)3 


0.999778 


Upper  Bound  =  1  -  (1  -  0.9599)3  _  0.999936 
Component  Level  Aggregated  to  System  DT; 

(99. 9841%, 99. 9919%) 

Lower  Bound  =  1  -  (1  -  LBi)  (1  -  LBz)  (1  -  LB3) 

=  1  -  (1  -  0.9493) (1  -  0.9420) (1  -  0.9460)  =  0.999841 
Upper  Bound  =  1  -  (1  -  UBi)  (1  -  UBz)  (1  -  UB3) 

=  1  -  (1  -  0.9599) (1  -  0.9533) (1  -  0.9569)  =  0.999919 
Component  Level  Aggregated  to  System  OT: 

(99. 9688%, 99. 9965%) 

Lower  Bound  =  1  -  (1  -  LBi)  (1  -  LB2)  (1  -  LB3) 

=  1  -  (1  -  0.9425)  (1  -  0.9240)  (1  -  0.9286)  =  0.999688 

Upper  Bound  =  1  -  (1  -  UBi)  (1  -  UB2)  (1  -  UB3) 

=  1  -  (1  -  0.9744) (1  -  0.9613) (1  -  0.9646)  =  0.999965. 

The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior. 

For  the  component  aggregated  data,  from  Table  FI: 

E{K)  =  T  =  0.04993  and  V(A,)  =  =  0.00002. 


These  yielded  the  parameter  estimates 


f  0.049932^  ca  ,  n 

“  =  I  0.00002  )  =  139-59  and  p  = 


f  0.04993^ 
V  0.000027 


=  2795.66. 


For  Component  1  data,  from  Table  F3: 

E{X)  =  T  =  0.04645  and  V{X)  =  Sx^  =  0.00002. 


These  yielded  the  parameter  estimates 
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a  =  ( 


0.04645^ 

0.00002 


)  -  92.19  and  p  -  (q!  00002 )  1984.67. 


For  Component  2  data,  from  Table  F5: 


E{X}  =  A  =  0.05372  and  V{X)  =  =  0.00018. 


These  yielded  the  parameter  estimates 


f  0.05372^^ 

V  0.00018  J 


16.15  and  p  = 


f  0.05372^ 

V  o.oooiey 


300.71. 


For  Component  3  data,  from  Table  F7 : 


E{X)  =  A  =  0.04962  and  V{X)  =  =  0.00012. 


These  yielded  the  parameter  estimates 


f  0.04962^^ 

I  0.00012  J 


19.97  and 


0.04962^ 

0.00012A 


402.48. 


Bayesian  Exponential  Posterior  Analysis, 

Aporecrated  Component t  A“  V.  719  884“2795  SsJ  *  0.05023. 


Component  1 : 

Component  2 : 

Component  3: 


A  =  1 

(  10+92.19  ^ 

s  0.04594. 

1 239.96+1984.67; 

A  =  1 

r  14+16.15  ^ 

=  0.05576. 

k  239.96  +  300.71; 

A  =  1 

f  13+19.97  ^ 

=  0.05132. 

^  239.96  +  402.48; 

The  80%  two-sided  probability  intervals  (TBPI)  are 
Aggregated  Component: 

f  ;fi2/2(2x37+2xl39.59)  2/2(2x37+2x139. 59)^ 

V  2(719.88+2795.66)  '  2(719.88+2795.66)  J 

~  V  7031.08  '  7031.08 )  ~  (0.04543,0.05511). 
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f  yg. 2/2(2x10+2x92. 19)  Xi-o.2/2(2x10+2x92.19)^ 

Component  1:  2(239.96+1984,67)  '  2(239.96+1984.67)  J 


(  178.58  230^28  ^ 

~  1^  4449.26  '  4449.26  J  ~  (0.04014,0.05176). 

(  ;cg.2/2(2xl4+2xl6.15)  Zi-o, 2/2(2x14+2x16. 15)^ 
Component  2;  2(239.96+300.71)  '  2(239.96+300.71)  ) 

-  (  10si!L  >  wkSi  )  -  (0. 04296, 0.0S880). 

(  Ar^.2/2(2xl3+2xl9.97)  ;cf-o. 2/2(2x13+2x19. 97)^ 
Component  3:  2(239.96+402.48)  '  2(239.96+402.48)  ) 


(  50.88  _79i97_^  «  Az-ooiii 

“  1^  1284.88  '  1284.88  J  ~  (0-03960,0.06224). 

The  Bayesian  estimators  of  R(t)  for  t  =  1  hour  are 

Aggregated  Component:  R(t)  =  95.10% 

(  719.88  ^  .  T 
I  2795.66; 


R{t)  = 


X  139.59+37 


(  719.88  >1  ( _ 1 _ ^  T 

V  2795.66;  219b. 


Component  1 :  R{t)  =  95.51% 


R{t)  = 


(  239.96  ^  .  T 
^  1984.67; 


\  92.19+10 


(  239.96  ^  J 
V  1984.67;  1984.67; 


+  1 


Component  2:  R[t)  =  94.58% 

r  96\ 


R{t)  = 


(  239. 
V  300. 


7i; 


+  1 


f2.39.,9_6.-\  ( _ 1 _ A  . 

V  300. 7i;  300. 7i; 


16.15+14 


Component  3;  R(t)  =  95.00% 

(  239.96^  .  T 
I  402.48; 


R(t)  = 


N  19.97+13 


f  239.96^  ; _ 1 _ >1 

V  V  402.48;  402. 48>^ 


+  1 
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99.9882% 


Aggregated  Component  to  System:  -R(t)  = 

R{t)  =  1  -  (1  -  Rc)^  =  1  -  (1  -  95. 10%)^  =  0.999882 
Component  Level  Aggregated  to  System;  R{t)  =  99.9878% 

R{t)  =  1  -  (1  -  0.9551) (1  -  0.9458) (1  -  0.9500) 

The  80%  two-sided  probability  intervals  (TBPI)  are 

Aggregated  Component:  (94.64%  <  A(l)  <  95.56%) 
exp(-0. 05511X1)  <  R(l)  <  exp  (-0 . 04543X1 )  . 

Component  1:  (94.96%  <  ^^(1)  <  96.07%) 
exp(-0. 05176X1)  ^  R(l)  <  exp{-0. 04014X1)  . 

Component  2:  (93.35%  <  R(l)  <  95.80%) 
exp{-0. 06880X1)  <  R{1)  <  exp{-0. 04296X1)  . 

Component  3:  (93.97%  <  R(l)  <  96.12%) 

exp(-0. 06224X1)  <  R{1)  <  exp  (-0. 03960X1)  . 

Aggregated  Component  to  System:  (99. 9846%, 99. 9913%) 
Lower  Bound  =  1  -  (1  -  0.9464)^  =  0.999846 
Upper  Bound  =  1  -  (1  -  0.9556)^  =  0.999913 
Component  Level  Aggregated  to  System; 

(99. 9798%, 99. 9936%) 

Lower  Bound  =  1  -  (1  -  LBi)  (1  -  LBa)  (1  -  LB3) 

=  1  -  (1  -  0.9496) (1  -  0.9335) (1  -  0.9397)  =  0.999798 
Upper  Bound  =  1  -  (1  -  UBi)  (1  -  UBz)  (1  -  UB3) 

=  1  -  (1  -  0.9607) (1  -  0.9580) (1  -  0.9612)  =  0.999936. 


108 


Appendix  6 . 


Simulated  Data  For  Parallel  System  With 


Different  Independent  Components. 


DT  Data  Run  #1 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

6.03 

Cl 

177.73 

C2 

328.02 

Cl 

9.36 

Cl 

180.72 

C2 

328.93 

Cl 

21.91 

Cl 

184.44 

Cl 

337.22 

C2 

24.08 

C3 

187.58 

Cl 

351.00 

Cl 

28.79 

Cl 

189.41 

Cl 

354.85 

C2 

33.23 

C2 

189.68 

Cl 

365.39 

C2 

33.42 

C2 

201.89 

Cl 

367.84 

C2 

38.45 

C2 

209.57 

C2 

371.18 

Cl 

41.87 

C3 

226.86 

Cl 

376.72 

Cl 

43.37 

Cl 

240.59 

Cl 

377.86 

Cl 

49.15 

C3 

243.74 

Cl 

383.68 

C2 

65.41 

C2 

249.69 

C2 

384.25 

C2 

74.71 

C2 

251.30 

Cl 

386.30 

C3 

76.85 

C2 

261.09 

Cl 

390.86 

C2 

102.23 

Cl 

262.25 

Cl 

397.32 

C2 

102.40 

Cl 

264.60 

Cl 

404.24 

Cl 

102.70 

Cl 

267.10 

Cl 

405.30 

Cl 

106.51 

Cl 

274.40 

C2 

418.69 

Cl 

108.57 

Cl 

276.84 

Cl 

424.05 

Cl 

110.18 

Cl 

279.60 

C3 

427.08 

Cl 

133.90 

Cl 

284.25 

Cl 

429.30 

Cl 

136.90 

Cl 

284.26 

C2 

434.15 

Cl 

140.04 

C3 

285.73 

Cl 

437.05 

Cl 

151.32 

Cl 

290.06 

C3 

437.24 

C2 

156.12 

C2 

296.62 

Cl 

443.94 

Cl 

161.19 

C3 

301.54 

Cl 

444.49 

C2 

165.34 

Cl 

304.56 

C2 

446.56 

Cl 

171.86 

Cl 

305.03 

C3 

448.34 

Cl 

172.97 

C2 

306.18 

C3 

457.56 

Cl 

173.20 

Cl 

322.17 

C2 

488.03 

Cl 

177.59 

Cl 

323.68 

C2 
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DT  Data  Run  #2 


Time 

Component 

Failure 

Time 

Component 

Failure 

3.90 

C2 

196.14 

C3 

7.23 

C3 

201.56 

Cl 

7.61 

C2 

203.70 

C2 

7.82 

Cl 

213.24 

Cl 

14.65 

Cl 

217.15 

Cl 

17.42 

Cl 

219.19 

C2 

18.91 

C2 

236.49 

C2 

22.44 

C2 

242.73 

Cl 

34.50 

Cl 

243.28 

Cl 

46.64 

Cl 

249.64 

C2 

46.79 

Cl 

253.23 

C2 

52.25 

Cl 

256.43 

C3 

56.58 

Cl 

257.41 

Cl 

61.16 

Cl 

275.39 

Cl 

84.34 

Cl 

277.64 

Cl 

85.06 

Cl 

283.34 

C2 

85.10 

'C2 

283.80 

Cl 

95.32 

C3 

296.04 

C2 

112.01 

Cl 

296.34 

Cl 

115.58 

C3 

298.74 

C2 

118.78 

Cl 

299.96 

C2 

124.60 

C2 

302.78 

Cl 

130.93 

Cl 

306.89 

Cl 

140.08 

C3 

309.74 

Cl 

156.14 

Cl 

322.13 

Cl 

164.75 

C2 

323.67 

Cl 

166.15 

C2 

324.09 

Cl 

166.29 

Cl 

326.37 

Cl 

177.87 

Cl 

328.53 

Cl 

186.36 

C2 

328.59 

C2 

191.28 

C3 

330.33 

C3 
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DT  Data  Run  #3 


Time 

Component 

Failure 

Time 

Component 

Failure 

2.47 

C2 

274.65 

Cl 

3.89 

C2 

279.05 

C3 

6.14 

C3 

283.69 

Cl 

12.44 

Cl 

288.17 

Cl 

19.74 

C2 

289.58 

Cl 

24.37 

C2 

292.01 

Cl 

30.43 

C2 

296.44 

C2 

37.23 

C2 

306.05 

Cl 

37.56 

Cl 

311.49 

C2 

52.17 

Cl 

313.14 

C3 

54.17 

Cl 

330.60 

C3 

57.07 

Cl 

330.62 

Cl 

58.99 

C3 

332.40 

C3 

64.06 

C2 

343.77 

Cl 

78.26 

Cl 

347.58 

Cl 

78.47 

Cl 

356.96 

Cl 

79.43 

Cl 

363.66 

C3 

100.52 

Cl 

370.35 

Cl 

100.66 

Cl 

371.06 

C2 

107.65 

Cl 

387.09 

C3 

108.07 

Cl 

392.95 

Cl 

112.62 

Cl 

394.49 

Cl 

117.15 

Cl 

394.51 

Cl 

125.21 

C2 

395.85 

Cl 

148.12 

C2 

401.96 

C2 

149.47 

C2 

416.91 

Cl 

152.94 

Cl 

419.24 

Cl 

153.45 

C2 

421.87 

C2 

155.09 

Cl 

425.12 

Cl 

155.69 

C2 

427.82 

Cl 

183.32 

Cl 

434.34 

Cl 

190.35 

C2 

435.53 

C2 

203.84 

C3 

441.33 

Cl 

228.08 

Cl 

446.11 

Cl 

228.64 

C3 

455.95 

Cl 

231.15 

C3 

464.17 

C2 

236.29 

C2 

474.43 

C2 

236.76 

C3 

481.60 

Cl 

245.00 

Cl 

488.36 

C2 

269.60 

C2 

498.39 

Cl 
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DT  Data  Run  #4 


Time 

Component 

Failure 

1.07 

C2 

9.19 

Cl 

10.02 

C2 

11.54 

C3 

11.68 

C2 

14.80 

C2 

17.36 

C2 

33.90 

Cl 

43.33 

Cl 

50.04 

C3 

50.40 

C2 

51.32 

Cl 

51.73 

Cl 

59.03 

Cl 

67.39 

C3 

85.41 

C3 

89.71 

Cl 

95.02 

Cl 

106.10 

Cl 

108.87 

C2 

110.09 

Cl 

113.31 

Cl 

123.58 

Cl 

126.61 

C2 

130.21 

C2 

141.32 

Cl 

144.89 

C2 

157.24 

Cl 

161.29 

C2 

164.69 

Cl 

208.39 

Cl 

208.97 

Cl 

215.92 

Cl 

216.81 

Cl 

234.96 

C2 

239.33 

C2 

254.11 


259.43 


260.10 


266.24 


275.53 


280.58 


283.27 


288.07 


288.67 


296.32 


297.52 


300.02 


300.80 


302.24 


309.28 


323.91 


336.46 


357.20 


381.58 


389.83 


406.69 


411.97 


415.72 


417.64 


435.36 


438.08 


441.96 


445.27 


448.37 


452.14 


454.89 


463.03 


467.94 


Component 

Failure 


Cl 


Cl 


Cl 


C2 


Cl 


C2 


C2 


C2 


Cl 


C2 


C2 


C3 


Cl 


C3 


Cl 


Cl 


C3 


C2 


C2 


C3 


C2 


Cl 


Cl 


Cl 


Cl 


Cl 


Cl 


Cl 


C3 


C2 
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DT  Data  Run  #5 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

0.17 

Cl 

136.39 

Cl 

311.20 

Cl 

11.94 

C3 

139.84 

Cl 

313.51 

Cl 

12.00 

C2 

141.80 

Cl 

319.84 

C2 

15.77 

C2 

151.48 

Cl 

319.96 

C2 

17.13 

Cl 

156.22 

Cl 

320.31 

Cl 

33.77 

C2 

159.43 

Cl 

328.99 

C2 

34.84 

Cl 

166.95 

C2 

334.68 

C2 

37.31 

Cl 

167.06 

Cl 

336.05 

Cl 

43.99 

Cl 

168.03 

Cl 

339.66 

C2 

47.63 

Cl 

168.91 

Cl 

344.45 

C2 

50.27 

C2 

176.16 

Cl 

345.21 

C2 

52.59 

Cl 

180.09 

Cl 

353.47 

Cl 

53.61 

C3 

181.48 

Cl 

364.79 

C3 

59.96 

C2 

198.49 

Cl 

368.50 

Cl 

62.89 

Cl 

202.00 

C2 

371.28 

C2 

67.11 

C2 

211.50 

Cl 

384.22 

C2 

71.75 

Cl 

216.34 

Cl 

400.70 

C2 

83.43 

Cl 

226.76 

C2 

404.40 

Cl 

90.72 

Cl 

235.59 

Cl 

409.07 

C2 

91.24 

Cl 

236.75 

C2 

414.83 

Cl 

93.38 

Cl 

239.91 

Cl 

417.70 

Cl 

103.97 

Cl 

247.21 

Cl 

425.99 

C3 

104.68 

Cl 

256.07 

C2 

426.65 

Cl 

110.20 

C2 

257.26 

Cl 

442.94 

Cl 

110.35 

C3 

257.53 

Cl 

449.73 

Cl 

111.37 

C2 

258.07 

C2 

451.98 

Cl 

112.83 

Cl 

264.20 

C2 

451.99 

Cl 

114.60 

Cl 

271.22 

C2 

452.57 

Cl 

117.70 

C2 

273.63 

Cl 

474.32 

'  C3 

121.11 

Cl 

281.84 

C2 

479.76 

Cl 

123.04 

C3 

282.96 

C3 

479.79 

C2 

126.10 

C2 

283.17 

C2 

480.37 

Cl 

131.07 

C2 

283.24 

Cl 

491.62 

Cl 

133.41 

C3 

310.37 

Cl 

496.36 

Cl 
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OT  Data  Runs  #1  and  #2 


OT  Data  Run  #1 

OT  Data  Run  #2 

Time 

Component 

Failure 

Time 

Component 

Failure 

7.82 

Cl 

5.58 

Cl 

8.11 

C3 

12.85 

C2 

8.27 

C2 

14.31 

C3 

22.23 

C2 

27.61 

C2 

22.99 

C2 

27.78 

Cl 

24.08 

C3 

32.47 

Cl 

27.60 

Cl 

36.89 

Cl 

36.39 

Cl 

47.72 

C2 

36.81 

Cl 

60.70 

Cl 

46.17 

C2 

62.34 

C2 

51.33 

Cl 

70.50 

C2 

51.35 

C2 

74.63 

C2 

66.15 

C2 

85.33 

Cl 

66.33 

C2 

86.59 

C2 

70.65 

Cl 

89.67 

Cl 

72.85 

Cl 

91.98 

C2 

80.58 

Cl 

105.08 

C3 

85.55 

C2 

109.31 

C3 

96.34 

Cl 

113.13 

Cl 

96.94 

C3 

113.29 

C3 

100.62 

Cl 

106.24 

Cl 

108.08 

C2 

112.81 

Cl 

122.22 

Cl 
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Appendix  H.  Calculation  For  Parallel  System  Reliability 
With  Different  Independent  Cos^onents. 

Consider  a  system  composed  of  k  =  3  different 

components  in  parallel  configuration.  The  time  to  failure 
of  each  component  is  assumed  to  be  exponentially 
distributed.  Let  =  0.1,  X2  -  0.05,  and  Xz  =  0.025  be  the 
respective  constant  failure  rates  (failure  per  hour)  of  each 
component . 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  and  3  are 

Ri  =  exp(-O.lOO)  =  90.48% 

R2  =  exp (-0.050)  =  95.12% 

R3  =  exp (-0.025)  =  97.53%. 

The  reliability  of  the  system  can  be  calculated  as 

Rs  =  1  -  (1  -  Ri)  (1  -  R2)  (1  -  R3) 

=  1  -  (1  -  0.9048)  (1  -  0.9512)  (1  -  0.9753)  =  99.9885%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  G) : 


Table  Hi:  Component  1  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

488.03 

55 

0.11270 

2 

498.15 

53 

0.10639 

3 

498.39 

44 

0.08828 

4 

496.55 

36 

0.07250 

5 

496.36 

60 

0.12088 
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Table  H2 :  Component  1  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

122.22 

13 

0.10637 

2 

113.29 

8 

0.07062 

Table  H3:  Component  2  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

488.03 

27 

0.05532 

2 

498.15 

30 

0.06022 

3 

498.39 

24 

0.04816 

4 

496.55 

26 

0.05236 

5 

496.36 

33 

0.06648 

Table  H4 :  Component  2  01 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

A7  =  Xi/Ti 

1 

122.22 

9 

0.07364 

2 

113.29 

8 

0.07062 

Table  H5:  Component  3  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

488.03 

10 

0.02049 

2 

498.15 

8 

0.01606 

3 

498.39 

12 

0.02408 

4 

496.55 

9 

0.01813 

5 

496.36 

9 

0.01813 

Table  H6:  Component  3  01 

Data 

Run  i 

Ti  Hours 

Xi  Failures 

A7  =  Xi/Ti 

1 

122.22 

3 

0.02455 

2 

113.29 

4 

0.03531 

The  Classical  Exponential  Method. 

Assuming  the  time  to  failure  of  each  component  has  the 
exponential  distribution,  point  estimates  for  the  failure 
rate  X  and  80%  confidence  intervals  about  k  will  be 
computed . 


Component 

1  DT: 

A  = 

-  1 

r  248  ^ 

=  0.10010 

ktJ  -  ' 

1  2477.48/ 

Component 

1  OT: 

A  = 

(■^)  -  1 

f  21  ^ 

=  0.08917 

1  r2  A  ~  ' 

L  235.51/ 
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Component  2  DT:  I  =  (-^)  =  (  2477  ^Ts) 
Component  2  OT:  k  -  =  {  235^5r} 


0.05651 


=  0.07218 

Component  3  DT:  I  =  ~  (  2 All  48 )  ~  0.01938 

Component  3  OT:  k  =  “  (  235^  51 )  ~  0.02972 

where  fi  =  248  is  the  number  of  failures  and  Ti  =  2477.48  is 
the  total  test  time  from  Table  HI,  fi  =  21  and  T2  =  235.51 
from  Table  H2,  f  3  =  140  and  T3  =  2477.48  from  Table  H3,  fi  = 
17  and  T4  =  235.51  from  Table  H4,  fb  -  48  and  T5  =  2477.48 
from  Table  H5,  fe  =  1  and  Te  =  235.51  from  Table  H6. 

The  80%  two-sided  confidence  intervals  (TCI)  are 
Component  1  DT:  (0.09205,0.10833) 


( 


^§.2/2(2x248)  ;f!-o.2/2(2x248)>| 

2x2477.48  '  2x2477.48  J 


456.09  536 


4954.96  '  4954 


.77  ^ 
1.96  J 


Component  1  OT:  (0.06532,0.11484) 


( 


^0. 272(2x2 1)  ^1-0.2/2(2x21) 

2x235.51  '  2x235.51  J 


f  30.77  54.09  ^ 

“  V  471.02  '  471.02  J 


Component  2  DT:  (0.05048,0.06271) 


f  ;^§.2/2(2x140) 

^1-0.2/2(2x140) 

1 

^  250.13 

310.72  ^ 

2x2477.48  ' 

2x2477.48  J 

=  1 

i  4954.96  f 

4954.96  J 

Component  2  OT:  (0.05085,0.09533) 

f  ^0.2/2(2x17)  /i-o. 2/2(2x17) ^  ^  23.95  44.9 

2x235.51  f  2x235.51  J  =  ^  471.02  '  471.02  J 

Component  3  DT:  (0.01589,0.02303) 
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(  ;yo. 2/2(2x48)  ;<:i-o.2/2(2x48)^  /  78.73  114.13  ^ 

2x2477.48  '  2x2  All  .AB  J  =  ^  4954.96  '  4954. 96y 

Component  3  OT:  (0.01654,0.04472) 

(  ^0.2/2(2x7)  ;fi-o.2/2(2x7)^  /  7.79  21.06  ^ 

2x235.51  f  2x235.51  /  “  471.02  '  471.02  j 

The  reliability  estimators  of  R{t)  for  t  =  1  hour  are 

Component  1  DT :  Ri {t}  =  exp(-0. 10010X1)  =  90.48% 

Component  1  QT:  Ri{t)  =  exp(-0. 08917X1)  =  91.47% 

Component  2  DT;  Ri  {t)  =  exp (-0.05651X1)  =  94.51% 

Component  2  OT:  R2{t)  =  exp(-0. 07218X1)  =  93.04% 

Component  3  DT:  Rs  it)  =  exp (-0 . 01938X1 )  =  98.08% 

Component  3  OT:  R3{t)  -  exp (-0. 02972X1)  =  97.07% 

Component  Level  Aggregated  to  System  DT;  R{t)  =  99.99% 

R{t)  =  1  -  (1  -  0.9048) (1  -  0.9451) (1  -  0.9808) 
Component  Level  Aggregated  to  System  OT: 

R(t)  =  1  -  (1  -  0.9147)  (1  -  0.9304)  (1  -  0.9707) 

=  99.9826%. 

The  80%  two-sided  confidence  intervals  (TCI)  are 
Component  1  DT:  (89. 73%, 91. 21%) 

exp (-0.10833X1)  <  R(l)  <  exp(-0. 09205X1) 

Component  1  OT;  (89. 15%, 93. 68%) 

exp(-0. 11484X1)  <  R(l)  <  exp(-0. 06532X1) 

Component  2  DT;  (93. 92%, 95. 08%) 
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exp(-0. 06271X1)  <  ^(1)  <  exp  (-0 . 05048X1 ) 

Component  2  OT;  (90. 91%, 95. 04%) 

exp (-0.09533X1)  <  ^(1)  <  exp (-0. 05085X1) 

Component  3  DT;  (97.72%, 98.42%) 

exp(-0. 02303X1)  <  R(l)  <  exp(-0. 01589X1) 

Component  3  OT;  (95. 63%, 98. 36%) 
exp(-0. 04472X1)  <  R{1)  <  exp (-0. 01654X1) 

Component  Level  Aggregated  to  System  DT: 

(99. 9858%, 99. 9932%) 

Lower  Bound  =  1  -  (1  -  0.8973) (1  -  0.9392) (1  -  0.9772) 

Upper  Bound  =  1  -  (1  -  0.9121) (1  -  0.9508) (1  -  0.9842) 

Component  Level  Aggregated  to  System  OT; 

(99. 9569%, 99. 9949%) 

Lower  Bound  =  1  -  (1  -  0.8915) (1  -  0.9091) (1  -  0.9563) 

Upper  Bound  =  1  -  (1  -  0.9368)  (1  -  0.9504)  (1  -  0.9836) . 

The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior. 

For  Component  1  data,  from  Table  HI: 

E{X)  =  T  =  0.10015  and  V(X)  =  sx^  =  0.00038. 


These  yielded  the  parameter  estimates 


a 


0.10015^ 

0.00038 


=  26.40  and 


r  0.10015^ 
V  0.00038>^ 


263.55. 


For  Component  2  data,  from  Table  H3: 
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E{X)  -  A  =  0.05651  and  V{X)  =  Sx^  =  0.00005. 

These  yielded  the  parameter  estimates 

_ f  0.05651^^  _  OT  ft  -  f  0-05651^  _  « 

“  ”  I  0.00005  )  ~  and  P  -  O.OOOOsJ  1130.2. 

For  Component  3  data,  from  Table  H5: 

E(X)  =  T  =  0.01938  and  V(X)  =  =  0.00001. 


These  yielded  the  parameter  estimates 


a 


f  0.01938^^ 

V  0.00001  J 


=  37.56  and 


0.01938^ 

o.ooooiy 


=  1938. 


Bayesian  Exponential  Posterior  Analysis. 

Component  ,1 ;  1  =  (  sssTsiIMtsF )  “  0.09496. 

CQmponent  2 :  2  =  (  23^11^^)  “  0.05922. 

Component,,  3:  2  =  (■2357^+1938')  “  0.02050. 

The  80%  two-sided  probability  intervals  (TBPI)  are 

Comppnent  1;  2(235.51+263.55)  '  2(235.51+263.55)  J 

“  998.12  '  998.12  J  “  (0.07706,0.11216). 

f  ;irg.2/2(2xl7+2x63.87)  X  1-0.2/ 2(2x11+2x63 .81)  \ 

Cpmp,o,nent„  2 :  2(235.51+1130.2)  '  2(235.51+1130.2)  J 

”  1^  2731.42  '  2731. 42y  “  (0.05070,0.06750). 

f  X^.2/2(^X^+1‘X^1 .36i)  ^ 

CpmppPent,  3;  2(235.51+1938)  '  2(235.51+1938)  ) 

(  72.39  106.47  ^  a,  a  aa^i-a^ 

~  1^  4347.02  '  4347. 02y  ~  (0.01665,0.02449). 
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The  Bayesian  estimators  of  R(t)  for  t  =  1  hour  are 
Component  1 ;  R{t)  =  90.95% 


^(t) 


235.51 

263.55 


235.51 

263.55 


263.55. 


+  1 


26.4+21 


Component  2:  R[t)  =  94.25% 


R[t) 


( 


235.51^  ,  1 
1130.27 


235.51 

1130.2 


1 

1130.2 


63.87+17 


Component  3 :  R{t)  =  97.97% 


^(t) 


(  235. 51> 
L  1938  J 

+  1 

r  235 

.51\  . 

f-  1  -"1 

i+ij 

^  1938  J  ' 

1  1938J 

37.56+7 


Component  Level  Aggregated  to  System;  R(t)  =  99.9894% 


R(t)  =  1  -  (1  -  0.9095) (1  -  0.9425) (1  -  0.9797). 


The  80%  two-sided  probability  intervals  (TBPI)  are 
Component  1:  (89.39%  <  R[l)  <  92.58%) 
exp(-0. 11216X1)  <  R(l)  <  exp  (-0. 07706X1)  . 

Component  2:  (93.47%  <  R{1)  <  95.06%) 
exp(-0. 06750X1)  <  R{1)  ^  exp (-0.05070X1)  . 

Component  3;  (97.58%  <  R(l)  ^  98.35%) 

exp(-0. 02449X1)  ^  R{1)  <  exp  (-0. 01665X1)  . 

Component  Level  Aggregated  to  System; (99. 9832%, 99. 994%) 
Lower  Bound  =  1  -  (1  -  0.8939) (1  -  0.9347) (1  -  0.9758) 
Upper  Bound  =  1  -  (1  -  0.9258) (1  -  0.9506) (1  -  0.9835). 
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Appendix  I . 


Simulated  Data  For  A  Bridge  System  With 


Independent  And  Identical  Cosponents. 

DT  Data  Run  #1 


Time 

Fail 

Time 

Fail 

Time 

Fail 

Time 

Fail 

6 . 66 

C5 

139.79 

C2 

260.39 

C3 

384.83 

C5 

12.04 

C3 

140.33 

C2 

264.41 

C2 

386.05 

C2 

12.06 

Cl 

142.74 

C3 

266.69 

C2 

391.94 

C4 

20.96 

Cl 

149.49 

C5 

270.53 

Cl 

394.66 

C3 

25.80 

C3 

158.44 

Cl 

271.10 

Cl 

401.82 

C4 

25.99 

C3 

164.77 

C2 

271.47 

C3 

403.74 

Cl 

29.44 

C4 

167.16 

C4 

275.73 

C5 

407.75 

C3 

31.02 

C3 

171.06 

C2 

277.59 

C4 

407.96 

C4 

31.76 

C5 

182.30 

C4 

277.71 

Cl 

420.43 

C3 

33.23 

C2 

185.90 

Cl 

279.82 

Cl 

424.30 

C4 

34.66 

C3 

188.21 

Cl 

282.18 

C5 

428.60 

C3 

43.96 

C3 

189.61 

C5 

292.90 

C5 

430.08 

C5 

46.09 

C3 

192.67 

C3 

298.96 

C5 

433.07 

C3 

50.10 

Cl 

192.92 

Cl 

303.40 

C5 

437.72 

C5 

50.43 

Cl 

194.60 

C5 

306.61 

Cl 

441.39 

Cl 

51.04 

Cl 

195.77 

C2 

312.40 

Cl 

445.33 

C3 

56.41 

C4 

201.01 

C2 

313.09 

C5 

447.24 

C2 

58.66 

Cl 

201.87 

C4 

317.51 

C4 

447.41 

C5 

60.53 

C4 

202.79 

Cl 

319.45 

C2 

450.87 

C5 

61.88 

Cl 

203.96 

C2 

320.35 

C5 

464.21 

C4 

67.88 

Cl 

209.42 

C5 

322.77 

C4 

466.10 

C3 

78.68 

C2 

211.44 

C2 

323.89 

C5 

466.32 

C3 

89.25 

C2 

212.15 

C3 

326.18 

Cl 

468.33 

C5 

96.73 

Cl 

219.26 

C5 

332.49 

C3 

485.42 

C4 

107.96 

C4 

219.83 

C5 

336.40 

C3 

486.10 

C4 

113.57 

Cl 

222.17 

C4 

341.22 

C4 

490.04 

Cl 

117.30 

C2 

223.68 

C4 

354.70 

Cl 

490.37 

Cl 

119.98 

C2 

224.57 

Cl 

360.58 

C3 

491.49 

C4 

124.74 

C2 

226.39 

Cl 

361.05 

Cl 

492.86 

C2 

125.37 

C3 

229.76 

C3 

365.46 

C4 

494.52 

C3 

126.60 

Cl 

237.22 

C4 

372.13 

C3 

495.03 

C4 

129.59 

Cl 

247.40 

C3 

380.68 

Cl 

495.29 

C3 

133.52 

C2 

257.93 

C3 

381.88 

C3 

495.50 

C3 

139.08 

C3 

259.89 

C5 

383.66 

C4 

498.06 

C3 

499.37 

C3 

DT  Data  Run  #2 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

4.36 

C2 

131.29 

C2 

344.14 

C5 

4.51 

Cl 

134.68 

C3 

345.53 

Cl 

6.77 

C4 

146.27 

C4 

351.96 

C2 

7.22 

Cl 

147.67 

C5 

353.05 

C3 

8.43 

Cl 

150.14 

C5 

359.08 

C2 

12.33 

C3 

151.84 

Cl 

368.48 

Cl 

12.69 

C5 

157.55 

C4 

372.65 

C5 

18.40 

C5 

162.11 

C5 

372.68 

C2 

19.66 

C4 

165.01 

C4 

379.50 

C4 

20.53 

C3 

165.05 

Cl 

383.91 

C3 

21.38 

C3 

166.16 

C3 

387.88 

C2 

22.73 

C4 

167.51 

Cl 

390.46 

C5 

25.95 

C3 

168.83 

C2 

398.52 

C5 

27.04 

C4 

173.24 

C5 

409.40 

Cl 

29.43 

C2 

174.81 

C2 

413.68 

C3 

37.07 

Cl 

176.95 

C2 

416.10 

C2 

43.19 

C5 

195.96 

Cl 

416.42 

C4 

50.09 

C3 

196.87 

C4 

418.07 

Cl 

52.26 

C5 

198.22 

C4 

418.60 

C5 

64.90 

C4 

199.84 

C3 

424.48 

Cl 

76.03 

C5 

202.09 

Cl 

426.34 

C4 

79.50 

C4 

215.72 

C5 

427.27 

C2 

83.22 

C3 

228.84 

C5 

433.87 

C3 

86.45 

Cl 

228.97 

C3 

436.88 

C5 

89.05 

Cl 

229.03 

Cl 

437.44 

C3 

90.71 

C3 

229.58 

Cl 

439.65 

C2 

94.35 

C5 

230.82 

Cl 

442.77 

C2 

95.79 

Cl 

241.36 

Cl 

443.83 

C2 

95.87 

C5 

271.43 

Cl 

443.93 

C3 

96.05 

C3 

274.81 

C5 

446.19 

C4 

97.80 

C5 

282.22 

Cl 

451.99 

C4 

109.38 

C3 

282.96 

Cl 

456.79 

C4 

110.57 

C4 

287.80 

C4 

457.42 

C4 

111.47 

Cl 

289.59 

C4 

467.93 

C3 

117.46 

C3 

290.93 

C5 

473.89 

Cl 

117.48 

C2 

308.02 

C5 

477.44 

C3 

119.12 

Cl 

310.90 

C3 

477.97 

C5 

119.12 

C4 

313.61 

Cl 

481.49 

C4 

123.30 

C4 

335.89 

C3 

486.24 

C4 

124.37 

C2 

337.22 

C4 

490.82 

C5 

126.57 

C5 

337.25 

C4 

497.27 

C5 
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DT  Data  Run  #3 


2.96 


4.04 


5.75 


12.68 


16.24 


33.59 


37.24 


38.44 


48.96 


55.34 


68.16 


69.77 


79.30 


81.39 


85.81 


86.40 


88.56 


90.60 


97.14 


97.65 


103.04 


103.96 


137.42 


140.36 


147.29 


148.69 


149.72 


149.80 


152.45 


153.95 


163.35 


164.02 


166.31 


173.01 


178.30 


186.75 


192.89 


196.39 


198.78 


201.63 


203.02 


Component 

Failure 


C2 


C2 


C4 


C2 


C2 


Cl 


Cl 


Cl 


C2 


C3 


C2 


C2 


C4 


Cl 


C4 


C2 


C4 


C5 


C5 


C3 


C4 


C3 


Cl 


C5 


Cl 


C3 


C2 


C3 


C3 


Cl 


C4 


Cl 


C3 


C2 


C2 


C3 


Cl 


Cl 


Cl 


C5 


Cl 


209.14 


215.08 


220.14 


221.98 


227.31 


231.32 


238.67 


240.93 


241.09 


241.26 


244.01 


244.84 


245.51 


247.55 


262.06 


263.28 


267.35 


268.55 


268.64 


273.88 


275.38 


275.81 


277.54 


278.22 


278.93 


279.98 


286.78 


291.53 


300.53 


301.34 


310.11 


326.98 


327.32 


329.63 


333.29 


336.91 


340.68 


350.62 


354.13 


358.68 


363.56 


Component 

Failure 


C3 


C4 


C3 


Cl 


C3 


C5 


C2 


C5 


C5 


C2 


C4 


C2 


Cl 


C4 


C5 


C5 


C3 


C4 


C5 


C5 


Cl 


C2 


C4 


Cl 


C2 


C3 


C2 


C5 


C3 


C5 


C2 


C3 


C3 


C2 


C3 


C5 


C5 


C4 


C3 


C5 


C2 


Time 

Component 

Failure 

363.63 

C5 

370.63 

C3 

370.89 

C5 

376.93 

C2 

380.55 

C3 

380.57 

C5 

385.48 

Cl 

386.03 

C4 

387.73 

C5 

397.54 

C2 

403.20 

Cl 

404.23 

Cl 

408.53 

Cl 

408.92 

C3 

410.34 

C3 

411.08 

C5 

412.12 

C2 

412.25 

C5 

418.46 

C2 

418.58 

C5 

422.00 

C2 

426.65 

C3 

429.11 

C4 

429.70 

Cl 

430.40 

C2 

431.59 

C5 

434.30 

C4 

434.67 

Cl 

438.59 

Cl 

441.21 

C4 

450.70 

C4 

457.12 

C4 

457.21 

C3 

457.94 

Cl 

459.89 

Cl 

461.63 

Cl 

466.28 

C2 

472.38 

C4 

474.14 

C2 

475.16 

C4 

476.14 

Cl 

492.26 

C3 
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DT  Data  Run  #4 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

9.99 

C4 

179.08 

C5 

341.36 

C3 

14.23 

C5 

182.59 

C2 

349.36 

C5 

14.24 

C5 

184.75 

C3 

351.70 

Cl 

23.47 

C3 

186.46 

C4 

357.70 

C4 

24.44 

C5 

192.17 

C3 

366.46 

C4 

25.32 

Cl 

192.27 

C2 

371.11 

C2 

28.88 

C5 

196.37 

Cl 

382.98 

C3 

30.50 

Cl 

209.78 

C3 

384.09 

C5 

32.08 

C3 

218.92 

C2 

390.94 

C5 

32.86 

C4 

222.44 

C2 

391.66 

Cl 

38.49 

C2 

224.64 

C5 

394.03 

C3 

40.20 

C4 

226.48 

Cl 

397.11 

C2 

41.35 

C4 

227.60 

C4 

397.88 

C2 

41.83 

C4 

231.42 

Cl 

400.00 

C4 

42.61 

C2 

233.02 

C5 

406.06 

C3 

44.88 

Cl 

243.07 

Cl 

419.63 

C5 

47.21 

C2 

243.08 

Cl 

419.94 

C3 

59.71 

Cl 

247.54 

C3 

425.22 

C3 

63.46 

C5 

249.07 

C3 

428.58 

Cl 

64.38 

Cl 

259.52 

C3 

432.91 

C2 

80.78 

C2 

260.99 

Cl 

434.18 

Cl 

84.46 

C2 

262.70 

C5 

435.45 

C5 

91.49 

C4 

276.34 

C4 

436.52 

Cl 

99.18 

C2 

286.19 

C2 

436.58 

C5 

103.74 

C2 

286.85 

C2 

447.31 

C4 

108.38 

C4 

291.65 

Cl 

448.34 

C3 

109.49 

Cl 

292.46 

C2 

452.06 

Cl 

109.73 

C2 

292.62 

C3 

453.11 

Cl 

112.68 

C5 

292.92 

C4 

457.29 

C4 

118.74 

Cl 

298.54 

C5 

457.46 

C4 

124.77 

C5 

298.93 

C5 

462.20 

C2 

129.61 

Cl 

300.74 

Cl 

462.78 

C4 

130.63 

C2 

302.56 

C2 

464.59 

C5 

131.63 

C2 

302.72 

C2 

478.08 

C5 

140.39 

C3 

303.75 

C2 

480.56 

C5 

146.45 

C5 

311.16 

Cl 

486.78 

Cl 

151.21 

C3 

319.74 

C3 

487.90 

C3 

158.31 

C4 

320.36 

C2 

491.72 

C4 

163.06 

C5 

322.40 

C3 

493.38 

C3 

166.08 

C4 

326.34 

C4 

493.77 

C4 

176.23 

C2 

327.79 

C5 

498.61 

C5 

177.60 

Cl 

337.68 

C5 
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DT  Data  Run  #5 


Time 

Fail 

Time 

Fail 

Time 

Fail 

0.66 

Cl 

100.21 

C5 

230.54 

Cl 

2.46 

Cl 

101.18 

C3 

237.16 

C3 

9.64 

C3 

107.58 

Cl 

250.45 

C3 

11.34 

C2 

108.34 

C3 

259.24 

C5 

15.24 

C5 

110.96 

Cl 

261.30 

C3 

15.76 

C2 

113.82 

C5 

269.72 

C2 

16.99 

Cl  . 

116.24 

C5 

271.63 

C3 

18.92 

C4 

136.24 

C4 

283.72 

Cl 

19.35 

Cl 

136.89 

Cl 

284.86 

C3 

23.96 

Cl 

137.38 

C3 

287.11 

C3 

27.71 

C5  1 

143.38 

C2 

290.48 

Cl 

32.46 

C3 

144.89 

C2 

293.61 

Cl 

39.19 

C4 

146.83 

C2 

298.03 

C5 

41.84 

C4 

148.11 

C2 

298.94 

C4 

42.73 

C3 

148.17 

C4 

303.67 

C5 

46.14 

Cl 

150.91 

C2 

304.29 

C3 

51.54 

Cl 

153.87 

C4 

304.86 

C2 

53.34 

C4 

157.99 

Cl 

312.43 

C3 

54.58 

C5 

159.91 

C5 

312.70 

C5 

58.58 

C2 

165.49 

C3 

318.07 

Cl 

59.05 

Cl 

166.48 

C5 

330.95 

C2 

62.97 

C4 

170.87 

C4 

339.34 

C4 

63.89 

C3 

177.89 

C3 

343.35 

C4 

70.26 

C3 

181.20 

C3 

344.37 

C5 

70.93 

C2 

182.23 

C2 

345.96 

C3 

83.10 

Cl 

184.71 

C3 

350.78 

C5 

84.92 

C4 

199.68 

Cl 

351.37 

C4 

88.38 

C3 

201.85 

C3 

351.86 

Cl 

89.47 

C2 

214.69 

C4 

352.32 

C4 

89.86 

Cl 

218.13 

C4 

354.27 

C4 

98.77 

C2 

221.07 

C2 

354.44 

Cl 

98.82 

C3 

223.39 

Cl 

358.72 

C3 

99.94 

C4 

225.80 

C5 

360.41 

Cl 

Time 


361.06 


362.47 


362.49 


366.30 


372.01 


379.40 


387.26 


387.52 


394.09 


398.00 


402.96 


408.64 


412.57 


413.73 


418.63 


424.71 


429.67 


429.92 


432.24 


432.71 


449.61 


457.91 


462.61 


462.62 


465.37 


469.82 


485.61 


491.02 


492.67 


496.42 


499.68 


Fail 


C3 


C4 


C3 


C2 


C3 


C3 


C3 


C2 


Cl 


C2 


C2 


Cl 


C4 


C3 


C4 


C2 


C4 


C2 


C3 


C4 


C5 


C5 


C5 


C3 


Cl 


C5 


C5 


C4 


C2 


C5 


C5 
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OT  Data 

,  Run  #1 

Time 

Component 

Failure 

4.06 

C3 

8.27 

C2 

12.04 

C3 

12.80 

C3 

15.64 

Cl 

22.23 

C2 

35.98 

C3 

36.82 

C3 

39.57 

C5 

39.80 

C2 

44.76 

C5 

52.07 

Cl 

52.25 

Cl 

59.55 

C5 

62.05 

C4 

63.94 

C5 

65.86 

C3 

71.46 

Cl 

77.51 

C4 

78.45 

C2 

80.02 

Cl 

88.39 

C3 

89.69 

C2 

93.16 

Cl 

95.46 

C5 

104.13 

Cl 

108.31 

C5 

108.52 

C2 

114.59 

C4 

123.77  C3 


OT  Data 

Run  #2 

Time 

Component 

Failure 

0.87 

C2 

8.85 

C3 

9.38 

Cl 

15.50 

C2 

17.54 

Cl 

19.63 

C2 

28.30 

C2 

29.50 

Cl 

33.68 

C2 

40.37 

C4 

42.48 

C4 

44.47 

C4 

47.61 

C5 

58.12 

C3 

69.32 

C5 

76.44 

Cl  .  1 

76.90 

C5 

79.17 

C3 

87.49 

C5 

90.62 

C5 

98.37 

C3 

103.75 

Cl 

116.16 

C3 

120.20 

C4 

120.63 

C3 

120.82 

Cl 

122.92 

Cl 

Appendix  J.  Calculation  For  A  Bridge  System  Reliability 
With  Independent  And  Identical  Conponents. 

Consider  a  system  composed  of  5  independent  and 
identical  components  given  in  Figure  1  (page  30) .  The  time 
to  failure  of  each  component  is  assumed  to  be  exponentially 
distributed.  Let  X  =  0.05  failure  per  hour  be  the 

respective  constant  failure  rates  of  each  component. 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  3,  4,  and  5  are 
Rc  =  Ri  =  R2  =  Ra  =  R4  =  R5  =  exp(-0.05)  =  95.12%. 

The  reliability  of  the  system  can  be  calculated  as 
Rs  =  R1R4  +  R2R5  +  R1R3R5  +  R2R3R4  ~  R1R2R3R4  “  R1R2R3R5  “ 
R1R2R4R5  “  R1R3R4R5  “  R2R3R4R5  2R1R2R3R4R5  ~  99.50%. 

Using  a  computer  program  called  RAPTOR,  which  is 
developed  by  AFOTEC,  to  simulate  the  true  system,  the 
following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  I) ; 


Table  Jl:  Aggregated  Components  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

2,496.85 

136 

0.05447 

2 

2,486.35 

123 

0.04947 

3 

2,461.3 

124 

0.05038 

4 

2,493.05 

125 

0.05014 

5 

2,498.4 

130 

0.05203 

Table  J2;  Aggregated  Components  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

618.85 

30 

0.04848 

2 

604.1 

26 

0.04304 
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Table  J3:  Component  1  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

499.37 

32 

0.06408 

2 

497.27 

28 

0.05631 

3 

492.26 

27 

0.05485 

4 

498.61 

26 

0.05215 

5 

499.68 

27 

0.05404 

Table  J4:  Component  1  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

123.77 

7 

0.05656 

2 

120.82 

6 

0.04966 

Table  J5;  Component  2  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

499.37 

21 

0.04205 

2 

497.27 

17 

0.03419 

3 

492.26 

28 

0.05688 

4 

498.61 

27 

0.05415 

5 

499.68 

23 

0.04603 

Table  J6:  Component  2  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

123.77 

6 

0.04848 

2 

120.82 

5 

0.04138 

Table  J7 ;  Component  3  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

499.37 

35 

0.07009 

2 

497.27 

24 

0.04826 

3 

492.26 

25 

0.05079 

4 

498.61 

22 

0.04412 

5 

499.68 

33 

0.06604 

Table  J8:  Component  3  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

123.77 

8 

0.06464 

2 

120.82 

6 

0.04966 
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Table  J9:  Component  4  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

499.37 

25 

0.05006 

2 

497.27 

27 

0.05430 

3 

492.26 

20 

0.04063 

4 

498.61 

23 

0.04613 

5 

499.68 

25 

0.05003 

Table  JIO:  Com] 

ponent  4  OT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

Ai  =  Xi/Ti 

1 

123.77 

3 

0.02424 

2 

120.82 

4 

0.03311 

Table  Jll:  Com] 

ponent  5  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

499.37 

23 

0.04606 

2 

497.27 

27 

0.05430 

3 

492.26 

24 

0.04876 

4 

498.61 

27 

0.05415 

5 

499.68 

22 

0.04403 

Table  J12:  Com] 

ponent  5  OT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

123.77 

6 

0.04848 

2 

120.82 

5 

0.04138 

The  Classical  Exponential  Method. 

Assuming  the  time  to  failure  of  each  component  has  the 
exponential  distribution,  point  estimates  for  the  failure 
rate  X  and  80%  confidence  intervals  about  X  will  be 
computed. 

Component  Aggregated  DT:  1  =  (  12/3^5^  95)  0.05130 

Component  Aggregated  OT:  2  -  (  122^2^  95)  0.04579 
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Component 

*1 

DT: 

X 

_  / 

^  140  ^ 

0.05629 

1 

1  T3J 

-  \ 

<  2487.197 

Component 

•1 

OT: 

X 

(h.) 

^  13  ^ 

0.05315 

1 

\  T4J 

<  244.597 

Component 

0 

DT: 

X 

_  / 

^  116 

0.04664 

z 

[  TsJ 

-  \ 

<  2487.197 

Component 

2 

OT: 

X 

(II) 

"  11  ^ 

0.04497 

1  TeJ 

^  244.597 

Component 

3 

DT: 

X 

fh) 

_  1 

^  139  ^ 

0.05589 

-  \ 

<  2487.197 

Component 

3 

OT: 

1 

(II) 

_  1 

^  14  ^ 

0.05724 

-  \ 

244.597 

Component 

4 

DT: 

1 

— 

(II) 

_  / 

^  120  ^ 

0.04825 

{  TsJ 

-  \ 

.  2487.197 

Component 

4 

OT: 

X 

=  1 

f  fio' 

]  _ 

r  7  "j 

0.02862 

{  Tio. 

1 244.597 

Component 

5 

DT: 

X 

=  1 

f  fii' 
<  Tn. 

)- 

(*9) 

= 

0.04945 

Component 

5 

OT: 

X 

=  1 

f  fx2'' 
1  Ti2. 

)  = 

( 244.59) 

0.04497 

where  fi  -  638  is  the  number  of  failures  and  Ti  =  12435.95 
is  the  total  test  time  from  Table  Jl,  fz  =  56  and  Tz  = 
1222.95  from  Table  J2,  fa  =  140  and  T3  =  2487.19  from  Table 
J3,  fi  =  13  and  T4  =  244.59  from  Table  J4,  fs  =  116  and  T5  = 
2487.19  from  Table  J5,  £&  =  11  and  Te  =  244.59  from  Table 

J6,  fi  =  139  and  T7  =  2487.19  from  Table  J7,  fe  =  14  and  Ts  = 
244.59  from  Table  J8,  jfg  =  120  and  T9  =  2487.19  from  Table 
J9,  jfio  =  7  and  Tio  =  244.59  from  Table  JIO,  fn  =  123  and  Tn 
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=  2487.19  from  Table  Jll,  fiz  =  11  and  T12  =  244.59 
Table  J12. 

The  80%  two-sided  confidence  intervals  (TCI)  are 
Component  Aggregated  DT:  (0.04872,0.05392) 

f  7o.2/2(^^638)  ;fi-o.2/2(2^^38) ^  ^  1211.7  1341. 15'\ 

y  2x12435.95  '  2x12435.95  J  =  24871.9  '  24871.9  J 

Component  Aggregated  OT;  (0.03815,0.05379) 

f  ^0. 272(^^56)  A(i-o. 272(^^56) ^  ^  93.3  131. 56^ 

^  2x1222.95  '  2x1222.95  J  =  ^  2445.9  '  2445.9  J 


Component  1  DT :  (0.05028,0.06247) 


( 


/0.272(^^^^0)  ^i-0.272(^^^^0)  A 

2X2487.19  '  2x2487.19  J 


=  ( 


250.13 


310 


4974.38  '  4974 


.72  ^ 
i.38  J 


Component  1  OT:  (0.03535,0.07270) 


( 


;^§.272(2X13)  ;y!-o.272(2Xl3)>| 

2x244.59  '  2x244.59  J 


(  17.29  35.56  ^ 

~  489.18  '  489.18  J 


Component  2  DT:  (0.04118,0.05227) 

Xo. 272(2x116)  272(2x116) f  204.86 


( 


2x2487.19  '  2x2487.19  )  = 


4974.38  '  4974 

Component  2  OT;  (0.02870,0.06299) 

(  yp. 272(2x11)  1-0. 272 (2x11)  f  14.04  30.81  ^ 

2x244.59  '  2x244.59  )  =  489.18  '  489.18  ) 


260  ^ 
)74.38  j 


Component  3  DT:  (0.04990,0.06204) 

■272(2x139)  j(i-o. 272(2x139) ^  f  248.24 


(4 


2x2487.19  '  2X2487.19  7  “  ( 


308 


4974.38  '  4974 

Component  3  OT:  (0.03872,0.07751) 


.61  ^ 

.38  J 


from 
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Component  4  DT :  (0.04270,0.05397) 


2x2487.19  '  2x2487.19 


Component  4  OT:  (0.01592,0.04306) 


1 

'  212.39 

268.47  A 

I  =  1 

k  4974.38  ' 

4974.38  j 

^0.2/2(^^^) 

2x244.59 


;yi-o.2/2(2x7)'\  /  7  79  _21,;_^^ 

.59  f  2x244.59  j  ~  \  489.18  '  489.18  y 


Component  5  DT;  (0.04383,0.05525) 


/o.2/2(^^^^3)  ;ri^0.2/2(^ 

2x2487.19  '  2X2487 


_  (  218.04  274.82  ^ 

.19  y  -  4974.38  '  4974.38  y 


Component  5  OT;  (0.02870,0.06912) 

/  ..2  /Ovii\  „2  z' 


>f0.2/2(^^^^)  1-0. 2/2 

2X244.59  '  2x244.59 


^  _  r  14.04  33.81  ^ 

)  ~  \  489.18  '  489.18  y  * 


The  reliability  estimators  of  R(t)  for  t  =  1  hour  are 
Component  Aggregated  DT; 

^(t)  =  exp(-0. 05130X1)  =  95.00% 

Component  Aggregated  OT: 

^(t)  =  exp  (-0.04  57  9X1)  =  95.52% 

Component  1  DT;  ^(t)  =  exp (-0.05629X1)  =  94.53% 

Component  1  OT;  ^ ( t)  =  exp (-0 . 05315X1 )  =  94.82% 

Component  2  DT:  ^ ( t)  =  exp(-0. 04664X1)  =  95.44% 

Component  2  OT:  %{t)  =  exp  (-0. 04497X1 )  =  95.60% 

Component  3  DT:  ^(t)  =  exp (-0. 05589X1 )  =  94.57% 
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94.44% 


Component  3  OT;  R3  [t)  =  exp (-0 . 05724X1 )  = 

Component  4  DT :  R^it)  =  exp ( -0 . 04825X1 )  =  95.29% 

Component  4  OT;  Ri (t)  -  exp (-0. 02862X1)  =  97.18% 

Component  5  DT:  Rsit)  =  exp (-0 . 04945X1 )  =  95.18% 

Component  5  OT:  R5  {t)  =  exp (-0 . 04497X1 )  =  95.60% 

Substitute  the  appropriate  values  into  -R(t)  =  R1R4  +  R2R5 
+  R1R3R5  +  R2R3R4  ~  R1R2R3R4  ”  R1R2R3R5  ”  R1R2R4R5  “  R1R3R4R5  “ 
R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 

Aggregated  Component  to  System  DT:  R(t)  =  99.48% 
Aggregated  Component  to  System  OT:  R{t)  =  99.58% 
Component  Level  Aggregated  to  System  DT:R(t)  =  99.50% 

Component  Level  Aggregated  to  System  OTiRlt)  =  99.63%. 
The  80%  two-sided  confidence  intervals  (TCI)  are 
Aggregated  Component  DT:  (94 . 75%, 95.25%) 

exp(-0. 05392X1)  <  R{1)  <  exp (-0 . 04872X1) 

Aggregated  Component  OT:  (94. 76%, 96. 26%) 
exp(-0. 05379X1)  <  R{1)  ^  exp  (-0. 03815X1) 

Component  1  DT:  (93. 95% , 95.10%) 

exp (-0.06247X1)  ^  R(l)  ^  exp  (-0.05028X1) 

Component  1  OT:  (92 . 99% , 96 .53%) 

exp(-0. 07270X1)  ^  R{1)  <  exp  (-0. 03535X1) 

Component  2  DT:  (94. 91%, 95. 97%) 
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exp(-0. 05227X1)  <  ^(1)  <  exp  (-0 . 04118X1) 

Component  2  OT:  (93. 90%, 97. 17%) 

exp(-0. 06299X1)  <  R(l)  <  exp  (-0 . 02870X1 ) 

Component  3  DT:  (93. 99%, 95. 13%) 

exp(-0. 06204X1)  <  R{1)  <  exp(-0. 04990X1) 

Component  3  OT;  (92. 54%, 96. 20%) 

exp(-0. 07751X1)  <  ^(1)  <  exp(-0. 03872X1) 

Component  4  DT:  (94 . 75% , 95 . 82%) 

exp  (-0.05397X1)  <  :R(1)  <  exp  (-0. 04270X1) 

Component  4  OT:  (95. 79%, 98. 42%) 

exp(-0. 04306X1)  ^  ^(1)  <  exp(-0. 01592X1) 

Component  5  DT:  (94. €3%, 95. 71%) 

exp(-0. 05525X1)  <  ^(1)  <  exp  (-0 . 04383X1) 

Component  5  OT:  (93. 32%, 97. 17%) 

exp(-0. 06912X1)  <  ^(1)  <  exp(-0. 02870X1) 

Substitute  the  appropriate  values  into  R(t)  =  R1R4  +  R2R5 
t  R1R3R5  t  R2R3R4  ~  R1R2R3R4  ~  R1R2R3R5  ~  R1R2R4R5  —  R1R3R4R5  — 
R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 

Aggregated  Component  to  System  DT:  (99. 42%, 99. 53%) 
Aggregated  Component  to  System  OT:  (99. 43%, 99. 71%) 
Component  Level  Aggregated  to  System  DT: 

(99.38%, 99. 61%) 
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Component  Level  Aggregated  to  System  OT: 

(99. 24%, 99. 85%) . 

The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior. 

For  the  component  aggregated  data,  from  Table  Jl: 


E{X)  =  X  =  0.051298  and  V(X)  =  Sx^  =  0.000004. 


These  yielded  the  parameter  estimates 


a 


0.051298^^ 
0.000004  J 


657.87  and  p  = 


f  0.051298^ 
V.  0.000004; 


12824.5. 


For  Component  1  data,  from  Table  J3: 

E{k)  =  T  =  0.05628  and  V[X)  =  Sx^  =  0.00002. 


These  yielded  the  parameter  estimates 

“  -  (  °o'.folf2  )  “  150.85  and  P  -  -  2680.14. 

For  Component  2  data,  from  Table  J5: 


E(X)  =  A  =  0.04666  and  V(>.)  =  sx'  =  0.00008. 


These  yielded  the  parameter  estimates 

f  0,04666^)  ^  ^  ;  0.04666^  ^  -5- 

®  V  0.00008  ;  25.92  and  p  V  0.00008;  555.48. 

For  Component  3  data,  from  Table  J7 : 

E(X)  =  T  =  0.05586  and  V(l}  =  sx'  =  0.00013. 


These  yielded  the  parameter  estimates 


a 


0.05586^^ 
0.00013  ; 


=  23.64  and 


0.05586A 

0.00013; 


423.19. 


For  Component  4  data,  from  Table  J9: 
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E{X)  =  I  =  0.04823  and  V Ck)  =  =  0.00003. 


These  yielded  the  parameter  estimates 


a 


(  0.04823^^  _ 
V  0.00003  )  - 


89.47  and 


0.04823^ 

0.000037 


=  1855. 


For  Component  5  data,  from  Table  Jll: 


E['k)  =  k  =  0.04946  and  V(^)  =  Sx^  =  0.00002. 


These  yielded  the  parameter  estimates 

“  -(  o'.ZVof)  -  and  P  =  ( I'-gb-Oof)  =  2248-09. 


Bayesian  Exponential  Posterior  Analysis, 

Aggregated  Component:  2=  (  1222^9^+12824  s)  *  0.05082. 


Component 

J,: 

1  =  1 

(  13  +  150.85 

0.05602. 

L  244. 59  +  2680. 14>^ 

Component 

J.: 

1  =  1 

(  11  +  25.92  >1 

0.04614. 

1  244.59  +  555.487 

Component 

A' 

1  =  1 

r  14+23.64  ^ 

0.05637. 

L  244.59  +  423.197 

Component 

J.: 

=  1 

r  7+89.47  ^ 

0.04595. 

^  244.59  +  18557 

Component 

_5: 

^  =  1 

<  11+111.19  ^ 

0.04902. 

^  244.59  +  2248.097 

The  80%  two-sided  probability  intervals  (TBPI)  are 
Aggregated  Component: 

f  0.2/2(2x56+2x6 57. 8 7)  /f_o. 2/2(2x56+2x657 . 87) \ 

V  2(1222.95+12824.5)  '  2(1222.95+12824.5)  ) 

(  1359.7  1496.64  ^ 

~  1^  28094.9  '  28094.9  )  ~  (0*04840,0 . 05327)  . 
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(  /o. 2/2(2x13+2x150. 85)  2/2(2x13+2x150 . 85) ^ 

Component  1:1^  2(244.59+2680.14)  '  2(244.59+2680.14)  ) 


(  294 . 68  360.17  ^  en 

“  1^  5849.46  '  5849. 46y  “  (0.05038,0.06157). 

(  ;^g.2/2(2xll+2x25.92)  Afx-o.2/2(2xll+2x25 . 92) 

Component...^,:  2(24  4.59+555.48)  '  2(244.59+555.48)  ) 

(  58.01  88.85  ^ 

“  1600.14  '  1600.14  )  ~  (0.03625,0.05553). 

(  ;cl2/2(2xl4+2x23.64)  yf-o.2/2(2xl4+2x23 . 64)^ 

Component  3;  2(244.59+423.19)  '  2(244.59+423.19)  ) 

(  59.8  _9L,^6_^  _  ft  ft^ft,«» 

-  1^  1335.56  '  1335.56  J  “  (0.04477,0.06818). 

(  ;(f^.2/2(2x7+2x89.47)  /f-o.2/2(2x7+2x89. 47) >> 

Component  ,  4 :  2(244.59+1855)  '  2(244.59+1855)  ) 


(  167.35  217.5  ^  ftftft„,.  ft  ft,.,«ft. 

“  ^  4199.18'  4199.18 )  ~  (0.03985,0.05180). 

(  ;c^.2/2(2x11+2x111.19)  ;f!-o.2/2(2xll+2xlll .  19)  \ 

CpmBonent,..5 :  2(244.59+2248.09)  '  2(244.59+2248.09)  J 

(  216.15  272.7  ^  ft.-ft,-  ft  ftir>i-.ftv 

~  4985.36  '  4985.36 )  ~  (0.04336,0.05470). 

The  Bayesian  estimators  of  R(t)  for  t  =  1  hour  are 


Aggregated  Component:  R{t)  =  95.05% 


^(t) 


(  1222.95^  . , 
V 12824.5; 


(  1222.95)  .(  .  1  )  .1 

V  12824.5;  12824.5; 


657.87+56 


Component  1:  R(t)  =  94.55% 


R(t) 


r  ' 

f  244. 
1 2680. 

59  ^ 

.14; 

1  +1 

li 

f  244. 

59  > 

1  +1 

(  1  ^ 

1  +1 

1 2680, 

.14; 

L  2680.14; 

150.85+13 
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Component  2 :  R{t)  =  95.49% 
f 

R{t)  = 


(  2.44_,59 ) 
V  555. 48/ 


\  25.92+11 


(  24, 4. .59^  _ 1 _ ^ 

I  V  555.487  555.487 


+  1 


Component  3:  R{t)  =  94.52% 

R(t)  = 


r  244. 
1  423. 

59> 

197 

+  1 

\ 

f  244. 

59> 

1  + 

r  1  .  > 

i+ij 

1  423. 

I9J 

l 423. I9J 

N  23.64+14 


Component  4 ;  R(t)  =  95.51% 

r  ,  T 

_ V  1855  7 


R{t)  = 


r  .2.44.59.^  (  ...1 . ^ 

V  1855  7  18557 


89.47+7 


Component  5:  R{t)  =  95.22% 

(  (  244.59 ^ 

'  _ V  2248.097 

(  244 


^(t)  = 


N  111.19+11 


(  244.59  ^  ( _ 1 _ ^1  T 

V  2248.097  2248.097 

Substitute  the  appropriate  values  into  R{t)  =  R1R4  +  R2R5 

+  R1R3R5  +  R2R3R4  ~  R1R2R3R4  ~  R1R2R3R5  ~  R1R2R4R5  ~  R1R3R4R5  — 

R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 

Aggregated  Component  to  System:  ^(t)  =  99.49% 

Component  Level  Aggregated  to  System:  R{t)  =  99.52%. 

The  80%  two-sided  probability  intervals  (TBPI)  are 
Aggregated  Component:  (94.81%  <  ■R(l)  <  95.28%) 


exp (-0.05327X1)  ^  R(l)  ^  exp (-0. 04840X1)  . 
Component  1:  (94.03%  <  ^(1)  <  95.09%) 
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exp(-0. 06157X1)  <  ^(1)  <  exp{-0. 05038X1)  . 

Component  2:  (94.60%  <  ^(1)  <  96.44%) 
exp  (-0.05553X1)  <  R{1)  ^  exp  (-0.03625X1)  . 

Component  3:  (93.41%  <  ^(1)  <  95.62%) 
exp(-0. 06818X1)  <  R[l}  <  exp  (-0. 04477X1)  . 

Component  4:  (94.95%  <  ^(1)  <  96.09%) 
exp(-0. 05180X1)  <  ^(1)  <  exp(-0. 03985X1)  . 

Component  5;  (94.68%  ^  ^(1)  <  95.76%) 
exp(-0. 05470X1)  <  R(l)  <  exp (-0. 04336X1)  . 

Substitute  the  appropriate  values  into  J?{t)  =  R1R4  +  R2R5 
+  R1R3R5  +  R2R3R4  —  R1R2R3R4  ~  R1R2R3R5  ~  R1R2R4R5  ~  R1R3R4R5  — 
R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 

Aggregated  Component  to  System;  (99. 44%, 99. 54%) 
Component  Level  Aggregated  to  System;  (99. 38%, 99. 65%) . 
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Appendix  K. 


Simulated  Data  For  A  Bridge  System  With 


Different  Independent  Cozoponents. 

DT  Data  Run  #1 


Time 

Fail 

Time 

Fail 

Time 

Fail 

Time 

Fail 

6.03 

Cl 

134.00 

Cl 

246.06 

C4 

368.17 

C5 

16.66 

C5 

135.49 

Cl 

250.50 

C3 

370.91 

C4 

18.06 

C3 

149.91 

Cl 

253.23 

C3 

378.72 

C4 

18.58 

Cl 

153.05 

Cl 

261.18 

Cl 

382.44 

Cl 

31.40 

C3 

153.22 

C3 

267.25 

C5 

382.74 

C3 

31.68 

C3 

154.03 

C3 

270.00 

Cl 

392.26 

C3 

33.15 

Cl 

154.61 

C2 

273.13 

C4 

394.56 

Cl 

33.23 

C2 

154.88 

Cl 

275.26 

Cl 

404.38 

Cl 

36.88 

C2 

168.43 

C4 

278.04 

C4 

409.59 

C3 

39.22 

C3 

174.94 

Cl 

283.04 

C2 

413.48 

Cl 

46.64 

Cl 

178.08 

Cl 

283.18 

Cl 

418.01 

C2 

51.07 

C5 

185.65 

Cl 

285.32 

C2 

424.22 

C3 

51.28 

Cl 

190.69 

C3 

285.89 

C2 

425.01 

Cl 

56.41 

C5 

198.01 

Cl 

292.50 

C2 

427.08 

C4 

57.23 

C5 

199.16 

Cl 

300.20 

C4 

430.78 

C2 

58.75 

C5 

201.52 

Cl 

309.56 

Cl 

436.65 

C3 

58.89 

C4 

204.01 

Cl 

313.11 

C4 

439.46 

C5 

67.14 

C4 

204.54 

C2 

317.32 

C4 

440.66 

C2 

73.58 

C4 

213.75 

Cl 

319.43 

C3 

447.63 

Cl 

77.80 

C5 

214.40 

C2 

329.52 

Cl 

453.97 

Cl 

90.92 

Cl 

219.63 

C2 

332.56 

Cl 

454.82 

C5 

92.80 

C5 

220.05 

C3 

334.78 

Cl 

456.28 

C3 

105.35 

Cl 

221.16 

Cl 

335.50 

C3 

462.14 

Cl 

107.39 

C3 

222.58 

C2 

339.62 

Cl 

464.38 

Cl 

119.24 

C5 

223.35 

C4 

344.19 

C3 

468.20 

Cl 

119.37 

Cl 

230.07 

C2 

346.51 

Cl 

474.33 

Cl 

125.88 

Cl 

232.06 

Cl 

349.14 

Cl 

475.25 

C5 

127,22 

Cl 

240.86 

Cl 

353.52 

C2 

478.31 

C2 

129.60 

Cl 

241.15 

Cl 

355.06 

C3 

479.17 

Cl 

132.66 

C3 

243.04 

C4 

357.07 

C2 

499.09 

C2 

DT  Data  Run  #2 


Time 

Component 

Failure 

Time 

Component 

Failure 

11.59 

Cl 

171.92 

Cl 

12.80 

Cl 

179.22 

Cl 

17.46 

C2 

181.59 

C2 

27.87 

Cl 

188.38 

Cl 

35.62 

Cl 

190.77 

C3 

39.53 

Cl 

192.12 

Cl 

40.82 

C2 

194.66 

C3 

42.43 

C4 

195.49 

Cl 

43.35 

C5 

196.26 

Cl 

46.11 

C5 

202.67 

C3 

58.12 

C2 

204.10 

Cl 

61.71 

C2 

205.55 

C3 

65.11 

Cl 

206.34 

C4 

68.73 

C4 

210.76  1 

Cl 

71.05 

C3 

211.40 

C5 

77.82 

C3 

212.66 

C2 

83.59 

Cl 

215.04 

Cl 

89.93 

Cl 

218.49 

Cl 

91.82 

C2 

222.49 

C4 

94.53 

C2 

229.88 

C2 

96.31 

C3 

230.52 

C5 

98.13 

C3 

230.85 

C4 

102.47 

Cl 

234.85 

Cl 

106.57 

Cl 

247.80 

C5 

107.42 

C2 

248.72 

C3 

109.43 

Cl 

252.85 

C2 

110.96 

Cl 

253.62 

Cl 

111.38 

Cl 

255.32 

C2 

113.66 

Cl 

259.60 

Cl 

115.82 

Cl 

263.34 

Cl 

116.73 

C5 

265.63 

C3 

127.89 

Cl 

268.53 

C2 

132.20 

C2 

268.90 

Cl 

140.65 

C4 

270.99 

C2 

141.09 

C3 

273.04 

C4 

158.80 

C4 

285.02 

C4 

Time 


285.74 


286.81 


299.44 


305.57 


306.92 


313.42 


326.50 


357.10 


369.98 


374.32 


376.79 


377.61 


379.47 


393.84 


395.29 


396.50 


411.47 


412.58 


412.62 


414.42 


415.28 


423.83 


441.89 


454.39 


458.56 


461.91 


461.98 


462.04 


469.00 


470.35 


473.91 


484.30 


489.34 


490.55 


Component 

Failure 


Cl 


Cl 


C2 


C2 


C2 


C3 


C5 


C3 


C4 


Cl 


C3 


C3 


C3 


C5 


C3 


C2 


C3 


C3 


C2 


C2 


Cl 


Cl 


Cl 


Cl 


C3 


C4 


C4 


C2 


C5 


Cl 


Cl 


C3 


Cl 


C2 
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DT  Data  Run  #3 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

1.41 

Cl 

160.90 

C2 

375.72 

C3 

7.29 

C3 

169.54 

C2 

380.98 

C3 

8.93 

Cl 

180.54 

Cl 

384.56 

C3 

15.20 

C2 

182.31 

Cl 

390.91 

C3 

23.93 

C2 

198.68 

Cl 

391.00 

Cl 

24.83 

C2 

200.50 

Cl 

391.32 

C2 

38.72 

Cl 

201.10 

Cl 

400.48 

Cl 

46.53 

Cl 

220.33 

C4 

402.32 

C2 

50.34 

Cl 

222.58 

Cl 

404.66 

C5 

51.13 

C2 

231.44 

C5 

408.31 

C4 

59.72 

Cl 

235.47 

C5 

409.49 

C2 

65.58 

Cl 

243.10 

C2 

414.94 

Cl 

70.22 

C5 

243.74  I 

Cl 

416.24 

Cl 

73.84 

C4 

249.61 

C2 

416.31 

Cl 

77.92 

C2 

252.36  ! 

C2 

419.10 

C2 

80.99 

C3 

258.72 

C4 

422.68 

C2 

80.99 

C2 

271.75 

Cl 

426.22 

C2 

81.05 

C3 

275.02 

Cl 

426.80 

Cl 

83.69 

C2 

277.05 

C5 

435.44 

C3 

88.17 

Cl 

287.67 

C4 

441.74 

Cl 

90.51 

Cl 

292.83 

C5 

442.35 

Cl 

96.39 

Cl 

295.49 

C3 

445.02 

Cl 

103.22 

Cl 

296.63 

Cl 

451.33 

Cl 

103.60 

C2 

310.30 

C3 

455.83 

Cl 

105.91 

Cl 

315.67 

C2 

457.19 

C2 

112.90 

Cl 

316.78 

C2 

458.45 

Cl 

116.65 

C2 

320.29 

C3 

459.87 

Cl 

126.22 

C2 

324.26 

C3 

461.43 

Cl 

127.22 

Cl 

327.27 

Cl 

463.48 

C5 

144.20 

C3 

332.30 

Cl 

466.94 

C3 

145.90 

C2 

340.06 

C2 

474.84 

C2 

147.47 

C5 

345.05 

C3 

478.73 

C3 

152.86 

Cl 

358.16 

Cl 

488.40 

C4 

159.59 

C3 

360.81 

Cl 

495.40 

C2 

159.82 

C2 

368.93 

C2 

497.97 

Cl 
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DT  Data  Run  #4 


Time 

Component 

Failure 

Time 

Component 

Failure 

0.22 

Cl 

147.31 

C2 

0.67 

C4 

152.25 

C4 

3.65 

Cl 

158.28 

C5 

16.32 

Cl 

167.07 

C5 

19.52 

C2 

169.02 

C3 

30.00 

C5 

174.94 

C2 

38.34 

Cl 

183.38 

Cl 

40.80 

C2 

186.99 

Cl 

46.32 

Cl 

187.80 

C4 

49.32 

C2 

189.02 

C5 

50.08 

C3 

189.94 

C4 

55.09 

Cl 

191.47 

Cl 

64.39 

C4 

199.83 

C2 

65.60 

Cl 

203.21 

Cl 

66.95 

Cl 

210.58 

C3 

67.51 

Cl 

210.59 

C3 

74.37 

C5 

214.07 

C2 

75.71 

Cl 

214.64 

Cl 

78.55 

Cl 

216.86 

Cl 

86.86 

Cl 

219.45 

Cl 

86.87 

Cl 

222.68 

C2 

91.08 

C5 

225.89 

C3 

91.80 

Cl 

236.74 

Cl 

92.30 

C2 

237.05 

C2 

96.83 

C4 

240.41 

Cl 

98.62 

Cl 

240.59 

C4 

100.03 

Cl 

240.99 

Cl 

106.39 

C4 

241.17 

C2 

106.40 

Cl 

241.55 

C4 

107.29 

Cl 

245.77 

C2 

116.25 

C3 

265.82 

Cl 

118.98 

Cl 

268.15 

Cl 

127.47 

C3 

271.20 

C4 

128.22 

Cl 

278.93 

C5 

130.74 

C3 

279.34 

C2 

131.84 

Cl 

288.12 

C5 

132.45 

C3 

292.76 

Cl 

141.67 

C2 

294.07 

C2 

144.92 

C2 

295.04 

Cl 

146.64 

C3 

300.05 

C2 

Time 


309.30 


320.00 


326.05 


330.21 


330.35 


331.48 


331.98 


351.89 


354.28 


361.43 


362.59 


362.72 


366.47 


376.66 


378.74 


379.84 


388.36 


397.50 


402.61 


402.89 


404.92 


423.19 


428.51 


429.30 


431.57 


434.60 


438.24 


440.71 


447.17 


450.31 


464.64 


464.67 


465.08 


465.85 


471.07 


487.62 


491.53 


495.31 


Component 

Failure 


C2 


Cl 


Cl 


C2 


C5 


Cl 


Cl 


C2 


Cl 


C4 


Cl 


C2 


Cl 


Cl 


C2 


Cl 


C3 


C2 


Cl 


C3 


C2 


Cl 


C4 


C3 


C2 


C3 


Cl 


Cl 


C3 


C5 


C3 


C3 


Cl 


Cl 


Cl 


Cl 


C3 


C2 
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DT  Data  Run  #5 


Time 

Component 

Failure 

Time 

Component 

Failure 

Time 

Component 

Failure 

0.67 

C2 

166.70 

C2 

322.14 

C2 

0.82 

C3 

175.46 

C4 

323.43 

C5 

7.46 

C2 

183.11 

C3 

325.20 

Cl 

9.39 

C2 

186.66 

C2 

329.85 

Cl 

10.75 

Cl 

192.37 

C5 

337.16 

C2 

11.20 

C4 

197.14 

Cl 

338.71 

Cl 

11.56 

C3 

197.29 

C3 

339.89 

Cl 

13.02 

C2 

197.48 

C5 

346.68 

C3 

17.47 

C5 

205.10 

C2 

349.52 

C5 

17.70 

Cl 

209.58 

C2 

358.04 

Cl 

30.01 

C5 

210.49 

Cl 

359.73 

Cl 

37.85 

C2 

213.39 

C3 

367.10 

C3 

44.40 

Cl 

223.49 

C4 

367.43 

C5 

48.36 

C4 

232.47 

C2 

373.47 

C5 

53.53 

Cl 

237.11 

Cl 

374.25 

Cl 

53.64 

Cl 

237.44 

Cl 

380.21 

Cl 

61.29 

C5 

239.24 

C3 

381.77 

C2 

61.65 

C4 

240.12 

C2 

385.44 

C4 

62.00 

Cl 

241.92 

C2 

388.47 

C4 

63.82 

C3 

242.26 

Cl 

390.76 

Cl 

66.33 

C2 

244.28 

C4 

391.40 

Cl 

69.25 

C3 

253.14 

C4 

392.34 

C4 

73.82 

Cl 

253.66 

Cl 

392.80 

Cl 

73.87 

C3 

256.45 

C2 

403.73 

C4 

74.55 

Cl 

258.81 

C2 

406.00 

C3 

74.64 

Cl 

262.10 

C3 

408.46 

Cl 

78.55 

C3 

267.99 

C5 

409.87 

C2 

78.99 

C5 

269.02 

C3 

411.74 

Cl 

85.24 

C5 

275.07 

Cl 

422.27 

C2 

106.44 

C5 

278.08 

C4 

1  ...  .... 

437.73 

C4 

108.23 

Cl 

279.08 

C2 

441.40 

Cl 

109.67 

Cl 

280.21 

Cl 

443.16 

Cl 

112.07 

!.  Cl 

283.38 

C4 

444.36 

C4 

114.33 

C3 

290.57 

C2 

462.58 

Cl 

118.03 

Cl 

290.79 

Cl 

466.10 

C2 

122.81 

C3 

294.19 

C4 

468.55 

C3 

149.15 

C2 

295.60 

Cl 

473.70 

C3 

149.39 

Cl 

298.09 

C2 

474.43 

Cl 

149.73 

C3 

301.78 

Cl 

478.64 

C4 

150.20 

C5 

309.33 

C3 

482.65 

C5 

151.99 

Cl 

312.76 

Cl 

492.94 

C4 

157.39 

Cl 

318.88 

C3 

498.76 

Cl 

161.74 

Cl 

321.82 

Cl 
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OT  Data  Runs  #1  and  #2 


I  OT  Data  Run  #1 

OT  Data  Run  #2 

Time 

Component 

Failure 

Time 

Component 

Failure 

6.08 

C3 

22.69 

Cl 

7.82 

Cl 

23.44 

C5 

8.27 

C2 

27.11 

Cl 

9.03 

C2 

30.16 

C3 

14.80 

Cl 

34.43 

Cl 

18.06 

C3 

37.63 

C2 

32.21 

C2 

38.51 

Cl 

33.01 

Cl 

41.76 

C2 

33.05 

C2 

42.56 

C4 

38.23 

C2 

44.49 

Cl 

44.42 

C3 

50.43 

C2 

47.53 

Cl 

53.33 

C4 

47.62 

Cl 

57.57 

C4 

57.23 

Cl 

61.55 

C4 

59.43 

Cl 

67.96 

Cl 

66.62 

C3 

78.81 

Cl 

67.16 

Cl 

89.34  1 

Cl 

76.87 

C2 

93.12 

Cl 

78.42 

Cl 

104.07 

C3 

85.43 

C2 

106.78 

Cl 

96.68 

C2 

116.38 

Cl 

96.96 

Cl 

117.95  ' 

Cl 

98.93 

C5 

119.96 

C3 

109.81 

C2 

113.89 

C3 

114.65 

Cl 

120.78 

C2 

124.09  ^ 

C4 
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Appendix  L.  Calcula'tion  For  A  Bridge  System  Reliability 
With  Different  Independent  Components. 

Consider  a  system  composed  of  5  different  independent 
components  given  in  Figure  1  (page  30) .  The  time  to  failure 
of  each  component  is  assumed  to  be  exponentially 

distributed.  Let  =  0.1,  A,2  =  0.05,  X3  =  0.033,  X4  =  0.025, 

and  I5  =  0.02  be  the  respective  constant  failure  rates 

(failure  per  hour)  of  each  component. 

For  the  true  system,  using  a  mission  time,  t  of  one 
hour,  the  reliability  of  components  1,  2,  3,  4,  and  5  are 
Ri  =  exp (-0.10)  =  90.48% 

R2  =  exp (-0.05)  =  95.12% 

R3  =  exp(-0.033)  =  96.75% 

R4  =  exp (-0.025)  =  97.53% 

Rs  =  exp (-0.02)  =  98.02%. 

The  reliability  of  the  system  can  be  calculated  as 
Rs  =  R1R4  +  R2R5  +  R1R3R5  +  R2R3R4  —  R1R2R3R4  ~  R1R2R3R5  “ 
R1R2R4R5  —  R1R3R4R5  —  R2R3R4R5  t  2R1R2R3R4R5  ~  99.48%. 

Using  a  computer  program  called  RAPTOR,  which  is 

developed  by  AFOTEC,  to  simulate  the  true  system,  the 

following  data  are  obtained  (to  review  the  simulated  data, 
see  Appendix  K) : 
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Table  LI:  Component  1  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

li  =  Xi/Ti 

1 

499.09 

51 

0.10219 

2 

490.55 

41 

0.08358 

3 

497.97 

44 

0.08836 

4 

495.31 

52 

0.10499 

5 

498.76 

48 

0.09624 

Table  L2:  Component  1  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

HSSQSH 

1 

124.09 

11 

0.08865 

2 

119.96 

12 

0.10003 

Table  L3;  Component  2  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

HISEBEESH 

1 

499.09 

19 

0.03807 

2 

490.55 

23 

0.04689 

3 

497.97 

30 

0.06025 

4 

495.31 

26 

0.05249 

5 

498.76 

25 

0.05012 

Table  L4 ;  Component  2  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

124.09 

10 

0.08059 

2 

119.96 

3 

0.02501 

Table  L5:  Component  3  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

499.09 

22 

0.04408 

2 

490.55 

21 

0.04281 

3 

497.97 

17 

0.03414 

4 

495.31 

18 

0.03634 

5 

498.76 

22 

0.04411 

Table  L6:  Component  3  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

1 

124.09 

5 

0.04029 

2 

119.96 

3 

0.02501 

148 


Table  L7 :  Component  4  DT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

499.09 

15 

0.03006 

2 

490.55 

12 

0.02446 

3 

491.91 

6 

0.01205 

4 

495.31 

12 

0.02423 

5 

498.76 

18 

0.03609 

Table  L8 :  Component  4  OT  Data 


Run  i 

Ti  Hours 

Xi  Failures 

Xi  =  Xi/Ti 

1 

124.09 

1 

0.00807 

2 

119.96 

4 

0.03334 

Table  L9:  Component  5  DT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

1 

499.09 

13 

0.02605 

2 

490.55 

9 

0.01835 

3 

497.97 

8 

0.01607 

4 

495.31 

10 

0.02019 

5 

498.76 

15 

0.03008 

Table  LIO:  Com 

ponent  5  OT  Data 

Run  i 

Ti  Hours 

Xi  Failures 

WKBSBEmKm 

1 

124.09 

1 

0.00806 

2 

119.96 

1 

0.00834 

The  Classical  Exponential  Method. 

Assuming  the  time  to  failure  of  each  component  has  the 
exponential  distribution,  point  estimates  for  the  failure 
rate  X  and  80%  confidence  intervals  about  X  will  be 
computed. 


Component 

1  DT: 

^  =  1 

-  1 

f  236  ^ 

=  0.09510 

L  r  J  -  1 

1  2481.687 

Component 

1  OT: 

1  =  1 

-  1 

r  .23  _  \ 

=  0.09424 

ITz)  -  ! 

K  244.057 
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Component  2  DT:  A  =  (  ^3)  =  (7491^)  "  0.04956 

Component  2  OT:  1  =  (-f^)  =  (^44^)  =  0.05327 

Component  3  DT:  ^  =  (fi")  =  (  2481^)  "  0.04030 

Component  3  OT:  1  =  (  244^. Os)  0.03278 

Component  4  DT:  A  =  =  (  248^1^. 68 )  “  0.02539 

Component  4  OT:  A  =  (^)  “  (  244^. 05 )  ^  0.02049 

Component  5  DT:  1  =  “  (  248^1^.  68)  ~  0.02216 

Component  5  OT:  2  =  (l^)  (  244^. 59)  0.00820 

where  fi  =  236  is  the  number  of  failures  and  Ti  =  2481.68  is 

the  total  test  time  from  Table  LI  and  similarly  fa  =  23  and 

Ta  =  244.05  from  Table  L2,  fs  =  123  and  T3  =  2481.68  from 

Table  L3,  f^  =  13  and  T4  =  244.05  from  Table  L4,  fs  =  100  and 

Ts  =  2481.68  from  Table  L5,  f  e  =  8  and  Tg  =  244.05  from  Table 

L6,  fj  =  63  and  T7  =  2481.68  from  Table  L7,  fs  =  5  and  Ts  = 

244.05  from  Table  L8,  fg  =  55  and  Tg  =  2481.68  from  Table 

L9,  fio  =  2  and  Tio  =  244.05  from  Table  LIO. 

The  80%  two-sided  confidence  intervals  (TCI)  are 

Component  1  DT :  (0.08726,0.10311) 

f  /o. 2/2(2x23 6)  2/2(2x236) \  f  433.08  511.78  ^ 

2x2481.68  '  2x2481.68  J  =  4963.36  '  4963. 36y 

Component  1  OT:  (0.07010,0.12014) 
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;t:Q.2/2(2x23)  ;>ri-o. 2/2(2x23)^  r  34.22  58.64  ^ 

2x244.05  f  2x244.05  )  “  ^  488.1  '  488. ly 


Component  2  DT;  (0.04393,0.05537) 


(  ^^272(2x123) 

;f?-o.272(2xl23)>^ 

(  218.04  274.82  ^ 

^  2x2481.68 

t  2x2481.68  ) 

“  4963.36  •  4963.36  ) 

Component  2 

OT:  (0.03543,0. 

07286) 

(  ;f"o.272(2Xl3) 

/l-0. 272(2x13)') 

(  17.29  35.56^ 

2x244.05  • 

2x244.05  )  '■ 

"  488.1  '  488.1 ) 

Component  3 

DT:  (0.03523,0. 

.04554) 

(  ;f§.272(2Xl00) 

/i-O.272(2X100) 

(  174.84  226.02  ^ 

^  2X2481.68 

'  2x2481.68  ) 

~  V  4963.36  '  4963.36 ) 

Component  3 

OT:  (0.01908,0. 

04823) 

(  ^0.272(2x8) 

/i-o. 272(2x8)') 

(  9.31  23.54  ^ 

V  2x244.05  ' 

2x244.05  )  =  ' 

^  488.1  '  488.1/ 

Component  4 

DT:  (0.02138,0. 

.02956) 

(  yE.272(2x63) 

jfi-o. 272(2x63) ') 

(  106.13  146.72 

^  2x2481.68  ' 

2x2481.68  )  = 

"  V  4963.36  •  4963.36 ) 

Component  4 

OT:  (0.00997,0. 

03275) 

(  ^0.272(2x5) 

/i-o. 272(2x5)^ 

(  4.87  15.99^ 

^  2x244.05  ' 

2x244.05  )  =  ' 

^  488.1  '  488.1 ) 

Component  5 

DT:  (0.01843,0. 

02607) 

(  ^0.272(2x55) 

^1-0.272(2x55) ') 

(  91.47  129.39 

^  2x2481.68  / 

2x2481.68  )  = 

=  4963.36  '  4963.36 ) 

Component  5 

OT:  (0.00218,0. 

.01594) 

y 0.2 72(2x2)  /i-o. 272(2x2)'^  f  1.06  7.78  ^ 

2x244.05  '  2x244.05  y  I  488.1  '  488.1 )  • 
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The  reliability  estimators  of  R(t)  for  t  =  1  hour  are 


Component 

1 

DT: 

^(t) 

= 

exp(-0. 09510X1) 

90.93% 

Component 

1 

OT: 

^(t) 

exp(-0. 09424X1) 

_ 

91.01% 

Component 

DT: 

^(£) 

exp(-0. 04956X1) 

_ 

95.17% 

Component 

OT: 

^(t) 

exp(-0. 05327X1) 

— 

94.81% 

Component 

3 

DT: 

^(£) 

exp (-0.04030X1) 

96.05% 

Component 

3 

OT: 

^(t) 

exp(-0. 03278X1) 

. 

96.78% 

Component 

DT: 

^(t) 

exp(-0. 02539X1) 

97.49% 

Component 

_4_ 

OT: 

^  ( t) 

exp(-0. 02049X1) 

_ 

97.97% 

Component 

DT: 

S(t) 

exp(-0. 02216X1) 

97.81% 

Component 

5 

OT: 

^(t) 

exp (-0 . 00820X1 ) 

99.18% 

Substitute  the  appropriate  values  into  -R(t)  =  R1R4  +  R2R5 
+  RiRaRs  +  R2R3R4  “  RiR2R3R4  ~  R1R2R3R5  “  R1R2R4R5  ”  R1R3R4R5  ~ 
R2R3R4R5  +  2R1R2R3R4R5  to  Obtained  the  following: 

Component  Level  Aggregated  to  System  DT:R(t)  =  99.50% 
Component  Level  Aggregated  to  System  0T:-R(t)  =  99.51%. 
The  80%  two-sided  confidence  intervals  (TCI)  are 
Component  1  DT:  (90.20%, 91.64%) 

exp(-0. 10311X1)  <  R(l)  <  exp(-0. 08726X1) 

Component  1  OT:  (88. 68%, 93.23%) 

exp(-0. 12014X1)  <  R{1)  ^  exp  (-0. 07010X1) 

Component  2  DT:  (94. 61%, 95. 70%) 
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exp(-0. 05537X1)  <  ^(1)  <  exp  (-0.04393X1) 

Component  2  OT;  (92 .97%, 96. 52%) 

exp(-0. 07286X1)  <  R(l)  <  exp (-0 . 03543X1 ) 

Component  3  DT:  (95. 55%, 96. 54%) 

exp(-0. 04554X1)  <  R(l)  <  exp  (-0. 03523X1) 

Component  3  OT;  (95. 29%, 98. 11%) 
exp(-0. 04823X1)  <  ^(1)  <  exp (-0. 01908X1) 

Component  4  DT :  (97. 09%, 97. 89%) 

exp(-0. 02956X1)  <  ^(1)  <  exp(-0. 02138X1) 

Component  4  OT;  (96. 78%, 99. 01%) 

exp (-0.03275X1)  <  ^(1)  ^  exp (-0 . 00997X1 ) 

Component  5  DT;  (97. 43%, 98. 17%) 

exp(-0. 02607X1)  <  ^(1)  <  exp  (-0.01843X1) 

Component  5  OT;  (98.42%, 99.78%) 

exp(-0. 01594X1)  ^  ^(1)  <  exp  (-0 . 00218X1 ) 

Substitute  the  appropriate  values  into  R{,t)  =  R1R4  +  R2R5 
+  R1R3R5  +  R2R3R4  —  R1R2R3R4  ~  R1R2R3R5  —  R1R2R4R5  ~  R1R3R4R5  — 
R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 


(99. 14%, 99. 76%) . 
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The  Bayesian  Exponential  Method. 

Bayesian  Exponential  Prior. 

For  Component  1  data,  from  Table  LI: 


E{l)  X  =  0.09507  and  V[X)  =  =  0.00008. 


These  yielded  the  parameter  estimates 


a 


(  0.09507^ 
^  0.00008 


=  112.98  and 


(  0.09507^ 
V  0.000087 


=  1188.38. 


For  Component  2  data,  from  Table  L3; 


E{X)  =  A  =  0.04956  and  V(A)  =  Sx^  =  0.00007. 


These  yielded  the  parameter  estimates 


a 


(  0.04956^^ 
V  0.00007  ) 


=  35.09  and 


(  0.04956^ 
V  0.000077 


708. 


For  Component  3  data,  from  Table  L5: 


E[X)  ^  k  =  0.04030  and  V{X)  =  sx^  =  0.00002. 


These  yielded  the  parameter  estimates 


a 


0.04030^ 

0.00002 


=  81.21  and 


f  0.04030^ 
V  0.00002>' 


=  2015. 


For  Component  4  data,  from  Table  L7: 

E(X)  =  T  =  0.02538  and  V(X)  =  Sx"  =  0.00008. 


These  yielded  the  parameter  estimates 


f  0.02538^^  -  „ 

“  =  I  0.00008  7  “  P 


f  0.02538^ 
V  0.000087 


317.25. 


For  Component  5  data,  from  Table  L9: 

E(X)  =  T  =  0.02215  and  V(X)  =  sx^  =  0.00003. 
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These  yielded  the  parameter  estimates 


r  0.02215^^ 
I  0.00003  J 


16.35  and  P  - 


738.33. 


Bayesian  Exponential  Posterior  Analysis. 


Component 

T  - 

1  =  1 

(  23+112.98  ^ 

0.09493. 

1 : 

^  244. 05  +  1188. 3U 

Component 

1  =  1 

(  13+35.09  ^ 

0.05051. 

Z : 

^  244.05  +  708/ 

Component 

2  =  1 

(  8+81.21  ^ 

0.03949. 

^  244.05  +  2015/ 

Component 

J.: 

A  =  1 

(  5+8.05  ^ 

0.02325. 

^  244.05  +  317.25/ 

Component 

_5 : 

A  =  1 

(  2+16.35  ^ 

0.01868. 

<  244.05  +  738.33/ 

The  80%  two-sided  probability  intervals  (TBPI)  are 

f  2/2(2x23+2x112 . 98)  2/2(2x23+2x112 . 98)  ^ 

Component  l;l^  2(244.05+1188.38)  '  2(244.05+1188.38)  J 

(  241.63  _301.;J1_'\  «  ../vcm 

-  2864.86  '  2864.86 )  ~  (0.08434,0.10515). 

C  Ifo . 2/2(2x1  3+2x3 5  •  0 9)  j|^^_Q^2/2(2xl3+2x35*09)^ 

Component  2:  ’  2(244.05+708)  '  2(244.05+708)  ) 

(  78.73  114. 13^ 

~  y  1904.1  '  1904.1  J  ~  (0.04135,0.05994) . 

(  ;if 0.2/2 (2x8+2x81. 21)  2/2(2x8+2x81. 21)  A 

Component  3:  2(244.05+2015)  '  2(244.05+2015)  ) 

(  154.29  202.57  ^ 

~  1^4518.1'  4518. ly  ~  (0.03415,0.04484). 

f  lfo.2/2(2x5+2x8 . 05)  Xi-o , 2/2(2x5+2x8 . 05)  'N 
Component  4:  (^  2(244.05+317.25)'  2(244.05+317.25)  J 

(  17.29  _35.i_^^  .. 

“  1^1122.6'  1122. 6y  “  (0.01540,0.03168). 
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f  ;Co.2/2(2x2+2x16.35)  ;ftQ.2/2(2x2+2xl6. 35) 

Component  5:  ^  2(244.05+738.33)  '  2(244.05+738.33) 


f  25.64 
“  1964.76  ' 


47.21 

1964.76 


=  (0.01305,0.02403). 


The  Bayesian  estimators  of  R(t)  for  t  =  1  hour  are 


;:omponent  1 :  R{t)  =  90.95% 


f  f  244.05  ^  .  T 

_  _ V  1188.38;  _ 

-  f  244.05  )  f _ 1 _ 

I  V  1188. 38^  1188. 38J 


112.98+23 


Component  2;  R{t)  =  95.08% 

f  (  244.^^  .  T 

_  V . 708. .J.:!::.. 

~  f  244.Q5S  V  1  ^ 

vv  708  y'^woe; 


35.09+13 


Component  3:  R(t)  =  96.13% 


R(t)  = 


(  244.05^  .  T 
J^2015jJjL 


r  244,05.^  +r_i_)  +1 
I  2015  )  2015^ 


Component  4 :  R(t)  =  97.70% 


R(t)  = 


f  244.05^  .  T 
V  317.25/ 

244.05^  ,f  1  ' 

317. 25^  317.25, 


Component  5:  !?(£)  =  98.15% 


R{t)  = 


r  244.05^  .  T 
^  138. 33J 

f  244.05S  _ 1 _ ^ 

V  738.33;  738.33; 


Substitute  the  appropriate  values  into  R{t)  =  R1R4  +  R2R5 
+  R1R3RS  +  R2R3R4  ~  R1R2R3R4  ~  R1R2R3R5  “  R1R2R4R5  ~  R1R3R4R5  ~ 


R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 
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Component  Level  Aggregated  to  System:  R{t)  =  99.50%. 

The  80%  two-sided  probability  intervals  (TBPI)  are 

Component  1 :  (90.02%  <  R{1)  ^  91.91%) 
exp(-0. 10515X1)  <  ^(1)  <  exp(-0. 08434X1)  . 

Component  2 :  (94.18%  <  -^(1)  ^  95.95%) 
exp(-0. 05994X1)  <  ^(1)  <  exp(-0. 04135X1)  . 

Component  3 ;  (95.62%  <  R{1)  <  96.64%) 
exp{-0. 04484X1)  <  ^(1)  <  exp (-0 . 03415X1 )  . 

Component  4;  (96.88%  ^  -^(1)  <  98.47%) 

exp(-0. 03168X1)  <  ^Rd)  <  exp  (-0. 01540X1)  . 

Component  5 :  (97.63%  <  ■R(l)  <  98.70%) 
exp(-0. 02403X1)  <  R(l)  <  exp  (-0 . 01315X1 )  . 

Substitute  the  appropriate  values  into  R{t)  =  R1R4  +  R2R5 
+  R1R3R5  +  R2R3R4  —  R1R2R3R4  ~  R1R2R3R5  —  R1R2R4R5  —  R1R3R4R5  — 
R2R3R4R5  +  2R1R2R3R4R5  to  obtained  the  following: 

Component  Level  Aggregated  to  System:  (99. 33%, 99. 65%) . 
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